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SUMMARY 


Nonequilibrium  effects  in  supersonic-nozzle  flows  are  important  in 
various  technical  applications  such  as  propulsion,  hypersonic  wind- 
tunnel  testing,  and  simulation  of  streamtube  flows  occurring  about  a  body 
in  flight.  This  report  discusses  such  effects  with  emphasis  on  survey¬ 
ing  the  research  of  recent  years.  The  scope  is  limited  to  homogeneous 
gas-phase  nozzle  flows  only,  and  to  departures  from  therirochemical 
equilibrium  arising  from  collisional  relaxation  of  internal  degrees  of 
molecular  excitation  and  from  chemical  reaction,  including  ionization. 
Translational-rotational  equilibrium  is  assumed.  Nonequilibrium  pheno¬ 
mena  due  to  classical  viscous  effects,  condensation,  rarefied  gas- 
dynamic  effects,  radiation,  and  translational  energy  difference  between 
free  electrons  and  heavy  particles  are  not  discussed. 

The  first  half  of  the  report  provides  a  background  on  the  high- 
temperature  phenomena  considered,  and  covers  equilibrium  properties  of 
gases  and  the  basic  nonequilibrium  rate  processes  involving  vibration 
and  electronic  excitation,  chemical  reactions  of  neutral  species,  and 
ionization.  The  last  half  of  the  report  discusses  the  nonequilibrium 
effects  arising  from  these  rate  processes  in  nozzle  flow  expansions. 
Bnphasis  is  on  diatomic  gases  and  air,  but  complex  mixtures  for  pro¬ 
pulsion  applications  are  also  discussed.  Bibliographies  are  included 
on  related  nonequi librium  flows,  as  well  as  on  nozzle  flows, 

The  review  indicates  great  progress  over  the  past  decade  in  increased 
knowledge  and  understanding  of  nonequilibrium  nozzle  flows  of  technical 
interest.  An  important  contribution  has  been  made  by  numerical  studies 
using  high  speed  machine  computers.  However,  much  remains  to  be  done  as 
regards  more  certain  and  complete  mathematical  description  of  various 
physical-chemical  rate  processes  which  are  important.  In  general,  much 
more  emphasis  is  now  needed  on  basic  studies,  particularly  on  experi¬ 
mentation,  if  a  significant  rate  of  progress  is  to  be  maintained. 
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RESUME 


Les  effets  de  non  dqulllbre  dans  les  dcouleaents  de  tuyeres  super- 
sonlques  sont  laportants  dans  dlffdrentes  applications  techniques, 
telles  que  la  propulsion,  les  essals  en  soufflerle  hypersonlque,  et  la 
slaulatlon  des  dcouleaents  en  tubes  de  courant  qul  se  produlsent  autour 
d'un  corps  en  vol.  Le  present  rapport  dtudle  ces  effets  en  tenant 
ccapte  en  partlculler  des  recherches  effectudes  au  cours  de  ces  dernldres 
amides  II  se  borne  4  n'  exaalner  que  les  dcouleaents  de  tuydre  hoaogdnes 
de  phase  gazeuse,  et  les  cas  de  non  dqulllbre  therao-chlaique  dfl  4  la 
detente  de  collision  de  degrds  Internes  de  1*  excitation  aoldculalre  et 
aux  reactions  chlaiques,  y  coaprls  1’  ionisation.  L’ dqulllbre  de  transla¬ 
tion  et  de  rotation  est  suppose.  Les  phdnoadnes  de  non  dqulllbre  ay  ant 
pour  orlglne  les  effets  vlsqueux  classlques,  la  condensation,  les  effets 
gaso-dynaalques  rardflds,  le  rayonneaent  et  la  difference  d’ dnergie  de 
translation  entre  les  electrons  libres  et  les  partlcules  lourdes  ne  sont 
pas  traltds. 

La  prealdre  partle  du  rapport,  qul  sert  de  fond  pour  les  phdnoadnes 
dtudlds,  porte  sur  les  proprldtds  en  dqulllbre  des  gaz  et  les  procddds 
fondaaentaux  de  vitesse  des  reactions  en  non  dqulllbre  lap 11 quant 
1' excitation  due  4  des  vibrations  et  electronlque,  des  reactions  chlaiques 
des  espdces  neutres.  et  1’  ionisation.  La  deuxldae  partle  dtudle  les 
effets  de  non  dqulllbre  dus  i  ces  procddds  de  vitesse  dans  dea  dilatations 
d’  dcouleaents  de  tuydre.  On  exaalne  surtout  les  gaz  diatoaiques  et  l’alr, 
aais  dtudle  dgaleaent  des  adlanges  coaplexes  pour  des  applications  4  la 
propulsion.  L’ Agardographie  coaporte  des  rdfdrences  blbllographlques  sur 
des  dcouleaents  en  non  dqulllbre  correspondants,  alnsl  que  sur  les 
dcouleaents  de  tuydre 

La  prdsente  dtude  lndique  les  progrds  laportants  rdallsds  au  cours  de 
ces  demldres  dlx  anndes,  »t  qul  se  traduisent  par  des  connalssances  et 
une  coaprdhenslon  aadllordes  des  dcouleaents  de  tuydre  dans  des  conditions 
de  non  dqulllbre  prdsentant  un  interdt  technique.  Des  dtudes  nuadrlques 
effectudes  d  1’ aide  de  calculateurs  4  grande  vitesse  ont  perm  Is  d’y 
apporter  une  contribution  laportante.  Ndanaolns  11  reste  beaucuup  4 
falre  en  ce  qul  concerne  une  description  plus  certalne  et  plus  coapldte 
de  dlffdrents  procddds  de  vitesse  physlco-chiaiques  qul  sont  lwortants. 

D*  une  faqon  gdndrale,  11  faudra  aettre  beaucoup  plus  d*  accent  sur  des 
dtudes  de  base,  en  partlculler  sur  des  expdriaentations,  si  l’on  veut 
qu’  une  cadence  d' avanceaent  lnportante  solt  aalntenue. 
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NOTATION 


a  sound  speed;  also  tan  f3,  where  P  =  asymptotic  half  angle  of  nozzle 

A  parameter  in  inter-particle  potential;  also  cross-section  area  of  nozzle 

C'  constant  in  rate  equation  for  Lighthill  gas 

D  dissociation  energy  per  mole 

Dq  dissociation  energy  per  molecule,  D/N0 

ey  vibrational  energy  per  unit  mass  of  mixture 

e j  energy  of  level  j 

h  enthalpy  per  unit  mass  of  mixture;  also  Planck  constant 

k  specific  reaction  rate  coefficient;  also  Boltanann  constant 

K  equilibrium  constant 

1  =  L/a  ,  where  a  =  tan  0 

L  radius  of  nozzle  throat 

m  mass  of  an  atom  or  molecule;  also  mass  flow  through  nozzle 

M  Mach  number;  also  any  species  of  mixture 

n.N  number  of  particles  per  unit  volume 

n  q 

Nfl  number  of  particles  per  mole  (6.03  x  10“"  / 

p  gas  pressure 

R,  Rg  intra-nuclear  coordinate 

RQ  universal  gas  constant 

Rex  free-gtream  unit  Reynolds  number 

S  entropy  per  unit  mass  of  mixture 

t  time 

T  temperature 

u  gas  velocity 


vi 


U,,,  test-flow  velocity  in  hypersonic  wind  tunnel 
V  intermolecular  potential 

x  distance  along  nozzle  from  throat 

a  mass  fraction  of  gas  in  dissociated  form 

y  mass  concentration  of  species,  raoles/gram  of  mixture;  also  vibrational  non¬ 

equilibrium  parameter 

p  gas  density,  mass  per  unit  volume  of  mixture 

v  characteristic  frequency  of  molecule 

r  relaxation  time 

4  molecular  weight 

6  half  angle  of  conical  portion  of  nozzle 

Subscripts 

i  denotes  reaction  i 

J  denotes  species  j 

o  denotes  reservoir  or  stagnation  condition 

EQ  denotes  equilibrium  state 

v  denotes  vibrational  mode 

F. R  denotes  forward  or  reverse  specific  rate  coefficients 
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SECTION  1 


INTRODUCTION 


Over  the  past  decade  a  great  amount  of  research  has  been  devoted  to  high-temperature, 
physical -chemical  phenomena  of  interest  in  technical  gasdynamics.  Much  of  this  effort 
has  centered  on  nonequilibrium  phenomena  occurring  in  shock  waves.  However,  consider¬ 
able  attention  has  also  been  devoted  to  another  elementary  flow  of  equal  technical 
importance,  namely,  expansion  flow  as  typified  by  flow  in  a  supersonic  nozzle.  The 
purpose  of  this  report  is  to  discuss  nonequilibrium  phenomena  in  such  expansion  flows 
with  emphasis  on  surveying  the  research  results  of  recent  years. 

Nonequilibrium  effects  in  nozzle  flows  have  long  been  of  interest  in  the  propulsion 
field  because  of  the  thrust  loss  resulting  from  chemical  freezing,  or  recombination 
lag,  occurring  in  the  nozzle  expansion  process.  For  example,  some  of  the  earliest 
theoretical  studies  of  chemical  nonequilibrium  effects  in  nozzle  flows  were  done 
around  1947  by  Penner,  Altman,  and  co-workers1,2  in  connection  with  rocket  nozzle 
performance.  In  recent  years  the  problem  of  chemical  freezing  in  propulsive  nozzles 
has  assumed  greater  importance  with  the  development  of  higher-performance  chemical 
rockets,  such  as  hydrogen-fluorine  system,  and  the  interest  in  development  of  advanced 
hypersonic  ram  jet  systems3'4. 

Nozzle-flow  nonequilibrium  is  also  of  much  concern  in  connection  with  hypersonic 
wind-tunnel  testing.  Current  test  devices  such  as  arc  jets  and  shock  tunnels  utilize 
nozzle  expansions  of  air  from  initial  high-temperature  conditions  where  the  air  may  be 
highly  dissociated.  Usually  some  degree  of  freezing  occurs  in  the  nozzle  expansion, 
which  produces  a  test  airflow  not  in  an  equilibrium  state5'6.  Most  often  this  is 
undesirable  as  it  complicates  interpretation  of  test  data  and  can  prevent  proper 
simulation. 

Still  a  further  significance  of  nonequilibrium  nozzle  flows  is  that  the  thermo- 
gasdynamic  environment  is  representative  of  that  occurring  along  certain  streamtubes 
undergoing  expansion  about  a  body  in  hypersonic  flight.  Thus  nozzle-flow  experiments 
provide  a  means  of  simulating  such  streamtube  flows  directly  in  order  to  study  non¬ 
equilibrium  phenomena  in  a  gasdynamic  environment  appropriate  to  the  flight  situation. 
The  importance  of  the  gasdynamic  environment  in  this  respect  is  emphasized  by  recent 
research  results  on  vibrational  relaxation  (discussed  in  Section  4.3)  which  indicate 
important  basic  differences  in  vibrational  behavior  between  shock -wave  and  expansion- 
flow  environments. 

Even  the  limitation  of  subject  to  nonequilibrium  effects  in  nozzle  expansion  flows 
still  leaves  a  very  broad  range  of  specific  phenomena  which  could  be  discussed.  The 
present  scope  is  arbitrarily  limited  to  consideration  of  homogeneous  gas-phase  flows 
only,  and  to  departures  from  the  thermochemical  equilibrium  arising  from  the  collisional 
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relaxation  of  internal  degrees  of  molecular  excitation  and  from  chemical  reaction, 
including  ionization.  However,  translational  and  rotational  degrees  of  freedom  are 
assumed  to  be  in  local  (Maxwell -Boltzmann)  equilibrium.  Thus  only  inviscid  nozzle 
flows  in  the  gas-phase  continuum  regime  are  considered,  and  nonequilibrium  phenomena 
associated  with  classical  viscous  effects,  condensation,  rarefied  gasdynamic  effects, 
radiation,  and  energy  balance  between  free  electrons  and  heavy  particles  are  not 
discussed. 

The  organization  of  the  material  covered  is  as  follows.  Section  2  reviews  the 
methods  available  for  calculating  the  equilibrium  properties  of  gases  at  high  tempera¬ 
tures.  Section  3  discusses  the  basic  rate  processes  of  interest,  which  concern  vibra¬ 
tion  and  electronic  excitation,  chemical  reactions  and  neutral  species,  and  ionization. 
Finally,  Section  4  considers  nonequilibrium  effects  arising  from  such  rate  processes 
in  nozzle  flow  expansions.  While  the  emphasis  is  placed  on  flows  of  diatomic  gases  and 
air,  nozzle  flows  of  complex  mixtures  of  interest  in  propulsion  applications  are  also 
discussed. 

For  each  section,  the  pertinent  auxiliary  material  of  tables,  figures,  and  a  combined 
bib1 iegraphy  and  reference  list  are  located  at  the  end  of  the  section.  The  combined 
bibliography  and  reference  lists  are  organized  under  several  appropriate  subheadings  as 
a  convenience  for  bibliographic  usage.  Although  plane  shock-wave  flows,  and  two-  or 
three-dimensional  inviscid  flows  are  generally  not  discussed,  bibliographies  covering 
nonequilibrium  effects  in  these  two  classes  of  flows  are  included  as  Appendices  A  and  B. 
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SECTION  2 

REVIEW  OF  THERMODYNAMICS  AND  STATISTICAL 
MECHANICS  OF  GASES 


There  are  several  separate  aspects  of  high-temperature  gas-flow  problems  that  require 
the  accurate  description  of  the  thermodynamic  state  of  the  gas  and  which  require  a 
detailed  knowledge  of  its  molecular  and  atomic  structure.  Foremost  is  the  proper 
description  of  the  internal  energy  of  the  gas  under  these  conditions  where  the  assump¬ 
tion  of  a  constant  specific  heat  can  no  longer  be  used.  This  assumption  begins  to  fail 
appreciably  for  air  at  about  1000°K,  when  vibrational  motion  of  the  oxygen  molecule 
becomes  important.  At  higher  temperatures,  when  molecular  dissociation  and  ionization 
become  involved,  these  chemical  reactions  supply  the  primary  contributions  to  the 
internal  energy  (Figure  2-1).  A  similar  motivation  is  the  effect  the  chemical  reactions 

have  on  the  equation  of  state  of  the  gas  simply  because  of  the  change  in  the  total 

number  of  particles  (i.e.  the  number  of  moles  of  gas).  Of  further  importance  is  the 

increase  of  chemical  activity  of  the  air  due  to  the  atomic  gases  that  are  present,  and 

the  increase  in  electrical  conductivity  from  the  free  electrons  that  are  produced.  At 
very  high  temperatures  the  transfer  of  energy  by  radiative  processes  can  be  important, 
and  can  be  calculated  only  with  detailed  knowledge  of  the  population  of  the  many  excited 
states  of  the  atomic  and  molecular  constituents.  For  all  of  these  reasons  it  was  early 
realized  that  hypersonic  flight  would  require  a  thorough  understanding  of  the  physical 
structure  of  the  air  molecules  and  the  chemistry  of  their  behavior  at  high  temperatures. 
The  resulting  investigations  have  provided  tables,  graphs  and  computer  programs  to 
describe  these  equilibrium,  high-temperature  properties  of  air.  The  use  of  the  shock 
tube  to  study  the  chemistry  of  other  gases  and  the  importance  of  ablative  chemistry 
during  reentry  has  motivated  the  calculation  of  high-temperature  properties  for  a 
variety  of  gas  mixtures. 

The  development  of  a  consistent  description  of  the  macroscopic  properties  of  gases 
from  the  microscopic  view  of  atomic  and  molecular  collisions  is  given  in  texts1’2  on 
statistical  mechanics.  Furthermore,  this  development  is  presented  within  the  context 
of  the  application  of  the  final  formulas  to  gasdynamic  problems  in  several  texts3'6 
and  reports7*8.  Within  the  present  section  it  will  not  be  attempted  to  repeat  these 
derivations,  but  instead  the  concepts  will  be  summarized  very  briefly  in  presenting 
the  basic  formulas.  A  more  detailed  discussion  is  then  presented  of  the  real-gas 
properties  that  have  been  obtained  from  these  relations  for  various  gas  mixtures  -  in 
particular  for  air  -  and  the  magnitude  and  importance  of  various  real -gas  effects. 

2.1  Equilibrium  Equations  for  a  Reacting  Gas 

For  a  given  gaseous  system,  the  equilibrium  number  of  particles  in  each  energy  state 
(of  molecules,  atoms  and  ions)  is  determined  by  the  number  of  stationary  states  each 
species  possesses,  the  energy  required  to  occupy  each  state,  and  the  energy  available 
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to  the  system.  As  long  as  the  density  of  the  gas  is  not  so  great  that  intermolecular 
forces  play  a  significant  role,  these  properties  are  sufficient  for  the  equilibrium 
description  and  are  summarized  for  each  species  by  the  partition  function1 


27TmjkT\3/2 
*2  ) 


EJn/kT 


(2-1) 


where  the  quantities  £jn  are  the  energies  of  the  stationary  states  of  species  j  , 
corresponding  to  various  internal  potential  and  kinetic  energies.  The  sum  is  taken 
over  all  energy  levels  (n),  and  the  terms  gjn  (the  degeneracy)  give  the  number  of 
states  with  the  same  energy  e._  .  Each  term  is  weighted  according  to  the  factor 

/k  T  J  •* 

e  Jn  ,  proportional  to  the  equilibrium  population  in  that  state  when  the  energy 
available  to  the  system  is  described  by  a  temperature  T  .  Then,  if  the  energy  levels 
and  degeneracies  are  known  for  a  given  species,  the  partition  function  for  that  species 
can  be  calculated  as  a  function  of  temperature.  With  the  molecular  mas  mj  in  grams, 
Boltzmann’ s  constant  k  in  ergs/°K  and  Planck’ s  constant  h  in  erg  secs,  Q  has 
units  of  cm'3. 


The  fact  that  the  energy  £jn  of  any  state  can  be  written  as  a  sum  of  rotational, 
vibrational  and  electronic  energy  permits  the  partition  for  any  species  to  be  written 
as  a  product  of  the  separate  functions  for  each  of  these  degrees  of  freedom.  The 
rotational  partition  function  is  given  accurately  as 


“rot  =  ST  •  (2'2> 

n 

where  6n  is  the  characteristic  temperature  for  rotation  for  the  molecule,  equal  to 
the  rotational  energy  constant  B  divided  by  the  Boltzmann  constant,  it  is  a  symmetry 
factor  that  is  equal  to  2  for  a  homonuclear  diatomic  molecule,  1  for  a  hetronuclear 
molecule.  The  vibrational  partition  function  is  given,  in  the  simple  harmonic  oscilla¬ 
tor  approximation,  as 

1 

^vib.11  -6  /T  -  <2'3> 

1  -  e 


where  6y  is  the  characteristic  temperature  for  vibration  for  the  molecule.  More 
accurate  expressions  for  the  vibrational -rotational  partition  function,  taking  into 
account  the  non-constant  energy-level  spacing  of  the  vibrator  and  the  coupling  between 
rotation  and  vibration,  are  used  in  the  tables  that  have  been  prepared  for  various 
gases.  No  approximate  formulas  are  available  for  the  electronic  partition  functions, 
so  that  these  sums  must  be  performed  numerically. 

For  any  chemical  reaction  that  can  occur,  the  species  which  are  involved  will  equili¬ 
brate  in  a  fixed  proportion  determined  by  their  partition  function  and  the  energy  of  the 
chemical  reaction.  Thus,  using  the  notation  of  Reference  3,  a  chemical  reaction 
(reaction  i)  is  specified  by  the  equation 
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where  there  are  s  species  and  r  reactions.  The  stoichiometric  coefficients  v{. 

th  1 J 

are  the  number  of  molecules  of  Mj  on  the  lefthand  side  of  the  i  reaction  and 
the  number  on  the  righthand  side.  Then  when  equilibrium  is  attained  at  a  given 
density  p  ,  the  species  concentration  (7j  in  moles/gm)  is  given  hy  a  set  of  coupled 
algebraic  equations 

\yj  =  p  J‘1  K±  ;  i  =  1  .  2,  ...  n  (2-5) 

.  =i 


and  by  equations  which  describe  the  conservation  of  the  nuclei  of  the  various  species 
and  the  conservation  of  electrical  charge.  The  quantities  are  the  equilibrium 
constants,  and  can  be  calculated  for  any  temperature  from  a  knowledge  of  the  appropriate 
partition  functions  through  the  relation 
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(R0T) 


e  1 


(2-6) 


In  this  equation  RQ  is  the  universal  gas  constant  and  N0  is  Avagadro’ s  number. 
The  number  of  species  which  must  be  included  in  any  particular  calculation  will  depend 
on  the  temperature  range  involved.  For  temperatures  beyond  about  10,000°K,  multiply- 
ionized  species  can  be  present  in  appreciable  quantities,  and  the  number  of  equations 
used  must  be  increased  accordingly.  Whatever  the  order  of  complexity  that  is  required, 
the  solution  of  this  set  of  simultaneous  equations  permits  calculation  of  the  chemical 
composition  of  the  high-temperature  gas.  The  thermodynamic  properties  of  the  gas  can 
then  be  determined  by  combining  the  expressions  for  the  properties  of  the  separate 
species  in  terms  of  their  partition  functions,  weighted  according  to  their  abundance. 
The  main  expressions  of  interest  for  gas  dynamic  calculations  are  the  equation  of 
state  and  the  expression  for  the  enthalpy  of  the  gas,  which  are  as  follows: 
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The  compressibility  factor  z  is  equal  to  the  ratio  of  the  total  number  of  particles 
to  the  original  number.  For  the  derivation  of  these  and  other  thermodynamic  expressions, 
reference  is  again  made  to  texts,  such  as  References  1  and  4. 


2.2  Numerical  Equilibrium  Calculations  for  Gases,  Especially  Air 

Equilibrium  equations  as  described  in  the  previous  sections  have  been  solved  for  a 
large  number  of  gas  mixtures,  and  results  have  been  published  in  tabular,  graphic,  and 
curve-fit  forms.  The  strong  interest  in  air  in  the  temperature  range  associated  with 
ICBM  reentry  led  to  the  publication  of  several  sets  of  tables  for  oxygen,  nitrogen  and 
air.  In  these,  equilibrium  constants  and  gas  compositions  have  been  tabulated  or 


graphed  for  various  air-species  reactions,  using  these  ideal-gas  assumptions  for  the 
component  species.  The  results  differ  slightly  depending  on  the  method  used  in  approxi¬ 
mating  the  vibrational  states,  on  the  exact  values  used  for  the  energy  levels,  and  on 
the  number  of  energy  states  used  in  the  calculations  of  the  partition  functions.  The 
effects  of  vibration -rotation  corrections  and  effects  of  the  use  of  slightly  different 
spectroscopic  constants  are  compared  by  Allison26  for  air  and  for  C02-N2-Ar  mixtures. 

The  most  important  experimental  parameters  that  are  required  are  the  heats  of  formation 
Dj  ,  the  dissociation  and  ionization  energies  of  the  component  species.  These  values 
have  been  accurately  measured  for  the  components  of  high -temperature  air  and  many  other 


The  question  of  the  total  number  of  levels  to  use  in  Equation  (2-1)  involves  the 
difficulty  that  there  are  an  infinite  number  of  bound  states  near  the  ionization  energy, 
so  that  the  sum  in  fact  diverges.  Different  criteria  for  terminating  the  sum  are  used 
in  different  calculations11*’15, 16  and  the  problem  is  discussed  in  some  detail  in 
Reference  6.  At  temperatures  in  the  reentry  range  all  these  methods  give  the  same 
results  for  the  partition  functions,  and  at  higher  temperatures  they  give  different 
results  only  for  species  which  are  almost,  completely  ionized16  and  therefore  present 
in  the  mixture  in  very  small  abundance. 

The  neglect  of  interparticle  forces  (i.e. ,  use  of  the  ideal-gas  assumption  for  the 
component  species)  leads  to  errors  in  the  calculations  at  high  densities.  In  addition 
to  the  usual  virial -coefficient  corrections19’17,  long-range  coulombic  forces  become 
important  at  the  high  temperatures  where  appreciable  ionization  takes  place18*27.  The 
compressibility  factor,  z  given  in  Equation  (2-7),  is  no  longer  equal  to  the  relative 
increase  in  the  number  of  moles  of  gas,  but  includes  the  effects  of  these  intermolecular 
forces.  For  nitrogen18  at  2000°K  and  10  times  normal  density,  z  is  increased  by  3%. 
and  at  100  times  normal,  by  15%.  At  15,000°K  the  effect  is  a  20%  increase  at  100  times 
normal  density.  Another  effect  of  these  forces  is  to  slightly  decrease  the  amount  of 
dissociation  and  ionization,  so  that  at  some  conditions  of  high  temperature  and  inter¬ 
mediate  density  the  effect  can  be  to  actually  decrease  the  value  of  z  . 

The  final  calculation  of  the  thermodynamic  properties  of  a  gas  from  equations  such 
as  (2-7)  and  (2-8)  are  relatively  simple  once  the  concentrations  jj  have  been  obtained. 
To  obtain  equations  for  the  specific  heats  and  the  speed  of  sound  requires  that  a  second 
derivative  of  the  partition  functions  be  obtained,  and  the  resulting  formulas  are 
lengthy,  but  tractable15’16.  The  equilibrium  speed  of  sound  for  air  is  tabulated  in 
References  15  and  28  and  given  graphically  in  the  Mollier  chart  which  accompanies 
Reference  23.  References  24  and  25  contain  these  properties  for  oxygen  and  nitrogen. 

A  Mollier  diagram  which  includes  the  effect  of  second  virial  coefficients  is  given  in 
Reference  29. 

2.3  Methods  for  Equilibrium  Calculations 

The  use  of  tables  of  thermodynamic  properties  provides  obvious  limitations  in  per¬ 
forming  gas-dynamic  calculations.  To  circumvent  this  difficulty  it  is  possible  to  use 
a  simplified  analytic  model  which  applies  to  a  dissociating  diatomic  gas,  or,  in 
machine  calculations,  to  use  a  sub-routine  capable  of  producing  the  thermodynamic 
properties  of  very  complicated  gas  mixtures. 
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2.3.1  Ideal  Dissociating  Gas 

The  equations  for  a  reacting  gas  simplify  considerably  for  the  special  case  of  a 
dissociating  diatomic  gas4.  They  can  be  simplified  still  further  by  the  approximation 
due  to  Lighthill30  that,  over  a  broad  range  of  temperatures,  the  vibrational  partition 
function  can  be  taken  as  proportional  to  the  square  root  of  the  temperature,  or 


(2-9) 


This  is  equivalent  to  the  assumption  tnat  the  vibrational  mode  is  "half -excited",  so 
that  it  contributes  &RT  to  the  enthalpy  per  mole.  An  approximation  to  the  equili 
brium  constants,  Equation  (2-6),  can  be  written  as 


„  -V.  -D'/kT 
=  CiT  l  e  i' 


(2-10) 


where  is  the  energy  of  the  reaction,  given  in  the  notation  of  Equation  (2-6)  as 


Di  _  j?i  ^iJ  “  viPDj 


(2-11) 


Hie  values  of  the  constants  C,  17  and  D  for  the  equilibrium  constants  for  high 
temperature  air  have  been  given  by  Wray27’4.  With  the  final  assumption  that  the 
electronic  partition  function  can  be  taken  as  a  constant  equal  to  the  degeneracy  of 
the  ground  electronic  state,  the  equilibrium  equations  can  be  written  in  a  much  simpli¬ 
fied  form  that  still  retain  the  important  physical  features  of  the  dissociating  gas- 
dynamics.  For  this  reason  the  Lighthill  model  has  been  used  in  many  analytic  approaches 
to  real -gas  problems. 


2.3.2  Methods  for  Sub-Routine  Use  with  Other  Calculations 

Because  of  the  great  calculation  speed  now  available  with  digital  computing  machines, 
it  is  possible  to  solve  numerical  gasdynamic  problems  using  sub-routines  to  calculate 
real-gas  properties  within  the  flow.  Several  such  programs  have  been  published  in 
detail31'34. 

In  some  of  the  programs  the  enthalpy  for  each  species  is  calculated  by  a  curve-fit 
formula  which  is  a  polynomial  in  the  temperature31'33,  rather  than  by  calculating  the 
derivatives  of  the  partition  function,  as  indicated  in  Equation  (2-8).  The  data  for 
the  curve-fit  can  be  obtained  from  .he  JANAF  tables10  for  many  chemical  species.  Hie 
conversion  of  the  JANAF  data  tor  thurmo-fit  purposes  is  described  in  Reference  33. 

Other  of  the  programs  use  approximate  formulas  for  the  internal  partition  functions 
and  calculate  the  enthalpy  directly.  These  calculations  are  then  used  as  the  starting 
point  to  solve  for  the  equilibrium  conditions35’36. 

Graphs  of  the  thermodynamic  properties  of  air  and  some  other  species  are  available 
for  graphical  solution  of  shock-wave37  or  expanding-flow23  problems. 
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SECTION  3 

BASIC  NONEQUILIBRIUM  RATE  PROCESSES 


The  conventional  treatment  of  aerodynamic  problems  requires  the  assumption  that  the 
air  can  be  treated  as  a  continuum,  and  that  the  state  parameters  vary  continuously 
from  point  to  point,  or  from  time  to  time.  It  is  understood,  of  course,  that  the  air 
is  composed  of  molecules  with  a  distribution  of  velocities,  and  that  collisions  among 
these  molecules  are  the  sine  qua  non  which  change  the  density  and  temperature  and  other 
state  parameters  of  the  gas.  The  assumption  of  a  continuum  is  the  assumption  that  on 
any  scale  of  distance  or  time  in  which  there  is  an  appreciable  change  in  the  state 
parameters,  many  of  these  collisions  occur,  so  that  for  most  of  the  molecules  a  Boltz¬ 
mann  distribution  of  velocities  is  always  present  and  the  state  parameters  can  be 
defined  at  each  point  in  the  flow.  The  number  of  collisions  required  for  this  local 
equilibration  of  translational  motion  to  be  established  cannot  be  specified  exactly, 
but  some  number  less  than  ten  would  suffice,  since  the  transfer  of  translational  energy 
among  the  molecules  is  an  efficient  process.  The  use  of  the  equilibrium  relations 
discussed  in  Section  2  requires  a  further  assumption,  and  that  is  that  so  many  of  these 
collisions  occur  that  all  important  chemical  processes  reach  an  equilibrium  distribution. 
This  might  easily  require  10s  or  more  collisions,  since  the  occurrence  of  many  chemical 
processes  upon  collision  requires  special  conditions  -  either  a  very  energetic  collision, 
or  some  special  geometry  of  the  impact,  or  a  sequence  of  several  processes  occurring  on 
different  collisions.  Thus,  there  are  many  aerodynamic  problems  of  interest  where  the 
first  assumption  is  satisfied,  so  ihat  the  flow  can  be  treated  as  a  continuum  and  a 
local  value  of  the  temperature  can  be  defined,  but  the  second  assumption  is  violated, 
so  that  equilibrium  relations  among  the  state  parameters  cannot  be  used.  It  is  the 
purpose  of  this  section  to  describe  the  "slow”  physical  processes  which  occur,  to  formu¬ 
late  the  nonequilibrium  equations  to  be  used  under  these  conditions  in  place  of  the 
equilibrium  thermodynamic  relations,  and  to  indicate  their  general  importance  in  aero¬ 
dynamic  problems. 

The  internal  motions  of  a  molecule  can  be  described  in  terms  of  rotation  and  vibration 
of  the  nuclei  and  orbital  motion  of  the  electrons.  Extreme  excitation  of  rotation  and 
vibration  leads  to  separation  of  the  nuclei,  or  dissociation,  while  extreme  excitation 
of  the  electronic  motion  results  in  ionization.  Upon  the  collision  of  some  atoms  and 
molecules  the  interchange  of  nuclei  occurs,  i.e. ,  there  is  a  chemical  reaction.  The 
energies  that  are  involved  with  the  various  processes  (for  air)  and  therefore  the 
temperature  ranges  where  they  are  of  interest,  have  been  discussed  in  Section  2.  These 
same  processes  are  no"'  to  be  discussed  again,  this  time  with  a  view  to  the  rate  at 
which  they  approach  equilibrium  in  a  given  physical  situation  -  that  is,  the  average 
number  of  collisions  required  for  each  process  to  obtain  its  equilibrium  value. 

On  every  collision  there  is  a  transfer  of  momentum  between  the  colliding  particles, 
and  this  transfer  is  really  what  is  meant  by  a  collision.  Except  in  a  case  where  the 
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masses  of  the  colliding  particles  are  very  different  (such  as  an  electron-atom  collision), 
this  momentum  transfer  is  accompanied  by  appreciable  energy  transfer.  As  mentioned 
before,  only  a  few  collisions  are  required,  on  the  average,  to  establish  a  Boltzmann 
distribution  of  translational  motion  at  a  new  temperature  when  the  overall  translational 
energy  of  the  gas  is  changed.  As  noted  in  the  Introduction,  the  present  scope  is 
limited  to  those  flows  where  the  translational  motion  remains  in  a  local  equilibrium 
as  the  state  properties  of  the  flow  change.  The  reader  who  wishes  to  investigate  the 
problem  of  translational  nonequilibrium  is  referred  to  Chapters  IX  and  X  of  Reference  1, 
and  the  many  references  given  therein. 

Since  translational  equilibration  is  all  that  is  required  before  a  temperature  and 
pressure  can  be  defined  in  the  usual  way,  the  equation  of  state  relating  p  ,  p  ,  and 
T  (Equation  (2-7))  can  be  applied  at  all  points  in  the  gas  under  consideration. 

3.  1  Relaxation  of  Internal  Modes 

The  adjustment  of  the  rotational  and  vibrational  energies  of  molecular  nucleii,  and 
the  adjustment  of  orbital  energies  of  the  electrons  of  either  atoms  or  molecules,  will 
be  considered  in  this  section.  The  relaxation  of  a  single  mode  can  never  be  considered 
completely  apart  from  the  other  motions  since  energy  is  transferred  between  modes,  very 
often  with  high  efficiency.  As  one  mode  approaches  equilibrium  it  affects  the  energy 
distribution  in  a  second  mode,  and  this  energy  distribution  affects  the  rate  at  which 
the  first  process  proceeds.  The  two  modes  are  coupled.  However,  under  many  conditions 
a  useful  separation  can  be  made  according  to  the  characteristic  times  for  equilibration 
of  the  separate  modes.  The  translational  and  rotational  relaxations  are  by  far  the 
fastest.  The  next  most  rapid  may  be  the  exchange  of  vibrational  energy  among  molecules, 
both  between  different  vibrational  states  of  a  given  molecular  species  and  between 
different  molecules.  Next  is  vibrational  relaxation  where  the  energy  is  transferred 
from  translational  motion,  and  often  at  the  same  rate  (or  faster)  is  the  excitation  of 
the  high  electronic  states.  The  processes  requiring  the  longest  time  to  equilibrate 
are  dissociation  and  ionization,  and  these  are  treated  separately  in  Section  3.2. 

The  rotational  motion  of  the  gas  is  known  to  be  coupled  very  closely  to  the  trans¬ 
lational  motion,  and  to  adjust  with  almost  equal  rapidity  to  a  Boltzmann  distribution 
at  the  temperature  defined  by  the  translational  motion1.  It  is  not  surprising  that 
this  adjustment  is  rapid,  since  the  energy  spacing  of  the  rotational  levels  is  small, 
and  energy  can  be  passed  from  translation  to  rotation  in  small  increments.  Thus, 
almost  any  collision  is  capable  of  this  energy  transfer,  and  it  is  a  common  assumption 
that  the  rotational  motion  relaxes  as  rapidly  as  the  translational  motion,  so  that  the 
two  modes  remain  in  local  equilibrium  at  all  times.  Experimental24’25  and  theoretical 
values  of  the  rotational  collision  number  range  from  about  4  collisions  (at  room  tem¬ 
perature)  to  about  10  collisions  (at  1500°K)  for  oxygen  and  nitrogen,  and  are  less  for 
polar  gases29.  The  assumption  will  be  made  here  that  the  translational  and  rotational 
motions  cannot  be  effectively  separated.  The  treatment  of  this  problem  is  surveyed  in 
Reference  2.  The  high  temperature  conditions,  where  rotational  relaxation  occurs  on 
the  same  time  scale  as  vibrational  relaxation,  electronic  excitation,  and  dissociation 
is  illustrated  in  Reference  30. 

3.1.1  Vibrational  Relaxation 


The  vibrational  motion  of  molecules  adjusts  relatively  slowly  to  an  equilibrium 
distribution  when  the  translational  and  rotational  motions  of  the  molecule  are  disturbed. 


15 


The  time  required  for  this  adjustment  can  easily  be  of  the  order  of  the  flow  time  in 
aerodynamic  problems,  since  104  to  107  collisions  can  be  required3  for  the  adjustment 
of  the  vibrational  energy.  At  final  equilibrium  of  diatomic  molecules,  this  vibrational 

motion  may  contain  as  much  as  2/7  of  the  total  energy  of  the  molecules;  thus,  the  ques¬ 

tion  of  whether  or  not  this  energy  has  been  absorbed  cannoted  be  neglected,  and  the 
correct  evaluation  of  the  rate  of  vibrational  equilibration  is  necessary  for  the  des¬ 
cription  of  air  flows  which  reach  temperatures  greater  than  about  1500°K.  In  this 
section  the  experimental  and  theoretical  treatments  of  vibrational  relaxation  for  dia¬ 
tomic  molecules  are  described.  Much  of  this  treatment  applies  as  well  to  more 
complicated  molecules  with  more  than  one  vibrational  mode,  but  the  discussion  is  limited 
to  diatoms. 

The  reason  for  the  slow  adjustment  of  vibrational  motion  is,  basically,  that  the 
natural  period  of  vibration  of  a  molecule  is  short  in  comparison  with  the  time  required 
for  most  molecular  impacts  to  occur.  During  the  impact,  both  atoms  of  the  molecule 

experience  similar  forces  during  the  collision  -  the  atom  nearer  to  the  impact  particle 

being  acted  on  directly,  and  the  farther  atom  being  acted  on  through  the  "spring”  forces 
of  the  molecule.  Since  the  spring  is  relatively  stiff,  the  contraction  of  the  spring 
can  follow  the  motion  of  the  impacting  particle  adiabatically,  and  as  the  particle 
slowly  recedes  the  spring  can  return  slowly  to  its  original  length,  without  being  left 
with  any  residual  vibrational  energy.  Thus,  although  there  may  be  a  transfer  of  trans¬ 
lational  energy  between  the  incident  particle  and  the  molecule  (as  there  would  be  between 
two  colliding  atoms),  there  is  no  transfer  of  energy  from  translation  to  vibration.  Only 
in  an  unusually  energetic  collision  does  one  atom  of  the  molecule  gain  considerable 
velocity  relative  to  the  other,  resulting  in  a  net  transfer  of  vibrational  energy  to 
the  molecule.  Thus,  vibrational  energy  is  transferred  only  by  the  very  energetic  collis¬ 
ions,  and  it  follows  that  on  the  average,  a  large  number  of  collisions  is  required  for 
the  transfer  of  a  small  amount  of  vibrational  energy.  Molecules  with  greater  spring 
constants  (higher  natural  frequencies)  require  more  collisions  for  vibrational  equili¬ 
bration.  As  a  corollary  it  follows  that  at  higher  temperatures,  where  there  are  more 
highly  energetic  collisions,  the  efficiency  of  energy  transfer  to  vibration  will  increase. 

The  analytical  description  of  relaxation  of  any  process  in  a  gas  wherein  a  local 
translational  velocity  distribution  can  be  defined,  can  be  conveniently  divided  into 
two  quite  separate  problems,  and  this  separation  has  been  used  in  the  presentation 
which  follows  for  vibrational  relaxation.  Hie  first  has  to  do  with  the  mechanical  prob¬ 
lem  already  discussed,  determining  the  amount  of  vibrational  energy  transferred  in  a 
single  collision  for  which  the  initial  conditions  are  defined.  These  transfer  rates 
can  then  be  averaged  over  all  collision  angles  and  over  the  Maxwellian  distribution  of 
translational  velocities  that  are  assumed  to  exist  in  the  gas,  to  obtain  the  average 
rate  of  energy  transfer  to  a  molecule  with  a  given  initial  vibrational  energy  as  a 
function  of  the  temperature  of  the  gas.  Given  these  rates,  the  second  problem  is  to 
describe  the  time  history  of  the  distribution  of  vibrational  energy  during  the  relaxa¬ 
tion  process.  Rach  of  these  problems  has  received  a  great  deal  of  attention,  and  each 
has  been  solved  with  certain  restrictive  approximations.  These  solutions  have  provided 
the  framework  within  which  it  has  been  possible  to  order  and  compare  experimental 
results,  and  to  perform  calculations  for  conditions  similar  to  those  that  have  been 
tested  by  experiment.  However,  the  approximations  limit  the  credibility  with  which 
the  theory  can  be  applied  to  cases  where  experimental  verification  is  not  available, 
or  where  the  experiments  are  difficult  to  interpret. 
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The  approximate  theories  have  been  described  in  several  texts1*4  and  so  will  be 
presented  only  in  summary  fashion  here.  The  basic  material  that  is  discussed  here  in 
more  detail  has  to  do  with  limitations  of  the  approximations  involved  and  the  possi¬ 
bility  of  extension  beyond  these  approximations.  Thus  the  emphasis  here  is,  in  a  some¬ 
what  unbalanced  measure,  on  the  restrictive  applicability  of  the  theory  and  the  care 
with  which  it  must  be  applied  to  aerodynamic  problems  where  direct  experimental 
measurements  are  not  available. 

Vibrational  energy  transfer  in  collision 

The  collision  of  two  molecules  or  of  a  molecule  and  an  atom  is  a  many-body  problem 
involving  all  the  electrons  of  the  colliding  partners,  as  well  as  their  nuclei.  As 
such,  the  problem  is  not  tractable,  either  classically  or  quantum  mechanically,  and  no 
discussion  of  the  difficulties  will  be  given  here.  Instead  the  problem  is  at  first 
simplified  by  neglecting  the  kineticsof  the  electron  motion,  and  replacing  the  electrons 
with  a  potential  field  dependent  only  on  the  separations  of  the  various  nuclei.  This 
(Born-Oppenheimer)  separation  of  the  electron  and  nuclear  motions  has  good  theoretical 
and  experimental  justification,  and  cannot  be  considered  a  serious  limitation  of  vibra¬ 
tional  energy-transfer  calculations.  However,  it  must  be  recognized  that  because  of 
this  extreme  simplification,  empiricism  enters  the  problem  at  a  very  early  stage;  a 
mathematical  form  of  the  potential  must  be  chosen,  and  the  parameters  which  enter  the 
form  must  be  determined  from  experiment. 


The  portion  of  the  potential  associated  with  the  intra-nuclear  coordinate  R  can 
be  obtained  with  good  accuracy,  because  of  the  mass  of  experimental  data  available  in 
the  form  of  molecular  spectra.  There  are  several  standard  analytic  forms  that  are 
employed  for  this  potential,  the  most  familiar  probably  being  the  Morse  potential 


=  °o  L1  - 


-/3(R-Re) 


which  has  two  parameters, 


and  J3  .  These  can  be  related,  for  example,  to  the  two 


experimental  data  of  the  dissociation  energy  (D0)  and  the  vibrational  frequency  of  the 
molecule 

..  _  1  lA"A 


The  equilibrium  value  of  the  intranuclear  separation,  Re  ,  does  not  enter  into  the 
collision  problem.  An  important  feature  of  all  such  potentials  is  that,  for  R-Re 
small,  they  reduce  to  the  simple-harmonic  potential.  For  example,  Equation  (3-1)  is 
written  in  this  limit  as 

V  =  D0/3z(R  -  Re) 2  (3-3 
and  many  of  the  calculations  discussed  later  are  made  within  this  approximation. 


The  potential  describing  the  interaction  of  the  incident  atom,  a  distance  r  from 
one  of  the  molecular  nuclei,  is  usually  taken  as  a  simple  repulsive  exponential 3 1-3 4 

V  =  Ae’0^  =  E0e  ^  .  (3-4) 
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The  two  forms  are  equivalent,  but  the  second  illustrates  that,  in  an  end-on  collision 
with  initial  kinetic  energy  EQ  ,  the  turn-around  point  will  be  at  r  =  rt  .  The 
exponential  form  is  much  more  tractable  in  collision  calculations  than  the  more  detailed 
potentials  which  include  the  effects  of  long-range  attractive  forces  and  it  is  the 
short-range  repulsive  forces  which  determine  the  rate  of  energy  transfer.  Values  of 
A  and  a  can  be  obtained  from  scattering  experiments  or  from  virial-coefficient  data. 

In  Reference  33,  a  is  obtained  for  various  collision  partners  from  tabulated  values 
of  the  constants  for  the  more  common  Lennard-Jones  potential  by  equating  the  expressions 
for  dV/dr  at  r  =  r.  ,  for  the  two  potentials.  The  result  can  be  written  approximately 
as33 

a  ^  (17.5)/r0  .  (3-5) 

where  rfl  is  the  constant  appearing  in  the  Lennard-Jones  potential 


The  problem  of  a  monatomic  species  colliding  with  a  diatom  can  iow  be  formulated  in 
terms  of  these  potentials.  The  general  formalism  is  reviewed  by  Takayunagi6,  including 
reference  to  other  intermolecular  potentials.  For  the  case  where  the  velocity  vectors 
are  in  the  plane  of  the  three  particles,  numerical  solutions  have  been  obtained  by 
Benson  and  Berend27.  For  an  analytic  solution,  however,  a  much  more  drastic  simplifi¬ 
cation  is  required.  The  problem  is  solved  for  the  very  special  case  of  end-on  collisions 
and  thus  rotational  excitation  is  excluded  from  the  start.  The  mechanical  problem  is 
then  a  function  of  the  two  coordinates  given  in  Equations  (3-1)  and  (3-4),  the  distance  R 
between  the  molecular  nucleii,  and  the  distance  r  from  the  closer  nucleus  of  the  mole¬ 
cule  to  the  incoming  nucleus. 

It  is  clear  that  once  this  simplification  is  made  the  mechanical  problem  becomes 
immensely  easier.  Using  a  potential  which  is  a  sum  of  Equations  (3-3)  and  (3-4),  that 
is,  the  end-on  collision  of  an  atom  with  a  simple-harmonic  oscillator,  the  problem  has 
been  solved  classically26,  3S,  quantum-raechanically  32,  33,  and  semi-classically 3<.  Each 
of  these  different  methods  contain  some  further  approximation  which  limits  its  range 
of  applicability,  as  discussed  briefly  below.  However,  for  conditions  where  these 
ranges  overlap,  the  different  methods  give  essentially  the  same  answer  for  the  proba¬ 
bility  of  energy  transfer  in  a  collision  of  given  velocity  and  molecular  parameters, 
and  the  reasons  for  this  are  discussed  in  References  35,  37  and  38.  Whether  the  ex¬ 
ponential  repulsion  of  Equation  (3-4)  is  used  or  a  linearized  form  is  employed  has 
little  effect  on  the  result  for  nitrogen  up  to  incident  velocities  cf  9  x  105  cm/sec 
(and  thus  for  temperatures  of  10,000°K).  This  is  demonstrated  in  Reference  38  by  com¬ 
parison  of  the  solution  with  a  linear  potential  (for  which  an  exact  solution  to  the 
semi-classical  problem  can  be  obtained)  with  the  numerical  solution  for  an  exponential 
intranuclear  potential39.  A  recent  numerical  quantum-mechanical  calculation40  of  this 
problem  gives  results  which  are  in  some  disagreement  with  these  earlier  calculations, 
and  this  discrepancy  has  not  yet  been  explained.  The  use  of  an  infinite-step  potential 
in  place  of  Equation  (3-4)  results  in  a  considerably  different  result,  as  demonstrated 
in  References  40  and  41.  This  disagreement  would  be  expected,  in  view  of  the  qualita¬ 
tively  different  results  obtained  in  classical  mechanical  problems  when  impulsive 
collisions  are  considered. 
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The  most  radical  (and  most  fruitful)  further  simplification  that  can  be  made  in  the 
quantum-mechanical  problem  is  the  application  of  first-order  perturbation  theory,  in 
which  case  the  results  are  much  more  easily  formulated.  The  general  conclusions  that 
can  be  drawn  from  this  collision  model  (simple  harmonic  oscillator,  end-on  collision, 
transition  probabilities  much  less  than  1)  are42'  3 


1.  The  probability  that  in  a  collision  a  molecule  in  the  first  vibrational  state 
will  make  a  transition  to  the  ground  vibrational  state,  giving  up  energy  hr1  to 
translation  is  given  by 


10 


~  e 


-  Urrv/va 


(3-7) 


2.  The  probability  that  in  a  collision  a  molecule  in  state  r  +  1  will  make  a 
transition  to  state  r  is  related  to  p1Q  by 

Pr+ 1,  r  =  <r  +  DPio  •  (3-8) 

3.  Transitions  occur  only  between  adjacent  levels 


Pj  j  =  0  for  i  i  j  ±  1  . 


(3-9) 


If  first-order  perturbation  theory  is  not  used,  but  instead  the  semi-classical  method 
is  employed35*37,  each  of  these  conclusions  is  modified.  For  very  high-energy  col¬ 
lisions  (in  nitrogen  negligible  below  5000°K)  the  value  of  p1Q  given  in  Equation  (3-7) 
approaches  and  exceeds  1,  so  that  the  true  value  is  less  than  this,  and  in  addition 
multiple-jump  transition  probabilities  are  not  zero.  These  effects  compensate  in  such 
a  way  that  the  rate  of  energy  absorption  in  transitions  from  the  ground  state  can  still 
be  correctly  calculated37' 38  from  Equations  (3-7)  and  (3-9).  The  results  that  have  been 
obtained  for  Pr  0  hy  different  methods  are  compared  in  Figure  3-1  for  the  case  of  col¬ 
lisions  of  two  nitrogen  molecules.  A  three-dimensional  solution  has  been  obtained34 
for  the  collision  problem  using  a  quantum-mechanical  perturbation  theory  (method  of 
partial  waves)  in  which  the  potential  is  nearly  spherical,  but  the  spherical  form  is 
multiplied  by  a  term  which  depends  on  the  oscillator  coordinate.  The  result  is  very 
similar  to  the  corresponding  one-dimensional  calculation33.  A  classical  two-dimensional 
calculation28,  including  effects  of  rotational  energy  transfer,  demonstrates  that  rota¬ 
tional  motion  can  be  effectively  transferred  to  vibration. 


The  calculation  for  a  Morse  internuclear  potential  -  this  is,  using  Equation  (3-1) 
instead  of  Equation  (3-3)  presents  several  difficulties  that  are  not  resolved43'44. 

The  spacing  between  the  vibrational  energy  states  decreases  for  higher  energy  states, 
and  the  probability  of  a  transition  between  these  states  increases  because  of  the 
smaller  energy  required  in  the  transition  (see  Equation  (3-7))  and  because  the  matrix 
elements  for  a  Morse  oscillator  increase  rapidly45.  It  has  been  shown46  that  if  a 
(in  Equation  (3-4))  is  equal  to  /3  (in  Equation  (3-l)(Ref. 1)).  Equation  (3-9)  is  correct 
within  the  limits  of  first-order  perturbation  theory.  In  summary,  Equations  (3-7)  and 
(3-9)  may  be  reasonably  accurate  for  the  end-on  collision  at  moderate  velocity  of  a 
Morse  oscillator,  but  Equation  (3-8)  is  completely  incorrect,  For  very  high-velocity 
collisions,  Equations  (3-7)  and  (3-9)  must  be  modified  as  well38’39. 

In  addition  to  the  transfer  of  energy  between  vibrational  and  translational  motion, 
it  is  also  of  interest  to  consider  the  transfer  of  vibrational  energy  between  molecules 
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with  similar  vibrational  energy  spacings.  This  resonant,  or  near-resonant,  exchange 
of  vibrational  energy  can  occur  with  high  probability  (~  10"1),  providing  that  the 
energy  defect  for  the  transition  is  small.  A  review  of  available  data  for  transitions 
of  this  type  impcrtant  in  the  C02-N2  laser  has  recently  been  published47.  The  methods 
of  calculation  and  the  potentials  involved  are  similar  to  those  previously  described. 
Both  classical  and  quantum  mechanical  approaches  have  been  used4*  33,  50. 


Master  equations  for  the  time-history  of  the  vibrational  energy 

When  the  rates  for  vibrational  energy  transfer  in  specific  collisions  have  been 
determined,  there  remains  the  problem  of  describing  the  change  of  the  vibrational  energy 
of  the  gas  when  the  translational  energy  undergoes  a  change.  The  molecules  are  first 
grouped  according  to  their  vibrational  energy  state,  and  then  a  differential  equation 
(a  rate  equation)  is  written  for  the  rate  of  change  of  molecules  in  a  given  energy 
state  (say  the  rth)  due  to  collisions  of  all  the  molecules.  Thus,  the  expression  for 
the  transition  probability  pr+1  r  ,  as  given  by  Equations  (3-7)  and  (3-8),  must  be 
averaged  over  all  velocities,  weighted  according  to  the  velocities  that  occur  in  a 
Boltzmann  distribution  of  translational  motion  at  a  temperature  T  .  The  rate  constants 
derived  in  this  way  are  then  functions  of  the  temperature  rather  than  the  velocity  of 
any  given  collision.  The  differential  equation  utilizing  this  rate  constant  contains 
a  separate  term  for  each  type  of  molecular  pair  that  can  collide  which  will  result  in 
the  production  (or  destruction)  of  a  molecule  in  the  rth  state,  with  the  other  molecule 
having  the  choice  of  remaining  in  its  original  state  or  making  a  transition  to  another 
state.  If  the  gas  being  considered  is  a  mixture  of  gases  A,  B,  C,  etc,  the  equation 
is  further  complicated  by  the  requirement  of  keeping  account  of  the  species  involved 
in  the  collision.  We  begin  by  writing  this  equation  in  some  detail,  after  which  it 
will  be  simplified  by  various  approximations  so  as  to  obtain  a  solution. 


Consider  a  mixture  of  diatomic  gases  A  and  B,  and  consider  transitions  where  only 
one-quantum  transitions  are  involved,  consistent  with  Equation  (3-9).  The  differential 
equation  for  the  number  density  of  molecules  of  Type  A  in  the  rth  state  is  then 


+ 


n  0'^tr+i+eR- r*r'*R)'/kT  N  1 

nr+iNR-i  "  e  nrNR 

m  ‘(er+€R'eR+rer+i)/kT  N 
nrNR  “  e  nr- inr+i  | 


(3-10) 
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A  term  such  as  nr+ins-i  describes  the  rate  at  which  molecules  enter  the  rth 

state  from  the  (r  +  i)  state,  due  to  collisions  of  molecules  in  the  (s  -  1)  state  with 
A-type  molecules  in  the  (r  +  1)  state,  where  the  same  collision  provides  a  transfer  of 
the  molecule  in  the  (s  -  1)  state  into  the  (s)  state.  Capital  letters  refer  to  the 
second-type  molecule.  The  equation  includes  the  substitution  of  the  detailed-balance 
relation,  namely,  that 


ps, s-i 
r.  r+i 


'(er+i+<s-i'er‘£s)/kTpS-i, s 
e  rr+i,r  • 


(3-11) 


In  writing  this  relation,  and  in  writing  the  rate  constants  as  a  function  of  the  trans¬ 
lational  temperature  alone,  there  is  the  assumption  that  all  other  modes  of  excitation 
are  in  equilibrium  at  the  translational  temperature.  If  this  were  not  the  case  then 
more  terms  would  be  required  in  Equation  (3-10)  accounting  for  the  exchange  of  vibra¬ 
tional  energy  with  the  other  modes,  and  a  separate  detailed-balance  relation  would  be 
required  for  each  of  these  separate  reactions.  With  gases  at  moderate  temperatures, 
where  a  clear-cut  distinction  can  be  made  between  translational-rotational  relaxation 
time  and  vibrational  relaxation  time,  this  restriction  is  not  a  serious  one.  At  higher 
temperature,  where  dissociation  and  chemical  reactions  are  involved,  as  discussed  pre¬ 
viously  this  separation  becomes  more  difficult,  and  a  complete  description  cannot  be 
obtained  without  taking  these  other  nonequilibrium  effects  into  account. 


Vibrational  relaxation  of  a  dilute  mixture  of  molecules  and  atoms 


A  great  simplification  of  Equation  (3-10)  is  obtained  in  treating  the  very  special 
case  of  vibrational  relaxation  of  a  diatomic  molecule  in  a  very  dilute  mixture  with  a 
monotomic  species.  If  A  is  the  diatom  and  B  is  the  monatom,  we  can  consider  a 
mixture  so  dilute  that  molecule-molecule  collisions  can  be  neglected,  and  Equation  (3-10) 
simplifies  to 


dnr 

dt  =  Wr»|>i  " 


r(er+rer)/kT 


T  -(£r'er-l>/kT 

nr  "  Qr, r-iN  nr  ’  e  nr-i  •  (3'12> 


The  next  step  in  the  solution  requires  the  specification  of  the  dependence  of  the  rate 
constants  Qr+1  r  on  the  energy  states  r  +  1  and  r  .  As  discussed  in  the  previous 
section,  this  dependence  can  be  obtained  for  the  case  where  the  molecule  is  properly 
described  as  a  simple-harmonic  oscillator  by  averaging  Equations  (3-7)  and  (3-8)  over 
a  Boltzmann  velocity  distribution,  resulting  in3’42 

Q10  ~  e"(5,t7r4l/Zs/kT)  1/3  ,  (3-13) 

<W  =  (r+1)Qio- 

where  m  is  the  reduced  mass  of  the  system,  giver,  for  an  incoming  particle  A  strik¬ 
ing  a  molecule  B-C  ,  as 

raA(mB  +  mc) 

m  -  -  • 

iD^  ^  nig  ^  iiIq 

Using  Equation  (3-14),  the  f irst-moment  equation  can  be  obtained42  from  Equation  (3-12) 
by  multiplying  both  sides  of  the  equation  by  r  (and  by  h v/p)  and  summing  over  all 
values  of  r  .  There  results 


ev  EQ  ev 


21 


fey 

dt 


(3-15) 


where  ey  is  the  vibrational  energy  per  unit  mass  of  gas  related  to  nr  by 

hi"  ® 

e  =  —  til.  (3-16) 

p  r»o 


with  p  equal  to  the  density  of  the  gas.  ey  EQ  is  the  vibrational  energy  per  unit 
mass  that  would  be  contained  by  the  gas  if  it  were  at  equilibrium  at  the  local 
temperature  T  .  ry  is  a  relaxation  time  that  is  related  to  Q10  through 


Qi0(l  -  e 


-£d/kT 


(3-17) 


These  extremely  simple  equations  have  found  widespread  use  in  spite  of  the  limitations 
required  for  their  derivation. 

A  detailed  solution  of  Equations  (3-12)  to  (3-14)  has  been  given51  for  the  case 
where  T  is  constant  during  the  relaxation.  This  limitation  is  consistent  with  the 
previous  assumption  that  the  diatomic  species  is  strongly  diluted  with  a  monatomic 
species.  If  the  gas  starts  in  an  initial  Boltzmann  distribution  at  temperature  TQ  , 
its  distribution  can  be  described  at  any  latter  time  by51 

-«r/kT  -rhv/kT 

nr  =  nQe  =  n0e  ,  (3-18) 


where  Ty  is  a  time -dependent  parameter  that  starts  at  a  value  T0  and  approaches  a 
value  T  as  the  relaxation  proceeds.  That  Equation  (3-18)  is  a  solution  of  Equations 
(3-12)  to  (3-14)  can  be  demonstrated  by  direct  substitution.  Thus,  the  vibrational 
population  remains  in  a  Boltzmann- like  distribution  throughout  the  relaxation,  and  Ty 
is  referred  to  as  the  vibrational  temperature.  It  is  related  to  the  vibrational  energy 
per  gram  by  the  usual  formula  for  the  energy  of  a  system  of  simple-harmonic  oscillators, 
viz. , 


n  hi" 

ev  =  “  ehv/kTv  _  j  ' 


(3-19) 


The  departure  of  real  molecular  oscillators  from  the  simple -harmonic  model  has  an 
important  effect  on  Equation  (3-14)  for  the  transition  probabilities  of  the  higher- 
level  transitions,  as  previously  described.  For  many  cases  of  aerodynamic  interest, 
almost  all  of  the  vibrational  energy  is  contained  in  the  first  few  energy  states,  so 
that  this  change  in  transition  probability  for  higher  states  has  little  effect52  on 
the  validity  of  Equations  (3-13),  (3-15)  and  (3-17).  However,  Equation  (3-18)  must  be 
applied  with  considerable  caution  for  higher  states,  and  in  problems  of  rapidly  ex¬ 
panding  flows  this  can  affect  the  validity  of  Equations  (3-15)  and  (3-17). 

Comparison  of  experimental  results  with  the  predictions  of  Equations  (3-13),  (3-15), 
(3-17),  and  (3-18)  is  given  in  Section  4.3. 
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Vibrational  relaxation  of  a  pure  gas 

To  describe  the  vibrational  relaxation  of  a  pure  gas  it  is  necessary  to  include 
terns  from  Equation  (3-10)  that  describe  both  the  transformation  of  translational 
energy  to  vibration,  and  the  exchange  of  vibrational  energy  between  molecules  in  a 
collision.  Equation  (3-10)  can  then  be  written 


A  unique  feature  of  this  equation  which  aids  in  its  approximate  solution  is  that  terms 

5Z  PrJ.i+r  “r+i1^  are  much  lar*er  than  the  terms  Pr+1  rnnr+1  .  That  is,  there  is  a 

much  higher  probability,  in  a  collision,  that  a  vibrational  quantum  will  be  exchanged 
between  the  colliding  molecules  than  there  is  that  the  vibrational  energy  will  be  trans¬ 
ferred  to  translational  motion.  One  molecule  makes  a  transition  to  the  next  lower 
vibrational  state,  while  the  other  molecule  makes  a  transition  to  the  next  higher  state. 
Transitions  of  this  type  do  not  change  the  total  number  of  vibrational  quanta  in  the 
gas,  defined  (per  unit  mass)  by 

U  =  -  Tmr  ,  (3-21) 

P 


but  they  do  rearrange  the  quanta  among  the  vibrational  states.  These  rapid  exchange 
collisions  attempt  to  establish  a  local -equilibrium  distribution  of  molecules  in  the 
vibrational  states  which  will  make  the  brackets  in  the  last  two  lines  of  Equation  (3-20) 
separately  equal  to  zero  for  a  given  number  of  vibrational  quanta.  It  can  be  seen  by 
direct  substitution  in  Equation  (3-20)  that  this  is  accomplished  by  the  distribution53’ 54 


=  noe’rre 


-«r/kT 


(3-22) 


where  y  is  determined  by  the  total  number  of  quanta  available,  evaluated  by  substi¬ 
tuting  (3-22)  into  (3-21).  Equation  (3-22)  reduces  to  a  Boltzmann-like  distribution 
at  a  temperature  Ty  for  a  simple  harmonic  oscillator  so  that  within  this  approxima¬ 
tion  the  exchange  of  vibrational  quanta  does  not  affect  relaxation,  and  Equations  (3-15) 
to  (3-17)  should  apply55’56.  For  an  anharraonic  oscillator,  the  departure  from  a 
Boltzmann-like  distribution  depends  on  the  number  of  vibrational  quanta  contained  by 
the  gas  relative  to  the  equilibrium  number  at  the  translational  temperature  T  .  If 
the  number  of  quanta  is  less  than  the  equilibrium  number  (as  during  vibrational  excita¬ 
tion),  the  upper  states  are  underpopulated,  and  if  the  number  of  quanta  is  greater  than 
the  equilibrium  number  (as  during  vibrational  deexcitation)  the  upper  states  are  over- 
populated. 
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For  the  topmost  quantum  numbers  the  distribution  given  by  Equation  (3-22)  becomes 
invalid  because  the  assumption  that  vibration-vibration  exchange  is  much  more  rapid 
than  vibration-translation  energy  transfer  becomes  incorrect  for  two  reasons;  these 
states  are  very  far  from  resonance  with  the  lower,  more  densely  populated  states,  thus 
decreasing  the  probability  of  vibrational  exchange,  and,  in  addition,  the  probability 
of  transfer  of  energy  from  vibration  to  translation  increases  to  an  appreciably  greater 
value  in  the  highest  states  because  of  the  closer  spacing  of  the  levels. 


With  these  restricMons  in  mind,  a  differential  equation  for  the  time  history  of  U 
can  be  obtained  by  multiplying  Equation  (3-20)  by  r/p  and  summing  over  all  states. 
The  terms  involving  exchange  do  not  contribute  to  the  sum,  since  these  exchanges  of 
vibrational  quanta  do  not  affect  the  total  number  of  quanta.  The  result  is5*' 57 


dl) 

dt 


-<Ver-i>/kT 


nr- 1  • 


(3-23) 


This  equation,  along  with  the  distribution  given  by  Equation  (3-22)  provides  a  closed 
set  which,  in  principle,  determines  l)  (or  y)  as  a  function  of  time.  The  solution 
requires,  however,  a  knowledge  of  all  the  values  of  Pr+1  r  ,  and  these  are  not  known 
with  any  certainty.  In  addition,  the  restrictions  at  high  quantum  number  on  the  dis¬ 
tribution  given  by  Equation  (3-22)  become  important  in  problems  involving  low  trans¬ 
lational  temperatures  and  high  vibrational  energies.  Some  attempts  at  the  solution  of 
Equation  (3-23)  are  given  in  References  54,  57,  58,  52,  and  56.  The  general  conclusions 
obtained  there  are  that  vibrational  excitation  of  a  pure  gas,  where  the  molecules  are 
described  as  anharmonic  oscillators  with  rapid  vibrational  exchange,  can  be  adequately 
described  by  assuming  a  Boltzmann-like  distribution  during  excitation  and  employing  the 
Landau-Teller  model.  However,  during  deexcitation,  the  overpopulation  of  upper  vibra¬ 
tional  levels  by  exchange  effects  could  increase  the  rate  at  which  quanta  are  lost, 
and  it  is  suggested  that  this  may  provide  the  reason  for  the  observed  rapid  vibrational 
deexcitation  in  nozzle  flows  (see  discussions  in  Section  4.3).  The  conditions  under 
which  vibrational  exchange  may  enhance  vibrational  deexcitation  are  given  in  References 
54  and  57  as: 

1.  The  translational  temperature  must  be  less  than  the  characteristic  vibrational 
temperature  of  the  gas  to  ensure  that  the  vibrational  exchange  mechanism  is 
dominant. 

2.  The  initial  vibrational  temperature  must  be  at  least  a  considerable  fraction  of 
the  characteristic  vibrational  temperature  if  there  is  to  be  any  significant 
influence  of  the  excited  vibrational  states. 

3.  The  ratio  of  the  initial  vibrational  temperature  to  the  translational  temperature 
must  be  considerably  greater  than  one,  for  the  anharmonicity  of  oxygen  or  nitrogen, 
so  that  Equation  (3-20)  represents  a  large  departure  from  a  Boltzmann-like 
distribution. 


These  conditions  are  considered  conservative,  in  the  sense  that  the  non-Boltzmann 
effect  may  also  become  noticeable  for  conditions  closer  to  equilibrium. 
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Numerical  solutions  for  Equations  (3-20)  are  given  for  12-  level  and  24-  level 
molecules  in  Reference  59,  with  assumed  values  for  P®|®+1  and  Pr+1  .  These 
solutions  demonstrate  the  time  dependence  of  the  vibrational-energy  distribution 
during  the  high-speed  vibrational-exchange  relaxation. 


Vibrational  relaxation  of  a  mixture  of  gases 

In  the  relaxation  of  two  or  more  diatomic  or  polyatomic  gases,  exchange  processes 
between  nearly  resonant  modes  can  be  of  primary  importance60,61.  Measurements  of  the 
simultaneous  relaxation  of  N2  and  NO  is  interpreted  in  terms  of  the  relative  im¬ 
portance  of  vibrational  exchange  and  vibration-translation  transitions  in  Reference  60. 
An  effective  relaxation  time  is  determined  from  the  separate  relaxation  rates4. 


The  equations  and  concepts  discussed  in  the  previous  section  can  be  applied  to  such 
mixtures  of  gases  with  the  same  type  of  restrictions  and  limitations  as  given  for  the 
case  of  a  pure  gas.  In  the  case  of  two  near-resonant  gases,  the  full  equation  given 
as  Equation  (3-10)  applies.  The  steady-state  distributions  established  by  exchange53,  54, 
as  given  in  Equation  (3-22),  applies  to  both  gases  with  the  same  value  of  y  in  the 
two  distributions.  Thus,  we  can  write 


nr  =  n0e'yre 


-er/kT 


Nr  =  N0e-**e 


(3-24) 


That  these  relations  reduce  all  the  exchange  terms  in  Equation  (3-10)  to  zero  can  be 
seen  by  direct  substitution.  As  for  a  single  gas,  the  simple-harmonic-oscillator 
approximation  yields  a  Boltzmann-like  distribution  with  some  vibrational  temperature 
Ty  .  However,  for  gases  with  different  fundamental  frequencies  their  vibrational  tem¬ 
peratures  will  not  be  equal.  These  distributions  are  discussed  further  in  References  53 
and  54.  The  numerical  calculation  of  the  time  dependence  of  these  distributions  is 
described  in  Reference  62. 


3.1.2  Electronic  Excitation 

The  orbital  motion  of  electrons  in  atoms  and  molecules  contributes  a  relatively 
small  amount  to  the  total  energy  of  a  high-temperature  gas.  This  is  demonstrated  for 
air  in  Figure  2-1.  The  importance  of  the  electronic  excitation  is  not  related  to  the 
total  energy  that  resides  there  at  any  one  time,  but  rather  to  the  many  ways  this 
energy  can  be  transferred  within  the  gas.  It  is  closely  coupled  with  vibrational 
excitation  in  molecules,  it  is  an  intermediary  step  in  the  ionization  process,  and  it 
is  the  starting  point  for  the  radiative  process.  The  fact  that  the  radiative  intensity 
is  proportional  to  the  population  of  the  excited  electronic  states  has  made  it  important 
to  understand  this  excitation  process  in  correlating  the  observed  radiation  from  high¬ 
speed  vehicles  with  the  physical  description  of  their  flight  trajectory. 

The  experimental  observation  of  excited  electronic  states  is  accomplished  directly 
by  observation  of  the  light  emitted  from  these  states.  However,  the  amount  of  data 
available  for  electronic  excitation  rates  in  the  high-temperature  range  is  much  less 
than  that  previously  discussed  for  vibrational  excitation.  The  chief  reason  for  this 
is  that  the  temperature  range  of  interest  is  typically  much  higher,  so  that  the  excita¬ 
tion  processes  occur  more  rapidly  and  are  less  accessible  to  shock-tube  instrumentation. 
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The  theoretical  interpretation  of  the  experimental  results  is  also  much  more 
difficult  than  for  the  vibrational  relaxation  case.  There  is  no  simplification  avail¬ 
able  for  low-lying  electronic  states  parallel  to  the  simple-harmonic-oscillator  approxi¬ 
mation  for  vibrational  states.  In  addition  there  are  many  vibrational  states  with  lower 
energy  than  even  the  low-lying  electronic  states,  and  so  coupling  with  these  and  other 
electronic  states  must  be  considered  in  the  excitation  process.  Finally,  chemical 
reactions,  particularly  dissociation  and  ionization,  occur  in  the  temperature  range 
where  electronic  excitation  is  observed  and  are  necessarily  involved  in  the  interpre¬ 
tation  of  the  experimental  results. 

Several  recent  reviews  of  electronic  excitation  rates  are  available7*8.  In  these 
reviews  it  is  emphasized  that  high-temperature-rate  data  for  the  excitation  of  elec¬ 
tronic  states  of  atoms  and  molecules  by  heavy-particle  impact  are  relatively  few,  and 
that  the  excitation  processes  pertinent  to  flight  problems  are  complicated  and  not 
well  understood.  Much  of  the  available  information  about  molecular  electronic  excita¬ 
tion  comes  from  fluorescence  and  electron-bombardment  experiments  performed  at  room 
temperature,  and  cannot  reliably  be  extrapolated  to  high  temperature  conditions.  Tne 
data  on  excitation  by  impact  of  free  electrons7*9  ars  obtained  with  molecules  at  room 
temperature,  and  thus  apply  to  low-lying  rotational  and  vibrational  states. 

The  limited  information  available  from  shock-tube  experiments  provides  very  little 
guidance  for  predictions  of  electronic  excitation  in  air63*64,  although  for  pure 
nitrogen  sufficiently  detailed  experiments  have  been  performed  to  correlate  them  with 
a  simple  model  for  excitation  of  the  A  and  B  states  of  N2  by  N2  and  N  collisions65*66. 
Similar  detailed  spectroscopic  investigations67  for  shock -heated  C02  involve  too 
many  unknown  chemical  reactions  to  permit  evaluation  of  mechanisms  or  rates  for 
excitation  of  electronic  states. 

Chemical  excitation  of  electronic  states  has  long  been  a  subject  of  fundamental 
interest  in  flame  spectroscopy10*68  and  other  experimental  investigations69*70,  and 
has  found  recent  application  to  chemical  lasers71.  Some  application  has  been  made  to 
reentry  physics  problems  in  connection  with  radiation  from  the  flow  field  of  blunt 
hypersonic  reentry  vehicles72. 

3.2  Chemical  Reactions 

Processes  which  involve  a  change  in  chemical  structure  of  the  molecules  of  the  gas 
usually  occur  at  a  rate  much  slower  than  those  involving  only  internal  readjustments. 

The  reactions  must  continue  until  an  equilibrium  distribution  as  described  by  Equations 
(2-5)  is  reached,  and  this  equilibration  requires  the  production  of  the  new  species  in 
all  its  permitted  equilibrium  energy  states.  These  products  will  be  produced  with 
various  efficiencies  from  the  many  permitted  states  of  the  parent  species.  The  kinetics 
of  attaining  these  final  states  can  involve  a  complicated  route  through  many  inter¬ 
mediate  states,  including  some  whose  equilibrium  population  is  so  small  that  the  states 
can  be  neglected  in  equilibrium  calculations.  Thus  a  complete  mathematical  description 
of  chemical  reactions  would  take  sets  of  equations  much  more  complicated  than  Equations 
(3-10)  that  were  used  to  describe  vibrational  relaxation. 

It  is  fortunate  for  the  subject  of  chemical  kinetics  that  it  is  not  necessary  to 
include  all  the  intermediate  steps  in  the  description  of  a  chemical  process.  Often 
only  the  starting  and  ending  species  need  be  included.  In  fact  much  of  the  experimental 
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and  theoretical  work  in  chemical  kinetics  is  directed  toward  finding  the  minimum 
mathematical  complexity  to  describe  the  observed  species  production.  It  must  be  under¬ 
stood,  however,  that  the  application  of  the  chemical  rates  and  mechanisms  to  environ¬ 
ments  much  different  from  those  where  experiments  were  performed  requires  caution. 

The  basic  equations  of  chemical  reactions  in  flow  systems  were  summarized  in  1955 
by  S.S.Penner11  and  they  have  since  been  discussed  in  several  aeronautical  texts1,5. 

No  attempt  is  made  here  to  repeat  the  derivations  of  the  basic  equations,  but  as  with 
the  basic  equations  of  Section  2  they  are  discussed  briefly  and  presented  for  reference 
later.  Several  items  of  special  interest  are  discussed  in  more  detail. 


3.2.1  Reaction  Rate  Equations 

Continuing  with  the  notation  of  Reference  11,  a  chemical  reaction  is  described  by 
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where  the  quantities  kpi  and  kRi  are  the  forward  and  reverse  rate  constants  for 
the  i  reaction  and  are  taken  to  be  functions  only  of  the  translational  temperature 
of  the  gas.  The  general  differential  equation  to  describe  the  rate  of  production  of 
species  j  is 
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There  is  in  these  equations  the  implication  that  the  collision  of  the  molecules  on 
the  left-hand  side  of  Equation  (3-25)  will  result  in  the  products  on  the  right-hand 
side,  the  probability  of  the  reaction  being  related  to  the  rate  constant  kpl  .  The 
reaction  does  not,  of  course,  occur  in  a  one-step  process,  but  instead  the  molecules 
absorb  rotational  and  vibrational  energy  in  a  series  of  collisions,  and  finally  the 
excited  molecules  are  involved  in  a  sufficiently  energetic  collision  for  the  reaction 
to  occur.  Likewise,  the  inverse  reaction  occurs  in  a  large  number  of  separate  steps. 
However,  a  description  that  is  adequate  for  most  purposes  can  be  formulated  by  ignoring 
the  intermediate  steps. 


Under  the  assumption  that  Equations  (3-25)  and  (3-26)  properly  describe  the  physical 
processes,  the  equilibrium  constants  Kj  discussed  in  Section  2,  can  be  related  to 
the  forward  and  reverse  rate  constants.  In  equilibrium  (dyj/dt)  =  0  and  from  the 
principle  of  detailed  balance,  each  of  the  terms  in  the  sum  on  the  right-hand  side  of 
Equation  (3-26)  equals  zero.  It  follows  that 
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The  accuracy  of  this  relation  when  the  details  of  intermediate  steps  in  the  reaction  are 
included  will  be  discussed  later. 
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Machine  programs  have  been  written  for  the  simultaneous  solution  of  Equations  (3-26) 
and  the  flow-field  equations  for  a  variety  of  flow  situations,  such  as  normal  shock 
waves,  quasi-one-dimensional  expansions,  external  flows  including  boundary  layer  effects 
and  ablation  effects,  and  wake  flows.  These  various  programs  will  not  be  discussed 
here,  but  references  are  given  in  the  Appendices.  Some  discussion  of  expansion-flow 
calculations  is  given  in  Section  4.  For  normal -shock  flows  a  large  number  (~  50)  of 
chemical  reactions  can  be  included  with  only  moderate  calculation  time  required. 

However,  for  more  complicated  flow  situations  the  calculation  time  can  increase  rapidly, 
so  that  the  amount  of  chemistry  which  can  be  included  must  be  reduced. 

The  temperature  dependence  of  the  forward  and  reverse  rate  constants  is  a  prime 
consideration  of  experimental  and  theoretical  investigation.  For  the  atomic  recombina¬ 
tion  reaction  to  occur,  there  is  no  appreciable  activation  energy  required.  A  reaction 
of  this  type  (which  we  write  as  a  reverse  reaction)  then  has  very  little  temperature 
dependence,  and  can  be  described  as 

kRi  =  CRiT7!R1  *  (3-28) 

where  77,  is  a  Suiall  constant.  The  forward  rate  constant  is  related  through  Equation 
(3-27).  Within  the  approximation  of  Equation  (2-10)  for  K,  kpi  is  then  of  the  form 

-D,/kT 

,  kFi  =  CpiT  “e  (3-29) 

Suitable  constants  for  these  relations  for  air  reactions  are  given  in  Reference  1. 


3.2.2  Dissociation 


The  application  of  Equations  (3-25)  and  (3-26)  to  the  dissociation  of  a  symmetric 
diatomic  molecule  provides  a  simple  example  which  has  been  studied  extensively  and 
which  involves  several  of  the  problems  associated  with  high-temperature-air  chemistry. 
Equation  (3-25)  becomes 


A,  +  A,  £  A.  +  2A 
k 

*R1 


A  +  A  J±  3A 

2  k 

KR2 


(3-30) 


and  Equation  (3-26)  is 

_ 


dt 


-  (-kFl^A2  +  +  ^~kF2^A2^A +kR2^A^  (3-31) 


These  equations  have  been  used  to  analyze  shock-tube  experiments  for  pure  02 
(Ref. 73)  and  N2  (Ref.  74)  an  well  as  many  other  gases  .  The  assumption  that  the 
rate  constants  can  be  taken  as  functions  of  temperature,  independent  of  the  state  of 
internal  excitation,  is  open  to  question  on  any  time  scale  of  order  of  the  internal 
relaxation  times,  as  discussed  below. 
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For  the  calculation  of  gas  flows  where  an  appreciable  dissociation  occurs,  the 
temperature  changes  considerably,  due  to  the  energy  required  for  the  dissociation,  so 
that  the  rate  constants  change  during  the  dissociation.  For  a  tractable  analysis  of 
aerodynamic  problems  without  use  of  large  machine  programs,  some  simplification  of 
Equation  (3-31)  is  required.  These  can  be  obtained  by  applying  the  approximations 
Introduced  by  Lighthill  for  the  ideal  dissociating  gas. 


Ideal  dissociating  g as 

The  approximations  for  an  ideal  dissociating  gas,  discussed  in  Section  2,  were 
formulated  by  Freeman75  to  describe  the  kinetics  of  dissociation.  He  included  the 
approximation  to  the  rate  constants  given  by  Equations  (3-28)  and  (3-29),  and  further 
neglected  the  change  in  dissociation  rate  associated  with  the  different  rate  constants 
for  atomic  and  molecular  collisions.  This  latter  simplification  is  obviously  a  poor  one 
if  the  degree  of  dissociation  is  large.  Changing  variables  to  the  mass  fraction  of 
dissociated  atoms 

a  =  yAMA  ,  (3-32) 

where  MA  is  the  atomic  weight  of  A  ,  the  conservation  of  mass  gives 

rA2(2MA)  =  (1-0.)  . 

The  approximations  then  reduce  Equation  (3-31)  to 


da 

dt 


C'T  "S/0 


(l-a)e 


-D/R0T 


(3-33) 


where  pD  is  a  constant  that  enters  from  the  Lighthill  approximation  to  the  partition 
functions  and  C'  and  s  are  related  to  Cf  and  77f  of  Equation  (3-29).  This 
formulation  was  used  by  Freeman75  to  discuss  dissociation  in  shock  waves,  and  by  Bray76 
to  discuss  recombination  in  expanding  flows.  It  has  since  been  used  by  several  authors 
to  investigate  the  magnitude  of  dissociation  effects  in  diatomic  gas  flows,  but  is  not 
suitable  for  careful  interpretation  of  experiments. 


Vibrat ional-dissociat ion  coupling 

The  dissociation  of  a  diatomic  molecule  by  collision  with  an  inert  atom  provides 
the  simplest  dissociation  reaction,  both  from  a  theoretical  and  experimental  point  of 
view.  For  example,  a  small  amount  of  oxygen  in  a  large  amount  of  argon  can  be  studied 
in  a  shock  tube  over  a  large  temperature  range77.  Theoretically,  there  is  no  need  to 
consider  temperature  changes  during  the  reaction,  and  vibrational -exchange  reactions 
and  dissociation  by  impact  with  the  small  amounts  of  02  or  0  in  the  experiment  are 
negligible.  Then  in  the  analysis  of  this  reaction  it  is  feasible  to  consider  the 
details  of  the  vibrational  excitation  of  the  molecule  as  the  intermediate  process  by 
which  the  molecule  proceeds  to  dissociation.  This  effect  is  of  special  importance  in 
applying  recombination  rates  to  problems  of  expanding  flows,  where  the  recombination 
rate  to  be  used  is  determined  from  the  measured  dissociation  rate  and  the  equilibrium 
constant,  through  Equation  (3-27). 


For  the  dissociation  of  a  small  amount  of  02  in  Ar  ,  Equation  (3-26)  is 


1  ^ 

p  dt 


-  kF>b2 ?Ar  +  kR?0?Ar  • 


(3-34) 
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However,  if  each  vibrational  level  of  the  molecule  is  considered  as  a  separate  species, 
with  concentration  y0  v  in  the  v  vibrational  level,  then  a  chemical  equation  of 
this  type  can  be  written  for  each  level,  along  with  terms  to  describe  the  change  in 
population  of  the  v  level  due  to  transitions  from  other  vibrational  levels.  The 
resulting  equations  are 


1  <fro2v 

p  dt 


®v+  1,  v^Ar 


70 
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■(ev  +  f  ev)/kT 
e  V  +  1  V  7o2v 
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-  e 


-(fv- Vi)/kT 


^o2,  v-l 


kPv70  v7Ar  +  kRv7o7Ar’  v  =  0,  1,2...  (3-35) 


where  kpv  and  kRv  are  the  forward  and  reverse  rates  for  dissociation  from  the 
vibrational  level  v  ,  and  the  remaining  terms  are  the  same  as  given  in  Equation  (3-12) 
for  vibrational  relaxation. 


These  equations  contain  two  approximations  that  limit  their  validity.  One  is  that 
the  translational  and  rotational  distributions  are  in  Boltzmann  distribution  at  tempera¬ 
ture  T  ,  and  are  unaffected  by  the  dissociation  process,  leading  to  rate  constants  in 
Equation  (3-35)  that  depend  only  on  T  .  In  particular,  the  forward  and  reverse  rates 
for  dissociation  from  the  v  level  are  taken  as  related  through  detailed  balance  as 


kFV 


S«vibe 


ey/kT 


(3-36) 


The  second  is  that  only  one-quantum  vibrational  transitions  are  included,  whereas 
multiquantum  transitions  certainly  occur  at  higher  vibrational  levels.  These  approxi¬ 
mations  are  included  to  simplify  the  description,  and  with  the  idea  that  to  some  extent 
the  values  of  the  rate  constants  could  be  modified  to  compensate  for  these  effects. 


The  equation  to  replace  (3-34)  is  obtained  by  summing  Equation  (3-35)  over  all 
vibrational  states  v  .  The  sums  for  internal  excitation  are  then  equal  to  zero,  and 
there  results 


P  dt 


-7Ar  SkFv7o2v  +  7Ar7| 


Ar  /0 


SIk 

v 


Rv 


(3-37) 


If  a  Boltzmann  distribution  at  temperature  T  is  assumed  for  the  vibrational  distribu¬ 
tion,  Equation  (3-34)  is  obtained  directly.  This  distribution  is  not  correct,  however, 
because  the  finite  rate  of  vibrational  relaxation  and  the  dissociation  process  both 
affect  the  distribution,  which  must  be  obtained  by  solving  Equation  (3-35). 


There  are  three  separate  regions  of  time  that  can  be  discussed  in  terms  of  Equation 
(3-35).  The  first  is  during  the  early  portion  of  the  dissociation  (or  recombination) 
process,  at  times  less  than  the  vibrational  relaxation  time.  For  dissociation 
immediately  behind  a  shock  wave,  there  are  essentially  no  molecules  in  upper  vibrational 
states  (from  which  most  dissociation  normally  occurs)  so  that  essentially  no  dissociation 
can  occur  until  a  time  of  order  of  the  vibrational  relaxation  time  (the  induction  time). 
Somewhat  after  this  time  the  collisional  excitation  of  vibrational  states  reaches  a 
steady  state  with  dissociational  destruction  of  these  states,  and  a  distribution  is 
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established  which  persists  as  long  as  recombination  C^Ar^o *s  unin|Portai,t.  In 
this  steady-state  distribution  the  upper  vibrational  levels  are  underpopulated,  so 
that  the  dissociation  rate  is  suppressed.  Finally,  when  enough  atoms  are  produced  to 
make  recombination  important,  the  vibrational  populations  approach  their  final,  Boltzmann 
distribution.  A  similar  sequence  can  be  described  for  a  recombination  experiment. 
Consider  an  experiment78  where  an  ozone-argon  mixture  is  shock-heated  to  moderate 
temperature.  The  ozone  quickly  dissociates  to  02  and  0  .  Most  of  the  0  atoms 
must  then  recombine  to  form  the  final  equilibrium  mixture.  At  very  short  times  the 
atoms  recombine  quickly  to  populate  the  upper  vibrational  states.  At  longer  times 
the  steady  state  is  reached  where  upper  vibrational  levels  are  overpopulated,  and 
recombination  to  upper  states,  vibrational  deexcitation  of  these  states,  and  thermal 
redissociation  of  these  states  are  balanced.  The  net  effect  (due  to  the  redissociation) 
is  a  depression  of  the  observed  recombination  rate.  The  third  region  is  again  the 
final,  Boltzmann  state. 

One  question  that  can  then  be  asked  of  Equation  (3-35)  is  under  what  conditions  will 
the  steady-state  dissociation  rate  and  the  steady-state  recombination  rate  be  related 
by  the  equilibrium  constant  as  in  Equation  (3-27).  This  question  has  been  addressed  in 
many  theoretical  papers  7  9  ,  8  4.  It  has  been  demonstrated  by  Snider80,  in  agreement  with 
Rice’s  original  discussion79,  that  the  forward  and  reverse  rates  measured  within  the 
steady-state  region  are  depressed  by  the  same  factor,  so  that  their  ratio  is  still  equal 
to  the  equilibrium  constant.  Detailed  consideration84  of  the  region  of  applicability 
of  this  statement  shows  that  it  applies  to  the  conditions  of  the  02  -  Ar  shocktube 
experiments.  Other  types  of  reactions,  where  Equation  (3-27)  can  be  applied,  are 
discussed  in  Reference  83,  and  it  is  clear  from  that  discussion  that  considerable 
caution  must  be  used  in  relating  forward  and  reverse  rate  constants  for  more  complicated 
reactions.  Very  few  experimental  results  are  available  where  both  the  forward  and 
reverse  rates  have  been  measured. 

The  complete  solution  of  Equations  (3-35)  requires  a  statement  of  all  the  values  of 
the  rate  constants  (Qv  +  1  v  and  kpv).  This  has  been  accomplished  for  the  steady-state 
region  using  various  approximations  for  these  rates88’86.  In  Reference  85  both  classical 
and  quantum  mechanical  models  are  investigated,  and  the  equations  for  the  classical 
model,  similar  to  Equation  (3-35),  are  reduced  to  a  diffusion-equation  form.  The 
calculation  of  rate  constants  from  a  classical  model  is  discussed  in  Reference  87  and 
in  other  papers  referenced  there. 

The  induction  time  for  02  dissociation  behind  strong  shocks  in  02  -  Ar  mixtures 
has  been  measured77  over  the  temperature  range  of  5,  500-18, 000°K.  The  measurements 
show  that  the  induction  time  is  of  the  order  of  the  vibrational  relaxation  time,  and 
appears  to  increase  more  rapidly  than  the  vibrational  relaxation  time  at  lower  tempera¬ 
tures.  The  experimental  method  (measurement  of  02  ultraviolet  absorption)  did  not 
permit  measurements  at  temperatures  below  5500°K.  However,  experiments  employing  the 
resonance -absorption  of  atomic  oxygen88  are  now  being  performed  to  extend  this  informa¬ 
tion  to  the  lower-temperature  range.  The  theoretical  problem  of  relating  the  induction 
time  to  the  relaxation  rates  is  given  in  Reference  89. 

In  shock-waves  at  very  high  temperature  the  times  for  vibrational  and  dissociational 
relaxation  become  comparable,  and  the  separation  into  specific  time  regions  is  no  longer 
possible.  High-temperature  experiments  where  these  regions  are  obviously  mixed  are 
described  in  References  90  and  91.  There  are  no  theoretical  treatments  available  for 
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flows  under  these  conditions,  but  basic  work  on  the  excitation  and  dissociation  from 
high-energy  collisions  is  a  subject  of  continuing  study92’93. 


Normal- shock  coupled  vibration-dissociation 


For  treatment  of  shock  calculations  fo”  conditions  where  vibration  and  dissociation 
are  coupled,  a  simplified  interpretation  of  Equation  (3-37)  was  given  in  References 
94-96.  It  was  assumed  that  the  solution  of  Equation  (3-35)  for  yQ  v  could  be  repre¬ 
sented  by  a  Boltzmann-like  distribution  at  an  unknown  temperature  Ty  .  Hie  value  of 
Ty  was  determined  from  a  vibrational  energy  relaxation  equation  which  included 
terms95’96  to  account  for  the  loss  of  vibrational  energy  due  to  dissociation.  Hie 
reverse  rate  constant  for  recombination  to  the  vibrational  level  was  taken  to  have  the 
form96 
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where  U  was  an  adjustable  parameter.  The  reverse  rate  was  related  to  the  forward 
rate  through  Equation  (3-36).  Hie  parameter  kU  was  taken  as  1/6  of  the  dissociation 
energy.  Calculations  with  this  model96  succeeded  in  duplicating  the  induction  time  of 
the  02  -  Ar  experiments77  and  resulted  in  a  steady-state  vibrational  temperature 
considerably  below  the  translational  temperature,  thus  depressing  the  dissociation  rate 
constant,  as  would  be  expected.  However,  the  assumption  concerning  the  vibrational 
temperature  cannot  be  justified,  and  numerical  calculations97 ■ 9B  show  that  the  upper 
levels  are  more  severely  depleted  than  this  assumption  would  imply. 


3.2.3  Other  Chemical  Reactions 

The  shock  tube  has  now  become  a  standard  method  of  studying  high-temperature  chemical 
reactions  and  has  been  applied  to  a  wide  variety  of  problems.  This  research  during  the 
past  several  years  has  been  comprehensively  reviewed  by  Bauer12,  and  the  wide  range  of 
chemical  kinetic  studies  being  pursued  with  the  shock  tube  is  clearly  demonstrated.  A 
review  of  the  very  recent  chemical  kinetic  work  is  given  in  Reference  13.  The  many 
reactions  of  both  clean-air  and  carbon  chemistry  involved  with  reentry  ablation  problems 
are  reviewed  by  Heicklen  in  Reference  14.  The  remarks  presented  in  the  present  section 
will  be  concerned  specifically  with  the  nitric  oxide  reactions  in  air,  in  order  to 
illustrate  the  level  of  complication  for  a  two-gas  mixture. 

The  primary  chemical  processes  in  shock-heated  air  other  than  the  dissociation 
reactions,  are  those  resulting  in  production  of  nitric  oxide.  They  were  studied  by 
Glick,  Klein  and  Squire99  in  the  chemical  shock  tube,  where  air  was  heated  to  high 
temperature  by  a  shock  wave  and  then  cooled  rapidly  after  a  fixed  length  of  time  by 
an  expansion  wave.  Subsequent  analysis  of  the  gas  in  the  shock  tube  permitted  evalua¬ 
tion  of  the  mechanisms  and  rate  constant.  It  was  concluded  that  the  mechanism  was  the 
one  previously  proposed  by  Zel’  dovich 10°,  namely 

(1)  02  +  M  j=!  20  +  M 

(2)  0+N2  zz  NO  +  N 

(3)  N  +  02  s  NO  +  0  . 

Thus  the  first  oxygen  atoms  are  produced101  by  dissociation  of  02  ,  and  then  two 
“shuffle”  reactions  become  active  in  producing  NO  .  The  importance  of  these  shuffle 


32 


reactions  in  expanding- flow  chemistry  is  discussed  in  Section  4.3.  In  shock-heated 
flows  the  high  speed  of  reactions  (2)  and  (3)  results  in  a  large  production  of  NO  , 
in  an  attempt  to  bring  the  NO  into  local  equilibrium  with  the  0?  and  N2  concentra¬ 
tions.  When  the  02  later  dissociates,  the  NO  concentration  decreases.  Such  NO 
overshoots  were  predicted  by  Duff  and  Davidson102,  and  by  Lin  and  Teare103,  and  have 
been  observed  in  shock-tube  flows, 

Recent  experiments 10  4  indicate  that  it  may  be  possible  to  investigate  the  NO-produ- 
cing  reactions  in  as  much  detail  as  has  been  obtained  for  dissociation  reactions.  Camac 
and  Feinberg104  have  observed  induction  times  in  the  formation  of  nitric  oxide  in  air, 
and  have  correlated  this  with  N2  vibrational  relaxation  times.  They  suggest  that  the 
vibrational  excitation  of  the  nitrogen  may  be  a  prerequisite  for  reaction  (2).  Again, 
sensitive  experiments  to  determine  atomic  oxygen  production88  would  add  considerably  to 
the  total  knowledge  about  this  reaction. 

The  formation  of  NO  from  the  reaction 

(4)  N2  +  02  -  2NO 

was  not  observed  in  the  work  reported  in  Reference  104,  leading  to  the  conclusion  that 
this  rate  is  at  least  10  times  slower  than  previously  reported.  The  reverse  reaction 
had  previously  been  measured105,15,  using  NO-absorbtivity  measurements  to  determine 
the  rate  at  which  NO  was  destroyed  behind  a  shock  wave  in  pure  NO  .  In  Reference  104 
it  was  determined  that  this  rapid  NO  destruction  was  caused  by  the  reaction 

(5)  2N0  -  N20  +  0 

rather  than  the  inverse  of  reaction  (4).  They  conclude  that  reaction  (4)  is  usually 
unimportant  for  air  kinetics.  This  simplification,  along  with  chemical  experiments  in 
expanding  flows  (Section  4)  may  permit  more  detailed  understanding  of  the  shuffle 
reactions. 

These  recent  considerations  of  the  nitric  oxide  reactions  constitute  the  major 
changes  in  the  knowledge  of  measured  reaction  rates  for  neutral  pure-air  chemistry 
since  the  review  presented  by  Wray15  in  1962.  The  more  recent  data  on  N2  dissocia¬ 
tion  rates  should  also  be  noted74. 

3.2.4  Ionization 

The  problems  of  ionizing  and  deionizing  reactions  in  air  have  received  a  considerable 
amount  of  attention  because  of  their  connection  with  upper-atmospheric  phenomena  and 
with  nuclear  reactions,  as  well  as  their  importance  in  reentry  phenomena.  The  well 
known  communications  black-out  during  reentry  results  from  a  high  density  of  free 
electrons  in  the  vehicle  flowfield,  and  free  electrons  in  the  trail  of  reentering 
missiles  provide  the  reflecting  material  for  radar  returns.  Thus,  there  has  been  a 
broad  interest  in  the  identification  of  ionic  reactions  and  the  measurement  of  rate 
constants  over  a  range  of  temperatures  from  atmospheric  up  to  many  times  that  associated 
with  reentry.  The  massive  literature  available  on  this  subject  is  under  constant  review. 
There  is  presented  here  a  brief  summary  of  the  review  articles  that  are  currently  avail¬ 
able,  with  emphasis  on  those  directed  toward  aerodynamic  problems. 
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The  simplest  models  of  high-temperature  air  flows  use  a  single  ionizing  reaction1,15: 

N  +  0  s  N0+  +  e  . 

For  a  good  description  of  atmospheric  reentry,  however,  a  host  of  other  reactions  must 
be  considered.  For  example,  the  forebody  flow  of  an  Apollo  vehicle  can  be  divided  into 
four  separate  regions106,  with  each  region  corresponding  to  a  different  type  of  ionizing 
reaction.  For  each  type  of  reaction,  many  separate  reactions  must  be  considered.  In 
References  16  and  106  this  problem  is  reviewed  for  Apollo-reentry  flows  for  pure-air 
reactions. 

For  missile  reentry,  the  review  of  gas  phase  chemistry  by  Heicklen14  includes  some 
discussion  of  ion -producing  reactions.  The  subject  of  pure-air  ion  chemistry  is 
reviewed  within  the  context  of  wake  flows  by  Sutton107  and  the  requirement  of  properly 
identifying  the  dominant  ionic  reactions  is  clearly  demonstrated. 

General  reviews  of  both  theory  and  experiment  are  given  in  several  texts9"17,18  and 
special  publications19.  Experimental  data  for  electron-impact  ionization  was  reviewed 
by  Kieffer  and  Dunn20  and  a  comprehensive  bibliography  of  both  theoretical  and  experi¬ 
mental  results  was  recently  prepared  by  Kieffer21.  A  current  bibliography  of  electron 
swarm  data  is  also  available22.  The  recent  publication  of  the  NASA  Reaction  Rate 
Handbook23  provides  a  summary  with  emphasis  on  the  ion  chemistry  of  high-altitude 
nuclear  bursts. 
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SECTION  4 

NONEQUILIBRIUM  FLOWS  IN  SUPERSONIC  NOZZLES 


The  classical  convergent-divergent  de  Laval  nozzle  accelerates  gas  from  subsonic 
to  supersonic  velocities  with  attendant  conversion  of  internal  or  thermal  energy  to 
the  kinetic  energy  of  bulk  motion.  The  isentropic  expansion  of  a  perfect  gas  in  such 
a  nozzle  is  a  familiar  flow  example  discussed  in  all  textbooks  on  gasdynamics  (e. g. 
Liepmann  and  Roshko1).  As  the  usual  analysis  based  on  pseudo-one-dimensional,  inviscid 
flow  shows,  local  sonic  or  critical  flow  conditions  are  attained  at  the  minimum  cross- 
section  area,  or  throat,  of  the  nozzle.  Downstream  of  the  throat  the  gas  pressure, 
density,  and  temperature  decrease  very  rapidly,  while  the  gas  velocity  increases; 
typical  maximum  cooling  rates  might  be  in  the  range  105  to  106  °K/sec.  For  large 
expansions  (large  ratio  of  nozzle  exit  area  to  throat  area)  the  exit  flow  Mach  number 
becomes  hypersonic,  and  the  exit  values  of  p  ,  p  ,  and  T  are  very  much  less  than 
the  pre-expansion  nozzle  reservoir  values. 

The  typical  situation  of  present  interest  is  where  the  inlet  or  reservoir  gas  is 
in  thermal-chemical  equilibrium  prior  to  expansion.  It  turns  out  that  over  a  wide 
range  of  reservoir  conditions  and  nozzle  scale  of  practical  interest,  the  rate  or 
relaxation  processes  discussed  in  Section  3  have  characteristic  times  of  the  same 
order  of  magnitude  as,  or  larger  than,  particle  residence  or  flow  times  in  the  nozzle. 

As  a  consequence,  departures  from  equilibrium  in  practice  can  be  quite  large.  If 
sufficient  energy  is  involved,  the  relaxing  mode  can  interact  strongly  with  the  local 
thermo-gasdynamic  state  of  the  nozzle  flow.  Typically,  energy  is  withheld  or  "frozen” 
in  the  relaxing  mode  during  expansion  (e.  g.  lagging  atom  recombination)  with  a  con¬ 
sequent  reduction  of  nozzle  flow  pressure,  velocity,  and  temperature  below  the  corres¬ 
ponding  equilibrium  flow  values. 

The  discussion  will  be  limited  to  nonequilibrium  effects  in  such  nozzle  flows  in¬ 
volving  either  the  vibrational,  chemical,  or  ionization  relaxation  processes  discussed 
in  Section  3.  While  departures  from  local  equilibrium  in  translational  as  well  as 
rotational  degrees  of  freedom  have  been  observed  in  steady  expansion  flows13'16  signifi¬ 
cant  nonequillbrium  effects  in  these  degrees  mainly  occur  at  lower  densities  than  of 
concern  here. 

The  pseudo-one-dimensional  formulation  of  the  inviscid,  nonequilibrium  nozzle-flow 
problem  will  be  briefly  outlined  first,  and  some  general  features  of  the  solution  and 
aspects  of  numerical  integration  of  the  governing  equations  will  be  discussed. 

Subsequent  sections  then  discuss  some  particular  theoretical  and  experimental  results 
for  nonequilibrium  nozzle  flows  involving  dissociation,  vibration,  and  reacting  gas 
mixtures.  Previous  review  articles  pertinent  to  the  scope  of  the  present  discussion 
include  the  excellent  review  papers  by  Bray8;  Penner,  Porter,  and  Kushida9;  and 
Olson10'  ll. 
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4.  1  Pseudo-One-Dimensional  Formulation  and 
Methods  of  Solution 

Most  analyses  of  nonequilibrium  effects  in  steady  supersonic  nozzle  flows  have  been 
simplified  by  assuming  the  flow  to  be  pseudo-one-dimensional  as  well  as  inviscid, 
adiabatic,  and  shock  free.  The  one-dimensional  flow  approximation  assumes  that  all 
flow  quantities  depend  only  on  axial  distance  x  along  the  nozzle;  gradients  and 
velocity  components  normal  to  x  are  neglected.  This  approximation  is  valid  provided 
the  nozzle  expansion  angle  is  sufficiently  small.  The  pseudo-one- dimensional  flow 
approximation  is  certainly  adequate  for  revealing  basic  effects,  and  is  not  inconsistent 
with  the  limited  state  of  knowledge  concerning  typical  nonequilibrium  rate  processes 
of  interest.  Even  disregarding  two-dimensional  effects  and  considering  only  the 
simplest  rate  processes,  such  as  dissociation-recombination  of  a  diatomic  gas,  analysis 
of  the  problem  is  still  far  from  trivial  . 

From  pseudo-one  dimensional  flow,  the  fluid-dynamic  equations  for  inviscid  steady 
flow  simplify  to 


puA  =  m  =  constant 


(4-1) 


du  dp 

pu  —  +  —  =  0  (4-2) 

dx  dx 


h 


constant  . 


(4-3) 


These  equations  express  conservation  of  mass,  momentum,  and  energy,  respectively.  The 
quantity  A  is  the  local  cross-section  area  of  the  nozzle,  assumed  to  be  a  prescribed 
function  of  distance  x  . 


If  we  regard  p  ,  p  ,  T  ,  and  u  as  the  primary  thermodynamic  dependent  variables 
of  interest,  then  three  additional  relations  are  needed  to  complete  the  formulation. 

These  are:  (i)  the  thermodynamic  equation  of  state  for  ideal-gas  mixtures,  i.e. 
s 

p  =  pRnT  5Z  i  (ii)  a  caloric  equation  of  state  expressing  the  specific  enthalpy 
j  =  i  J 

h  in  terms  of  the  translational  temperature  T  and  the  nonequilibrium  variables 
considered,  e.g.  vibrational  energy  eyj  or  species  mass  concentrations  y^  ;  and 
(iii)  the  appropriate  system  of  rate  equations  and  element  conservation  equations 
expressing  the  rates  of  change  of  the  nonequi librium  variables.  The  particular  expres¬ 
sions  for  h  and  typical  rate  equations  fer  vibrational  or  chemical  nonequilibrium 
have  been  given  in  Sections  2  and  3. 

The  necessary  boundary  conditions  can  take  various  forms  depending  on  the  particular 
problem  or  nozzle  geometry  considered.  Fo-  initially  supersonic  flow  expanding  in  a 
continuously  diverging  nozzle,  the  initial  state  may  be  specified.  The  case  of  super¬ 
critical  ilow  through  a  converging-diverging  nozzle  is  somewhat  more  difficult  as 
regards  boundary  conditions.  Commonly  in  this  case  the  thermodynamic  state  of  the 

•  Consideration  of  two-dimensional  effects  in  nonequilibrium  nozzle  flows  may  be  found  in  t 

References  119-122. 
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gas  effectively  far  upstream  of  the  nozzle  throat  is  specified,  i.e.  the  equilibrium 
reservoir  or  stagnation  conditions,  denoted  by  p 0  ,  P0  ,  and  T0  ,  say,  corresponding 
to  the  gas  initially  at  rest.  In  addition,  for  supercritical  flow  a  boundary  condition 
must  be  applied  at  the  throat  to  fix  or  determine  the  nozzle  mass  flow  m  (Equation 
(4-1))  as  the  maximum  or  critical  value.  The  problem  is  thus  a  two-point  boundary 
value  one.  In  general,  the  throat  boundary  condition  equates  the  local  flow  velocity 
to  the  appropriate  local  sound  speed,  as  further  discussed  below. 

Because  actual  reaction  rates  are  finite  and  streamwiso  gradients  exist,  the 
accelerating  flow  from  an  equilibrium  reservoir  must  progressively  depart  from  equili¬ 
brium,  as  is  evident  from  the  rate  equations  of  Section  3.  As  previously  mentioned, 
the  effect  of  nonequilibrium  is  to  reduce  the  pressure,  temperature,  and  velocity  from 
equilibrium  values,  whereas  the  local  density  is  increased.  Because  of  the  normal 
decrease  of  gas  density  and  temperature  with  flow  acceleration,  the  local  reaction 
rates  typically  decrease  strongly  (or  local  relaxation  lengths  increase  strongly)  as 
the  gas  expands  down  the  nozzle.  This  quenching  effect  becomes  particularly  strong 
immediately  downstream  of  the  throat  where  streamwise  gradients  are  large.  As  a 
consequence  of  departure  from  equilibrium  together  with  quenching,  the  expanding  flow 
typically  undergoes  a  transition  from  equilibrium  to  frozen  conditions.  The  location 
of  this  transition  in  the  nozzle  depends  on  nozzle  geometry  and  reservoir  pressure, 
as  well  as  on  other  parameters.  Of  course,  the  maintaining  of  equilibrium  flow  is 
favored  by  large  scale  and  high  pressures.  Over  a  wide  range  of  nozzle  scale  and 
reservoir  pressures  of  practical  interest,  significant  departures  from  equilibrium 
flow  occur  only  downstream  of  the  nozzle  throat. 

It  is  often  useful  to  have  the  limiting  solutions  for  equilibrium  and  frozen  nozzle 
flows,  corresponding  to  infinite  and  zero  reaction  rates,  respectively,  for  comparison 
with  experiment  or  corresponding  nonequilibrium  (finite- rate)  solutions.  Hie  calcula¬ 
tion  of  these  two  limiting  cases  is  relatively  straightforward  as  both  limiting  flows 
are  isentropic4,  and  their  solution  involves  only  algebraic  equations.  Thus  Equation 
(4-2)  may  be  replaced  in  both  cases  by  the  appropriate  entropy  relation  S(p,  T,7j)  = 

SQ  =  constant  .  The  rate  equations  may  be  replaced  by  the  corresponding  equilibrium 
Law  of  Mass  Action  in  the  equilibrium  case,  and  by  the  conditions  of  constant  vibra¬ 
tional  energies  or  species  mass  concentrations  in  the  frozen  case.  The  throat  boundary 
condition  to  determine  the  critical  mass  flow  equates  the  flow  velocity  at  the  geometric 
throat  to  the  local  equilibrium  sound  speed  ag*  for  equilibrium  flow  and  to  the 
local  frozen  sound  speed  ap*  for  frozen  flow.  Because  of  the  limiting  nature  of 
both  flows,  i.e.  infinite  or  zero  reaction  rates,  both  solutions  depend  only  on  the 
local  nozzle  area  ratio  A/A^  (A^  =  throat  area)  for  given  reservoir  conditions,  and 
not  on  the  nozzle  shape  or  scale.  Whereas  the  (completely)  frozen  solution  is  further 
independent  of  the  reservoir  state  (e.g.  when  normalized  as  p/p0,  etc.),  this  is  not 
so  for  the  equilibrium  solution.  Detailed  discussion  of  a  calculation  method  for 
equilibrium  and  frozen  nozzle  flows  is  given  by  Vincenti  and  Kruger4.  Specific  solutions 
for  such  nozzle  flows  of  air  are  given  in  References  94-96. 

In  contrast  to  the  equilibrium  or  frozen  flow  cases  where  the  solution  at  any  point 
depends  only  on  the  overall  area  change  involved,  the  solution  for  nonequilibrium  or 
finite-rate  nozzle  flow  depends  at  any  point  on  the  entire  history  of  the  flow  up  to 
that  point.  In  particular,  the  nonequilibrium  flow  solution  at  any  point  depends  on 
the  preceding  history  of  the  nozzle  shape  as  well  as  on  the  area  change  itself. 

Typically  in  the  nonequilibrium  case  the  governing  differential  equations,  i.e.  the 
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fluid  dynamic  and  coupled  nonlinear  rate  equations,  are  intractable  analytically,  and 
a  solution  must  therefore  be  obtained  by  numerical  or  step-by-step  forward  integration. 
Beginning  with  the  development  of  exact  numerical  methods  and  solutions  for  nonequili¬ 
brium  nozzle  flows  of  diatomic  gases17' 18,21,  much  effort  has  been  devoted  to  the 
successful  development  of  high-speed  computing  machine  methods  for  calculation  of  non¬ 
equilibrium  nozzle  flows  involving  complex  systems  of  coupled  reactions67'79. 

The  numerical  integration  of  the  nonequilibrium  nozzle  flow  equations  encounters  a 
basic  difficulty  when  the  starting  point  is  an  equilibrium  state  upstream  of  the  throat 
and  the  nozzle  geometry  is  specified.  This  difficulty  is  the  resulting  two-point 
boundary-value  problem.  The  critical  mass  flow  m  ,  which  is  determined  by  a  boundary 
condition  at  the  throat  region  equating  the  local  flow  velocity  to  the  local  frozen 
sound  speed,  cannot  be  determined  a  priori  since  the  flow  in  the  throat  region  depends 
on  the  entire  preceding  flow  history.  A  further  complication  is  the  fact  that  under 
nonequilibrium  conditions  the  sonic  condition  u  =  aF  for  supercritical  flow  does  not 
in  general  occur  at  the  geometric  nozzle  throat  but  an  unknown  distance  downstream  of 
that  point. 

To  illustrate  the  foregoing  remarks,  consider  the  familiar  relation 

1  du  (l/A)dA/dx 

=  - 2 - —  .  (4-4) 

u  dx  u 

(dp/dx)/(dp/dx) 

which  is  readily  derived  from  Equations  (4-1)  and  (4-2).  As  for  the  case  of  perfect 
gas  flow,  this  relation  states  that  at  the  geometric  throat.  (dA/dx)  =  0  ,  the  flow 
velocity  must  equal  [(dp/dx)/(dp/dx)] i  if  the  flow  is  to  be  continuously  accelerating, 
i.e.  du/dx  i  0  .  For  the  equilibrium  or  frozen  flow  limits,  we  have  [(dp/dx)/(dp/dx)] 
equal  to  ag  or  aF  ,  respectively.  However,  for  nonequilibrium  flow  the  quantity 
(dp/dx)/(dp/dx)  is  not  particularly  useful  since  it  is  no  longer  a  function  of  the 
thermodynamic  state  only,  but  depends  on  the  particular  flow  history  involved.  An 
alternative  relation  to  Equation  (4-4)  may  be  obtained  from  Equations  (4-1)  to  (4-3) 
by  considering  the  enthalpy  h  as  a  function  of  p  ,  p  ,  and  a  ,  say,  considering 
the  example  of  a  dissociating  diatomic  gas.  This  relation  is“ 

1  dA  ha  da 
1  du  A  dx  ph„  dx 

-  —  -  - 5 - * - ,  (4-5) 

u  dx  u 

H?  ~ 1 

aF 


(9h/3a). 


etc. ,  and  aF  is  the  frozen  sound  speed  defined  by 


In  general,  the  quantity  (ha/ph^) (da/dx)  in  Equation  (4-5)  is  positive.  Thus  the 
numerator  of  the  right  side  will  vanish  only  downstream  of  the  throat  where  dA/dx 
is  positive.  At  that  point  then,  the  condition  u  =  ap  must  simultaneously  apply 
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(to  make  the  denominator  vanish)  if  the  flow  is  to  be  supercritical.  In  practice,  it 
is  usually  found  that  the  point  where  u  =  is  displaced  only  a  small  distance 
downstream  of  the  geometric  throat. 

The  simplest  solution  of  the  difficulty  of  the  unknown  critical  mass  flow  is  to 
determine  its  value  by  trial  and  error.  That  is,  the  solution  is  computed  to  the 
throat  using  successive  guesses  for  the  mass  flow  until  supercritical  flow  downstream 
of  the  throat  is  obtained.  The  initial  guess  is  guided  by  the  equilibrium  and  frozen 
limits  which  do  not  differ  greatly  and  which  are  less  and  greater,  respectively,  than 
the  desired  value.  This  approach  has  been  used  by  Bray17,  and  by  Hall  and  Russo18, 
and  a  systematic  refinement  is  discussed  by  Emanuel70.  Another  solution  to  the  diffi¬ 
culty,  described  by  Eschenroeder,  Boyer,  and  Hall67,  involves  an  inverse  method.  In 
this  approach  the  starting  point  is  the  density  distribution  corresponding  to  equilibrium 
flow  in  the  subsonic  region  of  the  given  nozzle.  A  perturbation  procedure,  including 
perturbation  of  the  nozzle  geometry,  then  provides  an  exact  nonequilibrium  flow  for  a 
corresponding  new  nozzle  geometry  (subsonic)  matching  this  density  distribution.  Because 
of  the  general  insensitivity  of  flow  density  to  nonequilibrium,  the  new  subsonic 
nozzle  geometry  closely  matches  that  originally  given. 

A  second  difficulty  encountered  in  numerically  solving  the  nonequilibrium  nozzle 
flow  is  that  of  starting  the  integration  of  the  rate  equations  from  an  initial  equili¬ 
brium  state  wherein  the  initial  gradients  of  7j  and  eyj  given  by  the  rate  equations 
are  identically  zero.  As  pointed  out  by  Bloom  and  Ting9,  the  initial  departure  from 
equilibrium  is  mathematically  a  singular  perturbation  behavior,  analogous,  for  example, 
to  the  fluid  dynamic  boundary  layer.  Starting  methods  used  to  overcome  this  difficulty 
include  asymptotic  expansions17,18,  linear  perturbation  about  the  corresponding  equili¬ 
brium  flow67,71,  and  direct  starting  from  equilibrium  with  the  nonlinear  rate  equations68. 

Apart  from  the  special  problem  in  starting  the  integration  from  an  initial  equili¬ 
brium  state,  it  may  be  noted  that  the  usual  numerical  methods  for  integration  of  the 
exact  nonlinear  rate  equations,  such  as  the  fourth-order  Runge-Kutta  scheme,  for 
example,  generally  encounter  a  serious  instability  problem  for  near  equilibrium  flow 
conditions.  In  order  to  maintain  stability  near  equilibrium,  the  integration  step 
size  must  be  kept  very  small.  Thus  the  overall  computing  time  can  become  unreasonably 
long  if  an  appreciable  region  of  near  equilibrium  flow  is  involved.  The  underlying 
reasons  for  this  behavior  are  discussed  in  Reference  69.  In  the  past  few  years, 
methods  for  overcoming  this  problem  have  been  developed  by  several  authors73"77. 

In  spite  of  these  difficulties,  machine  methods  for  computing  nonequilibrium  nozzle 
flows  with  complex  systems  of  coupled  reactions  have  rapidly  been  brought  to  a  state 
of  considerable  efficiency  and  precision.  Unfortunately,  knowledge  of  pertinent 
reaction  mechanisms  and  rates  has  not  advanced  at  a  similar  pace.  Currently,  un¬ 
certainties  in  the  fundamental  kinetics  often  place  a  serious  limitation  on  the  signi¬ 
ficance  of  numerical  solutions. 


4.2  Dissociation  Nonequilibrium  in  Nozzle  Flows 
of  Diatomic  Gases 

This  section  discusses  an  important  example  of  chemical  nonequilibrium  in  steady 
nozzle  flow  involving  a  single  dominant  reaction,  namely,  that  of  symmetrical  diatomic 
gas  which  undergoes  dissociation  and  recombination  by  two-  and  three-body  collision 
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processes,  respectively,  during  the  flow  expansion.  The  appropriate  type  of  rate 
equation  for  this  example  can  be  taken  as  that  for  the  Lighthill  ideal  dissociating 
gas  discussed  in  Section  3,  i.e. 


Da 

Dt 


da 

u  —  =  C'pT 
dx 


(1  -  a)e 


-D/RjjT 


P  2 
—  a2 

Pd 


(4-6) 


In  terms  of  a  local  equilibrium  value  of  a  equal  to  ae  ,  defined  in  terms  of  density 
and  temperature  by  the  relation 


ae  _  -D/V 

1  -  ae  P 

Equation  (4-6)  can  be  written 

da  c'p2T's 

u  —  =  - 

dx  Pd 


l  -  a 
l  -  a„ 


ae  " 


For  a  «  1  ,  this  becomes 


da 


u 


dx 


C'p2  T~s 
Pd 


(4-7) 


Rate  equations  of  the  general  form  of  Equations  (4-6)  or  (4-7)  can  also  describe  non¬ 
equilibrium  ionization  where  electron- ion  recombination  occurs  by  either  a  two-  or 
three-body  collision  process  without  an  energy  of  activation.  In  this  case  a  re¬ 
presents  the  degree  of  ionization  and  D  is  the  ionization  energy  per  mole. 


4.2.1  Numerical  Solutions 

Some  of  the  earliest  theoretical  studies  of  nozzle- flow  nonequilibrium  governed 
by  a  single  chemical  reaction  were  done  by  Penner  and  his  co-workers  around  1947  (see 
References  2  and  9  for  surveys).  One  of  the  first  numerical  solutions  was  obtained 
by  Krieger20  for  the  finite-rate  dissociation  and  recombination  of  hydrogen  in  a 
rocket  nozzle  flow.  Numerical  solutions  for  hypersonic  nozzle  flows  of  dissociating 
diatomic  gases  were  obtained  by  Bray17  for  the  Lighthill  gas  and  by  Hall  and  Russo18 
for  oxygen.  The  case  of  nonequilibrium  ionization  was  studied  by  Smith66,  Eschenroeder61, 
and  Rosner62.  More  recently,  comprehensive  analytical  treatments  of  nozzle-flow  non¬ 
equilibrium  for  the  Lighthill  gas  have  been  given  by  Blythe22  and  by  Cheng  and  Lee23,24. 

As  previously  mentioned,  the  numerical  solutions  for  monotonically  diverging  nozzles 
typically  exhibit  the  characteristic  feature  of  the  “freezing-out”  of  the  atom  mass 
fraction  a  ,  and  thus  the  energy  in  dissociation,  at  sufficiently  large  area  or  dis¬ 
tance  down  the  nozzle.  Illustrative  results  for  the  hypersonic  nozzle  flow  of  dis¬ 
sociated  oxygen,  taken  from  Reference  18,  are  shown  in  Figures  4-1  to  4-5.  These 
figures  show  the  computed  variation  of  a, ,  T  ,  p  ,  u  ,  and  p  with  nozzle  area  ratio 
A/A^  for  the  equilibrium  reservoir  condition  T0  =  5900°K  ,  pQ  =  82  atm  .  At  this 
reservoir  condition  the  oxygen  is  about  69%  dissociated  (aQ  =  0.6931).  The  nozzle 
geometry  for  these  calculations  was  axisymmetri'’  of  the  form  (A/A^)  =  1  +  (ax/L)2  , 
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with  a  throat  radius  L  =  0.5  nun  and  asymptotic  half-cone  angle  tan'1  a  =  7.5°  for 
large  ax/L  .  The  oxygen  molecule  was  assumed  to  be  a  harmonic  oscillator,  in  vibra¬ 
tional  equilibrium  throughout  the  expansion.  The  dissociation  rate  equation  for  a 
was  of  the  general  form  of  Equation  (4-6);  rate  constant  data  from  the  shock  tube 
studies  of  Byron26  were  used  in  the  calculations. 

Several  distinct  regions  may  be  identified  in  the  variation  of  a  as  typified  by 
Figure  4-1.  Following  departure  from  the  equilibrium  reservoir  state,  near  equilibrium 
conditions  are  maintained  initially,  in  this  case  until  some  distance  downstream  of 
the  throat.  A  rather  rapid  departure  from  equilibrium  then  occurs  over  a  limited 
transition  region.  Downstream  of  this  region,  a  approaches  asymptotically  a  constant 
frozen  value*.  The  temperature  (Fig.  4-2)  is  strongly  affected  by  freezing  in  this 
case  because  of  the  significant  fraction  of  the  total  enthalpy  frozen  in  dissociation. 
When  the  transition  to  frozen  dissociation  commences,  the  temperature  decreases  rapidly 
from  the  equilibrium  limit  as  the  heat  release  from  recombination  decreases.  The 
static  pressure  (Fig. 4-3)  is  not  as  strongly  affected,  but  the  reduction  in  pressure 
below  the  corresponding  equilibrium  limit  becomes  quite  large  at  large  area  ratios. 

The  velocity  and  density  (Figs. 4-4,  4-5)  are  least  affected  by  freezing,  and  are  de¬ 
creased  and  increased,  respectively,  relative  to  equilibrium  by  roughly  similar  amounts 
since  the  mass  flow  is  only  weakly  affected  by  nonequilibrium. 

The  frozen  degree  of  dissociation  is  decreased  by  an  increase  in  reservoir  density 
p0  ,  recombination  rate  constant  kR  ,  or  nozzle  scale  L/a  ,  and  by  a  decrease  in 
reservoir  temperature  TQ  .  The  numerical  studies  cited  show  that  the  effects  of 
p0  and  Tq  in  this  respect  are  much  stronger  than  of  kR  or  L/a  .  The  strong 
effect  of  density  level  is  due  to  the  dependence  of  XQ  on  pa  as  well  as  the  de¬ 
pendence  of  Dci/Dt  on  p2  .  Typically,  a  change  of  a  factor  of  two  or  three  in  kR 
or  L/a  has  little  effect  on  the  frozen  dissociation  level  (kR  and  L/a  enter  the 
now!  nsional  rate  equation  as  a  product).  Also,  the  dissociation  is  only  weakly 

afft  I  by  the  exact  temperature  dependence  of  kR  ,  i.e.  by  weather  s  has  the 

value  0,  1,  or  2,  for  example,  in  the  relation  kR  a  T~s  assumed. 

4.2.2  Approximate  Methods  of  Solution 

Approximate  methods  for  the  analysis  of  dissociating  nozzle  flows  have  been  deve¬ 
loped  by  various  authors,  for  example  by  Penner2,  Bray17  and  Hall  and  Russo18,  as  well 
as  by  others.  The  method  due  to  Bray  for  determining  the  approximate  freezing  point 
and  frozen  dissociation  level  in  terms  of  a  corresponding  equilibrium  solution  will 
be  briefly  discussed  since  it  has  been  widely  used  in  recent  years.  The  method  of 
Reference  18  is  essentially  similar,  and  Penner  has  discussed  the  similarity  of  his 
criterion  for  near-equilibrium  flow  to  Bray’s  method9. 

Bray’s  procedure  is  based  on  the  following  observations  concerning  the  rate  equation 
(Equation  (4-6))  expressed  as 

DCt 

-  -  Rn  —  Rd  , 

Dt  D  R 

where  RQ  and  Rr  are  the  rates  of  dissocia  ion  and  recombination,  respectively. 
Upstream  of  the  transition  region  where  freezing  commences,  the  flow  is  close  to 


Under  certain  conditions  to  be  discussed,  freezing  at  a  nonzero  value  of  may  not  occur. 
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equilibrium  and  and  Rr  are  large  and  very  nearly  equal.  In  this  near-equili¬ 

brium  region,  Da/Dt  is  small  compared  with  either  RD  or  Rr  .  Once  appreciable 
departure  from  equilibrium  has  occurred,  RD  rapidly  becomes  very  small  because  of 
the  exponential  Arrhenius  factor  e‘D/R°T  .  Downstream  of  the  transition  region, 
where  tne  flow  tends  toward  a  frozen  state,  RD  is  much  smaller  than  Rr  and 
Da/Dt  ia  Rr  »  Rd  .  Thus  in  the  change  from  near-equilibrium  to  near-frozen  flow, 

Rd  varies  from  »  Da/Dt  to  «  Da/Dt  .  At  some  point  in  the  transition  region, 
Da/Dt  must  therefore  be  of  the  same  order  as  RD  ,  Bray17  postulates  a  narrow 
transition  region,  or  sudden  freezing,  with  the  freezing  point  located  where  the 
relation 


/ba\ 

!—)  =  K0*D>e  (4-8) 


is  satisfied.  Here  K  is  a  constant  of  order  unity  (set  equal  to  1  by  Bray)  and  the 
subscript  e  indicates  evaluation  from  the  equilibrium-flow  solution.  The  approximate 
solution  obtained  thus  consists  of  an  equilibrium-flow  solution  upstream  of  the  freezing 
point  and  a  frozen- flow  solution  downstream. 

Since  this  approximate  solution  consists  of  two  isentropic  branches,  it  neglects 
the  increase  in  entropy  which  occurs  in  the  actual  nonequilibrium  flow.  However,  the 
entropy  rise  which  occurs  in  such  nonequilibrium  expansions  is  generally  very  small 
compared  with  the  reservoir  entropy  itself,  as  discussed  in  detail  by  Eschenroeder27. 

As  a  further  extension  of  the  sudden  freezing  approximation,  Bray28  showed  that  the 
gas  properties  at  the  freezing  point  (ap,  Tp,  hp  etc. )  depend  only  on  the  reservoir 
entropy  SQ  (and  not  also  on  hQ.  for  example)  for  a  particular  nozzle  geometry.  Thus 
within  the  accuracy  of  the  sudden  freezing  concept,  nonequilibrium  flows  in  a  given 
nozzle  may  be  represented  as  isentropes  on  a  Mollier-type  diagram  of  h  versus  SQ  . 
This  usefulness  of  the  reservoir  entropy  as  a  correlating  parameter  for  nonequilibrium 
nozzle  flows  extends  to  more  complex  chemistry,  as  discussed  in  Section  4.4  in  con¬ 
nection  with  nozzle  flows  of  air. 

Comparisons  made  by  various  authors  (e.g.  References  17,  18,  29)  of  the  approximate, 
sudden  freezing  solution  with  exact  numerical  solutions  show  the  sudden  freezing 
approximation  to  be  surprisingly  accurate  over  a  wide  range  of  conditions  despite  its 
empirical  nature  (see  Figure  4-1).  Typically,  the  sudden  freeze  approximation  predicts 
the  frozen  dissociation  level  to  within  a  few  percent.  Several  authors  have  improved 
on  this  simple  sudden-freeze  solution  by  taking  account  of  the  variation  of  a  in  the 
asymptotic  recombination  region  downstream  of  the  freezing  point.  This  was  done  by 
Smith21,  Eschenroeder61,  and  Rosner62  for  the  case  of  weak  ionization  decoupled  from 
the  gasdynamic  behavior  (where  the  rate  equation  is  of  the  general  form  of  Equation 
(4-7)  and  by  Eschenroeder 2 7  for  atom  recombination.  These  authors  obtained  improved 
asymptotic  solutions  by  solving  the  equation  Da/Dt  —  -  Rr  applying  in  the  near- 
frozen,  recombination  region  and  matching  to  the  equilibrium  solution  at  the  freezing 
point  determined  by  Equation  (4-8). 


4.2. 3  Analytical  Treatments 

As  distinct  from  numerical  or  semi-empirical  treatments,  rigorous  analytical  treat¬ 
ments  of  nozzle-flow  nonequ) librium  governed  by  the  nonlinear  Lighthill-gas  rate  f 

(Equation  (4-6))  have  been  given  by  Blythe22  and  by  Cheng  and  Lee23,24.  Blythe,  who 


51 


also  treats  separately  the  linear  rate  equation  for  vibrational  relaxation  (see 
Section  4.3).  considers  the  case  of  very  small  dissociation  such  that  the  appropriate 
(nonlinear)  rate  equation  is  Equation  (4-7)  which  is  assumed  to  be  completely  un¬ 
coupled  from  the  fluid  dynamic  equations.  Blythe’s  analysis  of  this  limiting  case 
reveals  three  basic  regions  to  the  flow  structure.  As  evidenced  in  the  numerical 
studies  of  previous  workers,  these  are:  an  upstream  near-equilibrium  region,  a  rather 
narrow  transition  region  where  rapid  large  departure  from  equilibrium  occurs,  and  a 
downstream  asymptotic  region  of  near- frozen  flow  where  recombination  dominates. 

Correctly  matched  analytic  solutions  to  these  three  regions  are  obtained.  Blythe 
shows  that  the  simple  sudden- freeze  solution  discussed  here  is  mathematically  incorrect 
in  the  limiting  case  of  weak  dissociation.  The  error  is  primarily  in  using  the  wrong 
asymptotic  solution  downstream  of  the  freezing  point  (i.e.  completely  frozen  flow) 
rather  than  in  matching  this  solution  directly  to  the  equilibrium  solution  (i.e. 
without  the  transition  region  rigorously  required).  In  this  respect,  the  asymptotic 
recombination  solutions  of  References  21,  27,  61,  and  62,  which  are  also  matched 
directly  to  the  equilibrium  solution,  are  a  better  approximation  than  the  sudden 
freeze  model. 

The  assumption  that  the  rate  equation  (Equation  (4-7))  can  be  uncoupled  from  the 
remaining  gasdynamic  equation  limits  a  in  practice  to  exceedingly  small  values. 

For  negligible  coupling,  the  energy  in  dissociation  (aD  per  mole,  say)  must  be  a 
completely  negligible  part  of  the  enthalpy.  Any  effect  on  temperature  due  to  dis¬ 
sociation  can  have  a  very  strong  effect  on  Dct/Dt  due  to  the  dissociation  term  in 
the  rate  equation  containing  the  exponential  factor  e~D/R°T  .  Cheng  and  Lee”  show 
that  the  completely  uncoupled  case  applies  only  for  *fp  «  e2  .  where  subscript  FP 
denotes  values  at  the  freezing  point  given  by  Equation  (4-8)  and  the  quantity 
e  =  (RQTFP/D)  is  assumed  «  1  .  Under  the  assumption  of  small  e  ,  these  authors 
define  two  basic  flow  regimes  in  terms  of  the  ratio  (XFP/e  ,  which  is  roughly  the 
ratio  of  the  energy  in  dissociation  to  that  in  the  translational  degrees  of  freedom. 

For  regime  I  defined  by  <xFp/e  <  0(e)  ,  the  gasdynamics  are  only  weakly  influenced  by 
the  chemistry,  whereas  for  the  regime  II  defined  by  CxFp/e)  >  0(e)  ,  the  gasdynamics 
are  strongly  influenced.  In  regime  I  the  analysis  of  Reference  23  reveals  the  three 
basic  flow  regions  found  by  Blythe.  However,  in  regime  II  the  analysis  reveals  four 
distinct  flow  regions:  upstream  equilibrium,  incipient  transition,  rapid  transition, 
and  downstream  recombination23. 

Cheng  and  Lee  mainly  treat  regime  II  and  obtain  analytical  solutions  for  the  four  flow 
regions  by  means  of  matched  asymptotic  expansions  in  the  small  parameter  e  .  They 
also  assess  the  validity  of  previous  approximate  methods.  In  contrast  to  regime  I, 
the  simple  sudden  freeze  solution  of  Equation  (4-8)  is  shown  to  be  valid  and  quite 
accurate  for  velocity  in  regime  II,  and  also  valid  for  the  frozen  dissociation  level 
if  (^fP/e)  »  1  ,  provided  the  dissociation  level  does  ultimately  freeze.  However, 
the  sudden  freeze  solution  for  temperature  is  not  rigorously  valid,  although  the 
numerical  error  may  typically  be  small.  The  authors23  also  determine  the  conditions 
for  which  the  asymptotic  downstream  flow  will  freeze  (i.e.  conditions  for  which 
lim  a  =  t  0)  for  nozzle  geometries  of  the  form  A  oc  x2m  .  For  example,  it  is 

X-a> 

shown  that  the  dissociation  level  will  freeze  (a^  i  0)  if  m  >  1/4  .  and  will  be 
’’self  limiting"  (\,  =  0)  if  m  <  1/4  ,  provided  |s|  <  1  .  Here  s  is  the  exponent 
for  the  inverse  temperature  dependence  of  kR  .  For  s  >  1  ,  the  dissociation  will 
\  freeze  for  m  >  3/{ 2(7  -s)}  and  will  be  self  limiting  for  m  <  1/4  .  Subsequent  to 

their  initial  study,  Cheng  and  Lee  have  further  generalized  their  analysis  to  include 
nozzle  geometries  of  arbitrarily  diverging  shape24. 
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4.2.4  Experimental  Studies 

Despite  this  extensive  amount  of  theoretical  study,  and  the  pressing  need  for  more 
precise  chemical  kinetic  data  in  nozzle-flow  environments,  nonequilibrium  nozzle  flows 
involving  a  single  dominant  reaction  have  received  rather  limited  experimental  study. 

The  lack  of  experiment  is  due  in  part  to  the  difficulty  of  creating  well  understood 
flows  with  significant  dissociation  of  diatomic  gases  as  oxygen  and  nitrogen,  as  well 
as  the  generally  difficult  diagnostic  problem  involved  in  determining  local  chemical 
composition. 

One  of  the  earliest  definitive  nozzle-flow  experiments  involving  a  single  reaction 
for  dissociation  and  three-body  recombination  was  done  by  Wegener30'32  for  the  reaction 

N2  +  N20u  «=?  N2  +  2N02  , 

which  involves  the  dissociation  of  nitrogen  tetroxide.  The  experiments  were  done  in 
a  carefully  designed,  closely  one-dimensional  supersonic  nozzle.  Small  concentrations 
of  N204  were  used  in  a  large  excess  of  inert  N2  carrier  gas.  The  rate  constants 
for  this  reaction  are  well  known  in  the  temperature  range  of  the  experiments  (about 
250  to  400°K) .  The  local  N02  concentration  along  the  nozzle  was  determined  by 
measuring  the  absorption  of  blue  light  (4300  X)  to  which  N2  and  N204  are  transparent. 
Appreciable  departures  from  equilibrium  composition  were  observed.  Wegener’ s  results 
generally  confirm  the  theoretical  expectations,  and  show  good  qualitative  agreement 
with  the  simple  sudden  freeze  model.  Further  experimental  results  obtained  for  the 
same  chemical  system,  but  in  a  strongly  two-dimensional  nozzle,  are  reported  by 
Zupnik  et  aZ.119. 

Additional  experimental  studies  of  single- reaction  chemistry  in  nozzle  flows  have 
included,  for  example,  studies  of  dissociated  nitrogen  flows  by  nozzle  static  pressure 
measurements  (Nagamatsu  and  Sheer33);  dissociated  oxygen  flows  by  ultraviolet  light 
absorption  measurement  of  O2  (Wilson3 “)  and  by  catalytic-probe  heat-transfer  measure¬ 
ment  of  0  (Reddy35;  Bartz132);  dissociated  hydrogen  flows  by  static  pressure  measure¬ 
ments  (Widawsky,  Oswalt,  and  Harp36;  Russo,  Hall,  and  Lordi37);  and  ionized  air  flows, 
where  the  dissociative- recombination  reaction  N0+  +  e';z£N  +  0  dominates,  by  micro- 
wave  interferometer  measurements  of  electron  density  (Eschenroeder  et  aZf3).  Generally 
speaking,  all  such  experiments  have  given  results  which  agree  qualitatively,  at  least, 
with  the  theoretical  predictions  based  on  available  kinetic  data  and  the  flow  models 
discussed.  What  few  recombination  rate  data  have  been  determined  directly  in  these 
experiments  do  not  contradict  the  usual  assumption  that  kD/kR  =  Kg(T)  under  non¬ 
equilibrium  conditions.  However,  in  most  cases  the  precision  of  the  recombination 
rate  constant  determination  is  too  low  to  provide  an  accurate  assessment  of  this 
relation.  In  particular,  static  pressure  measurement,  although  perhaps  the  easiest 
type  of  measurement  to  make,  is  inherently  a  difficult  approach  to  rate  constant 
determination  because  of  the  relative  insensitivity  of  pressure  to  nonequilibrium 
effects.  Improved  diagnostic  methods  are  needed  to  exploit  nozzle  flows  as  a  chemical 
kinetic  tool  with  the  precision  desired. 

4.3  Vibrational  Nonequilibrium  in  Nozzle  Flows  of  Diatomic  Gases 

In  the  preceding  discussion  of  nozzle- flow  nonequilibrium  involving  dissociation 
and  recombination,  the  question  of  simultaneous  vibrational  nonequilibrium  was  not 
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considered.  The  Lighthill  gas  model  effectively  assumes  that  vibration  is  always 
excited  to  one-half  the  classical  equipartition  value;  the  oxygen  nozzle-flow  solutions 
discussed  were  obtained  assuming  vibrational  equilibrium.  As  discussed  in  Section  3, 
the  coupled  effects  of  nonequilibrium  vibration  and  dissociation  are  a  matter  of 
current  research.  Simple  models  for  such  coupling  in  nozzle  flows  have  been  considered 
recently  by  several  authors  (e.g.  Bray  and  Pratt36,  and  Tirumalesa40).  However,  because 
of  the  lack  of  a  firmly  established  theoretical  model  of  coupling  useful  for  flowfield 
analysis,  most  calculations  of  chemically  reacting  nozzle  flows  have  assumed  either 
vibrational  equilibrium  or  else  a  vibrational  relaxation  process  uncoupled  from 
dissociation. 


The  present  section  will  be  limited  to  discussion  of  the  simplest  example  of  vibra¬ 
tional  nonequilibrium  in  nozzle  expansions,  namely  that  where  chemical  change  does 
not  occur  so  that  vibrational  relaxation  is  the  only  rate  process  involved.  An  im¬ 
portant  example  is  the  expansion  of  a  diatomic  gas  from  initial  nozzle  reservoir 
conditions  of  equilibrium  state  such  that  vibration  is  excited  but  dissociation  is 
negligible.  Such  reservoir  conditions  are  employed,  for  example,  in  lower-temperature 
operation  of  hot-shot  and  shock  tunnels,  and  in  steady-state  heated  nitrogen  tunnels. 

Of  course,  the  energy  in  vibrational  excitation  never  becomes  nearly  as  large  a  fraction 
of  the  total  enthalpy  as  does  the  energy  in  dissociation  at  higher  temperatures. 
Nevertheless,  the  vibrational  energy  is  not  negligible,  being  approximately  10%  of  the 
total  internal  energy  in  air  at  2000°K  for  example.  This  amount  of  internal  energy  is 
enough  to  produce  nontrivial  gasdynamic  effects  if  appreciable  departure  from  equilibrium 
nozzle  flow  occurs.  Thus  vibrational  nonequilibrium  can  produce  effects  of  some  im¬ 
portance  on  nozzle  test-flow  conditions  in  these  types  of  wind-tunnel  facilities. 


lt.3.1  Analysis  for  Harmonic  Oscillator  Model 

As  discussed  previously  in  Section  3,  the  simplest  equation  which  approximately 
describes  the  exchange  of  energy  between  vibrational  and  translational  degrees  of 
freedom  in  a  pure  undissociated  diatomic  gas  is  the  linear  equation  due  to  Landau  and 
Teller: 


Dt 


(4-9) 


Here  ev  is  the  vibrational  energy  per  unit  mass  of  gas  and  ey  Eq  is  the  value  of 
ey  when  vibration  is  equilibrated  with  translation.  The  quantity  Ty  is  the  relaxa¬ 
tion  time  for  vibrational-translational  energy  exchange;  ry  is  inversely  propor¬ 
tional  to  collision  frequency  (or  density)  and  increases  strongly  with  decreasing 
translational  temperature.  As  pointed  out  in  Section  one  of  the  several  key  assump¬ 
tions  on  which  this  equation  is  based  is  that  the  molecules  behave  as  simple  harmonic 
oscillators  for  which  the  vibrational  energy  levels  are  equally  spaced.  In  this  case, 
the  exchange  of  vibrational  energy  between  the  vibrational  modes  of  two  colliding 
molecules  is  a  resonance  phenomenon  which  serves  only  to  redistribute  vibrational 
energy  without  exchange  with  translational  degrees  of  freedom.  If  a  Boltzmann  popula¬ 
tion  of  the  vibrational  energy  levels  is  maintained,  as  is  typically  the  case  for 
harmonic  oscillator  molecules,  the  vibrational  mode  can  be  characterized  by  a  vibrational 
temperature  Ty  .  For  the  harmonic  oscillator  model,  ey  is  then  given  in  terms  of 
Ty  by: 

e  -  _ ^ _ 

v  _  k/i  (ehl//kTv  -  1)  ' 
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Despite  its  simplifying  assumptions,  Equation  (4-9)  has  been  widely  used  to  analyze 
vibrational  relaxation  in  nozzle  flows.  Although  the  equation  is  linear,  numerical 
solution  is  still  generally  necessary  for  the  calculation  of  nozzle  flows,  where  both 
density  and  translational  temperature  vary  strongly.  However,  because  the  vibrational 
energy  is  typically  a  rather  small  fraction  •  f  the  total  enthalpy,  Equation  (4-9)  can 
often  be  decoupled  from  the  associated  fluid  dynamic  equations  without  introducing 
significant  error. 

Numerical  solutions  based  on  Equation  (4-9)  have  been  obtained  by  several  authors 
for  diatomic-gas  nozzle  flows  (e.g.  References  41-44).  Most  of  these  calculations 
have  used  values  for  the  relaxation  time  ry  deduced  from  shock-tube  experiments  on 
vibrational  relaxation  behind  plane  shock  waves.  The  calculations  of  Stollery  and 
Smith42  for  vibrational  relaxation  in  conical-nozzle  flows  of  undissociated  air  give 
typical  results.  As  the  gas  expands  down  the  nozzle  from  the  initial  equilibrium 
state,  the  vibrational  temperature  Ty  decreases  more  slowly  than  the  rapidly  de¬ 
creasing  translational  temperature.  At  some  point  in  the  diverging  nozzle,  depending 
on  the  initial  reservoir  conditions  and  the  nozzle  geometry,  the  vibrational  energy 
(or  Ty)  freezes  or  levels  off  at  a  constant  value  higher  than  the  local  translational 
temperature. 

This  freezing  of  vibrational  energy  is  analogous  to  the  freezing  of  molecular 
dissociation  discussed  previously.  The  gasdynamic  effects  are  qualitatively  the  same 
as  for  dissociation  freezing,  but  typically  much  smaller.  Of  course,  in  the  case  of 
vibrational  relaxation  only  binary  collisions  are  involved,  whereas  ternary  collisions 
are  also  important  in  chemical  nonequilibrium  involving  three-body  recombination.  It 
may  be  noted  that  as  vibrational  relaxation  involves  only  two-body  collisions,  binary 
scaling  can  be  applied  to  solutions  of  Equation  (4-9)  for  undissociated  flows  in 
nozzles  of  given  shape.  For  the  nozzle  geometry  mentioned  in  Section  4.2,  for  example, 
solutions  can  thus  be  correlated  or  scaled  in  terms  of  the  parameter  pQ  (L/a)  ,  where 
p0  is  the  initial  reservoir  density,  for  a  given  gas  and  initial  reservoir  temperature 
or  enthalpy. 

Rigorous  analytic  treatment  of  Equation  (4-9)  for  nozzle  flows  has  been  given  by 
Blythe22,47  for  the  limiting  case  where  the  vibrational  energy  is  sufficiently  small 
that  Equation  (4-9)  can  be  decoupled  from  the  fluid  dynamic  equations.  Blythe’s 
analysis  and  solutions  reveal  a  gross  mathematical  similarity  to  the  dissociation 
freezing  pattern,  in  that  three  basic  regions  (as  discussed  in  Section  4.2)  character¬ 
ize  the  flow  structure.  However,  the  detailed  asymptotic  behavior  for  vibrational 
freezing  differs  from  dissociation  due  to  the  linearity  of  Equation  (4-9)  as  opposed 
to  the  nonlinearity  of  Equation  (4-6). 

As  for  the  case  of  dissociation  freezing  in  nozzle  flows,  semi-empirical  approximate 
methods,  including  similar  sudden  freeze  approximations,  have  been  proposed  and  used 
to  estimate  vibrational  fieezing  from  Equation  (4-9)  (see  References  41,  43,  45,  46). 
Although  these  methods  typically  give  results  which  agree  quite  well  with  numerical 
solutions,  Blythe22  shows  them  to  be  mathematically  incorrect  as  in  the  case  of  weak 
dissociation  freezing  discussed  in  Section  4.2. 

4.3.2  Experimental  Studies 

In  contrast  to  the  relative  scarcity  of  experimental  data  on  dissociation  non¬ 
equilibrium  in  nozzle  flows,  considerable  effort  has  been  expended  in  recent  years  on 
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nozzle-flow  experiments  involving  vibrational  relaxation.  One  reason  for  this  activity 
is  the  results  of  initial  experiments  by  Hurle,  Russo,  anu  Hall48  in  1962-63  to  measure 
vibrational  temperature  relaxation  in  nozzle  flows  of  undissociated  nitrogen.  The 
vibrational  temperatures  determined  in  these  experiments  showed  only  qualitative  agree¬ 
ment  at  best  with  the  predictions  of  equation  (4-9)  based  on  values  of  ry  (=  ryg) 
from  shock  wave  experiments.  Quantitatively,  the  measured  values  of  Ty  were  con¬ 
siderably  lower  than  predicted,  which  suggested  faster  vibr»tional  de-excitation  than 
given  by  ry  =  ryg  These  initial  experiments  stimulated  further  experimentation 
(Refs. 51-54, 56-58) ,  and,  more  recently,  theoretical  studies  concerning  effects  not 
included  in  Equation  (4-9),  particularly  the  effects  of  anharmonicity5  9,  60. 

The  experiments  of  Reference  48  utilized  the  particular  spectrum-line  reversal 
method  due  to  Gaydon  and  Hurle5  49  as  the  diagnostic  technique  to  determine  the  vibra¬ 
tional  temperature  of  N2  at  different  stations  along  a  15°-conical  nozzle  driven  by 
a  shock  tube.  The  technique  involves  simultaneous  monitoring  of  the  emission  and 
absorption  intensities  of  a  suitable  spectrum  line  from  a  trace  species  (Na  and  Cr 
were  used)  added  to  the  test  gas.  Interpretation  of  a  vibrational  temperature  depends 
on  close  equilibration  of  the  electronic  transition  with  the  vibrational  mode.  Typical 
results  from  these  experiments  in  the  form  of  Ty/T0  versus  A/A#  are  shown  in 
Figure  4-6  for  both  Na  and  Cr  as  the  trace  additives  for  line-reversal  measurement. 
The  uppermost  curve  in  Figure  4-6  is  the  prediction  of  Equation  (4-9)  with  ry  equal 
to  Tyg  ,  the  value  deduced  from  shock-wave  experiments.  The  lowest  curve  is  for 
complete  equilibrium  flow,  i.e.  ry  =  0  .  The  experimental  points  are  fitted  by  the 
curve  shown,  which  is  calculated  from  Equation  (4-9)  with  ry  equal  to  (1/70)  ryg  . 
Thus  the  data  suggest  that  the  vibrational  de-excitation  of  N2  in  the  nozzle  flow 
proceeds  much  faster  than  predicted  by  the  Landau-Teller  relaxation  equation  on  the 
basis  of  the  shock-wave  relaxation  times*. 

In  view  of  the  unexpected  nature  of  these  results,  further  experimentation,  in¬ 
cluding  the  use  of  different  diagnostic  techniques,  has  been  done  on  vibrational 
relaxation  in  nozzle  flows  of  N.  .  In  total,  these  further  studies  support  the 
results  of  Reference  48.  Further  Na  line-reversal  experiments  include  studies  of 
the  influence  of  impurities  (Russo51),  and  the  behavior  of  dilute  N2-Argon  mixtures 
with  attention  to  the  effects  of  free  electrons  (Hurle  and  Russo59).  Different 
diagnostic  techniques  applied  include  the  electron  beam  method  (Petrie53  and 
Sebacher54' 55)  where  the  N?  is  ionized  and  electronically  excited  by  a  narrow  beam 
of  electrons  passing  through  the  flow  field.  The  vibrational  temperature  of  ground 
state  N2  is  deduced  from  the  ratios  of  intensities  of  different  vibrational  bands 
in  the  observed  emission  spectrum  (usually  of  N*).  A  different  reversal  spectroscopic 
technique  has  been  used  recently  by  Russo56  wherein  CO  is  added  in  trace  amounts  to 
the  N;  as  a  thermometric  element.  Band  reversal  spectroscopy  using  the  fundamental, 
infrared  vibration-rotation  band  of  CO  permits  determination  of  the  CO  vibrational 
temperature,  from  which  the  coupled  vibrational  temperature  of  N;  can  be  deduced*. 

*  The  factor  1/70  is  a  revision  froa  the  earlier  value  1/15  given  in  Reference  48.  The  value 
of  this  factor  depends  on  the  reference  values  chosen  for  ryg  .  The  revision  to  1/70  from 
1/15  is  based  on  shat  is  believed  to  be  a  better  correlation  of  data  for  ry8  ,  as  discussed 
In  Reference  50. 

1  In  addition  to  these  experiments  on  N2  vibration  in  steady  expanding  floss  generated  in  a 
nozzle,  Holbecbe  and  loodley58  have  used  the  Na  line  reversal  method  to  determine  vibrational 
temperatures  in  one -dimensional ,  unsteady  expansion  flows  generated  in  a  shock  tube.  Their 
results  show  a  similarly  fast  relaxation  for  N;  ,  but  not.  however,  for  CO  . 
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As  already  noted,  the  further  studies  of  N2  vibrational  relaxation  in  nozzle 
flows  generally  support  the  finding  of  Reference  48  that  the  effective  relaxation 
time  in  the  nozzle  flow  is  smaller  than  behind  the  shock  wave  at  tie  same  translational 
temperature  and  density.  The  nozzle-flow  vibrational  temperatures  measured  by  the 
different  techniques  and  experimenters  appear  to  be  in  good  quantitative  agreement,  as 
discussed  in  References  50  and  55  .  Based  on  the  experiments  of  References  48,  53, 
and  54  and  Equation  (4-9),  Sobacher55  has  deduced  an  empirical  expression  for  the 
variation  with  temperature  of  tne  vibrational  relaxation  time  of  N2  in  nozzle  flows, 
which  is  then  usable  with  Equation  (4-9)  for  the  purpose  of  calculations.  This  expres¬ 
sion  is 

ry p  =  3  x  10" 12  exp  (181  T" 1  /  3 )  atm  sec  , 

where  ry  is  the  relaxation  time  in  seconds,  and  p  the  local  static  pressure  in 
atmospheres.  Although  Equation  (4-9)  with  this  expression  for  ry  will  give  pre¬ 
dictions  for  Ty  in  agreement  with  the  above-mentioned  experiments,  it  must  be 
emphasized  that  such  a  procedure  is  empirical,  and  without  a  rational  basis  at  present. 
As  mentioned  at  the  end  of  this  section,  recent  theoretical  studies  suggest  that 
enharmonic  effects  may  be  partly  responsible  for  the  different  behavior  of  N2  vibra¬ 
tional  relaxation  in  nozzle  flow  and  shock  wave  environments.  Such  effects  are  not 
included  in  Equation  (4-9). 

Similar  nozzle-flow  experiments  with  CO  ,  again  using  the  Na  line- reversal  method 
to  measure  vibrational  temperatures,  were  done  by  Russo  and  Rich  at  the  Cornell  Aero¬ 
nautical  Laboratory50.  Qualitatively  similar  results  as  for  N„  were  observed  in 
that  the  measured  vibrational  temperatures  were  closer  to  equilibrium  values  than 
predicted  by  Equation  (4-9)  using  shock-wave  relaxation  times.  More  recently,  Russo57 
has  repeated  this  experiment  with  CO  but  using  the  infrared  band- reversal  method 
mentioned  above  to  measure  Ty  .  Similar  results  were  obtained  as  with  the  Na  line- 
reversal  method. 

In  Reference  48  it  was  speculated  that  the  apparent  faster  vibrational  relaxation 
in  the  N2  nozzle- flow  expansions  as  compared  to  shock  wave  flows  was  possibly  due 
in  part  to  differing  effects  of  anharmonicity  under  the  different  thermo-gasdynamic 
conditions  and  initial  states  of  molecular  excitation  existing  in  the  two  environments. 
In  the  shock  wave  flow,  where  the  net  relaxation  process  is  excitation,  the  vibrational 
temperature  is  initially  much  less  than  the  translational  temperature.  Also,  the 
translational  temperature  and  gas  density  remain  roughly  constant  through  the  relaxation 
process.  In  the  nozzle  expansion  flow,  however,  the  net  relaxation  process  is  de¬ 
excitation.  In  this  case  the  vibrational  temperature  is  initially  high,  and  rapidly 
becomes  substantially  higher  than  the  translational  temperature  as  the  latter  decreases 
strongly  and  a  vibrational  lag  develops  during  flow  expansion. 

Recent  theoretical  studies  suggest  that  anharmonic  effects  may  indeed  be  a  factor 
in  the  above  respect5  9,60.  As  discussed  in  Section  3,  with  anharmonic  oscillators 
the  vibration-vibration  exchange  process  can  produce  very  large  departures  from  a 
Boltzmann- like  population  of  the  vibrational  energy  levels  when  the  population  of  the 
low  lying  levels  exceeds  that  for  equilibrium  with  the  translational  temperature  The 
effect  is  to  overpopulate  the  higher  vibrational  levels,  and  for  sufficiently  extreme 

•  The  recent  work  of  Russo56  using  infrared  band  reversal  with  CO  as  a  trace  additive  provides 
additional  data  beyond  that  discussed  in  References  50  and  55. 
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conditions,  a  large  population  inversion  can  occur.  As  vibration-translation  energy 
exchange  becomes  more  efficient  at  higher  vibrational  levels,  this  effect  can  increase 
the  net  rate  of  vibrational  de-excitation.  The  low  vibrational  levels  can  become 
significantly  underpopulated  relative  to  predictions  of  the  linear  Landau-Teller 
relaxation  equation.  The  conditions  necessary  for  anharmonic  vibration-vibration 
exchange  to  affect  significantly  the  vibrational  population,  i.e.  high  vibrational 
energy  and  low  translational  energy,  tend  to  occur  in  the  nozzle  expansion-flow  environ¬ 
ment,  but  not  behind  a  shock  wave.  Rich  and  Rehm59,  and  Bray60  have  recently  considered 
theoretically  the  relaxation  of  anharmonic  oscillator  models  which  exhibit  such  effects. 
Bray's60  work  includes  calculations  of  anharmonic  relaxation  effects  for  specific 
comparison  with  the  N2  nozzle- flow  experiments  of  Reference  48.  Although  there  is 
necessarily  much  uncertainty  in  the  basic  transition  probabilities  assumed  in  the 
theoretical  model,  it  is  encouraging  that  the  calculated  results  show  qualitative 
agreement  with  experiment. 

4.4  Chemical  Nonequilibrium  in  Nozzle  Flows  of 
Coup led -React ion  Gas  Mixtures 

Sections  4.2  and  4.3  have  discussed  the  simplest  examples  of  nozzle-flow  non¬ 
equilibrium,  involving  only  single  rate  processes  governing  dissociation  or  vibrational 
relaxation  in  diatomic  gases.  Of  course,  in  technical  applications  involving  non¬ 
equilibrium  nozzle  flows  the  working  fluid  is  often  a  multi-species  gas  mixture  in¬ 
volving  coupled  chemical  reactions,  rather  than  a  pure  diatomic  gas.  This  is  the 
case  with  high-enthalpy  airflows  in  hypervelocity  test  facilities,  and  is  typically 
the  case  also  with  nozzle  flows  of  interest  in  the  propulsion  field.  This  section 
will  review  briefly  some  typical  studies  of  such  coupled-reaction  nozzle  flows.  The 
example  of  high-enthalpy  nozzle  airflow  is  emphasized  because  of  its  relative  im¬ 
portance;  some  considerations  of  nonequilibrium  effects  in  hypervelocity  wind-tunnel 
testing  are  included  (Section  4.4.2). 

As  discussed  in  Section  4.1,  machine  computational  methods  are  now  fairly  well 
developed  for  obtaining  "exact”  numerical  solutions  to  nonequilibrium  nozzle  flows  of 
complex  chemically  reacting  mixtures.  Typically,  present  machine  codes  for  this 
purpose  are  based  on  the  simplifying  gasdynamic  assumptions  of  pseudo-one-dimensional 
and  inviscid  flow.  In  particular,  the  one-dimensional  flow  assumption  does  not  re¬ 
present  a  fundamental  limitation,  but  is  more  a  matter  of  practical  convenience  as 
regards  the  costs  and  limitations  of  present  machine  computing  techniques.  It  may  be 
noted  that  what  studies  have  been  made  of  two-  or  three-dimensional  inviscid  nozzle 
flows  with  nonequilibrium  chemistry  suggest  that  departuies  from  a  one-dimensional 
inviscid-core  behavior  may  often  be  significant,  particularly  as  regards  the  chemical 
composition  (see,  for  example,  References  119-122).  Thus,  inclusion  of  two-dimensional 
effects  in  the  supersonic  inviscid  core  will  likely  become  more  common  (for  example, 
using  the  method  of  characteristics)  as  the  general  capabilities  of  computing  mechines, 
and  also  the  knowledge  of  reaction  rate-constant  data,  continue  to  improve. 

Of  course,  particularly  with  complex  mixtures  the  various  simplifications  and  un¬ 
certainties  ;n  the  phenomenological  description  of  the  physical-chemical  rate  processes 
involved  constitute  the  real  limitation  of  present  nozzle-flow  machine  codes.  With 
complex  kinetics,  a  description  of  the  coupling  of  chemical  reaction  rates  with  internal 
degrees  of  excitation  is  generally  not  attempted.  All  internal  modes,  including 
vibration,  are  usually  assumed  to  be  in  thermal  equilibrium  with  the  translational 
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temperature.  Thus  the  forward  and  reverse  specific  rate  constants  of  all  reactions 
are  related  through  the  equilibrium  constants  in  the  usual  way.  With  ionized  nozzle 
flows  the  problem  is  further  complicated  by  the  effects  of  Coulomb  forces  and  the 
possible  lack  of  equilibrium  between  the  light  free  electrons  and  the  translational 
energy  of  the  heavy  particles.  No  doubt  the  mathematical  description  of  the  kinetics 
of  complex  mixtures  will  continue  to  improve  as  knowledge  increases  of  important 
elementary  reactions  and  rate  constants,  as  well  as  of  the  influence  of  internal 
degrees  of  excitation. 

The  approximate  methods  of  analysis  described  in  Section  4.2  for  dissociating 
diatomic  gases  have  had  some  extension  to  multi-component  coupled-reaction  nozzle 
flows,  principally  for  gas  mixtures  of  interest  in  propulsion  applications  (e.g. 
References  8,  78,  81,  83,  84,  88.  89).  However,  the  general  applicability  of  the 
simple  sudden-freeze  approximation  to  complex  kinetics  systems  appears  rather  limited. 
Apart  from  the  question  of  mathematical  correctness,  which  exists  even  for  the  pure 
diatomic  gas  case  (Section  4.2),  the  detailed  manner  of  application  of  sudden-freeze 
criteria  becomes  quite  arbitrary  and  empirical  with  complex  systems.  The  basic 
difficulty  is  that  in  a  complex  kinetics  system,  neither  all  species  nor  reactions 
“freeze”  simultaneously  in  general.  Typically,  complex  mixtures  of  interest  involve 
fast  bimolecular  exchange-type  reactions  in  addition  to  coupled  dissociation-recombina¬ 
tion  reactions.  Such  bimolecular  reactions  do  not  tend  to  freeze,  and  can  provide 
such  a  fast  path  for  atom  removal  that  atom  freezing  does  not  necessarily  occur. 

Wilde88  gives  a  brief  but  informative  summary  of  various  approaches  taken  to  apply 
the  simple  sudden-freeze  concept  to  complex  mixtures.  These  approaches  include  applying 
the  sudden-freeze  criterion  of  reactions8- 814  to  species78,  and  to  the  average  molecular 
weight  of  the  mixture83.  Perhaps  the  main  usefulness  of  approximate  methods  for  complex 
kinetics  at  present  is  where  the  basic  chemical  kinetic  behavior  of  the  flow  is  already 
fairly  well  known  from  more  exact  methods  of  solution.  Certainly  where  applicable  and 
carefully  used,  approximate  methods  can  provide  worthwhile  reductions  in  computing 
effort.  This  is  illustrated,  for  example,  by  Migdal  and  Goldford122  who  greatly 
simplify  the  analysis  of  nozzle-flow  two-dimensional  effects  by  applying  a  sudden- 
freeze  criterion  along  individual  streamlines.  The  streamline  gasdynamic  histories 
and  geometries  are  initially  determined  by  the  method  of  characteristics  for  equilibrium 
or  frozen  flow. 

4.4.1  !\onequilibrium  \'ozzle  Flows  of  Propulsion  Mixtures 

The  past  five  or  siy  years  has  seen  intensive  application  of  the  machine  codes  to 
study  nonequilibrium  nozzle  flows  of  complex  mixtures  specifically  of  interest  in 
various  propulsion  applications,  as  well  as  of  air.  Of  course,  chemical  freezing  in 
rocket  exhaust  nozzles  has  long  been  of  interest  because  of  the  consequent  loss  in  rocket 
thrust  or  specific  impulse.  Usually,  however,  such  effects  are  not  of  major  importance 
with  solid  or  liquid  propellant  rockets  of  modest  performance.  The  limits  of  equili¬ 
brium  and  frozen  flow  might  typically  give  thrust  differences  of  several  percent,  so 
that  with  the  characteristically  high  pressures  of  rocket  chambers  and  modest  nozzle 
expansion  rates,  the  percentage  impulse  loss  due  to  nonequilibriupi  effects  tends  to 
be  quite  small.  Nonequilibrium  nozzle-flow  effects  can  assume  much  greater  importance 
with  high  performance  rocket  systems  and  hypersonic  ramjet  applications10,11,80.  In 
the  case  of  the  hypersonic  ramjet,  the  net  propulsive  thrust  represents  essen'  ’ rlly 
the  relatively  small  difference  between  the  exhaust  and  inlet  momentum  fluxes,  both 
of  which  are  large  quantities.  Thus  nonequilibrium  effects  of  a  magnitude  considered 
unimportant  for  rocket  performance  can  be  quite  important  for  hypersonic  ramjet 
performance. 


59 


A  comprehensive  review  of  research  in  this  area  up  to  1962  is  given  by  Olson11. 
Representative  numerical  studies  of  recent  years  are  those  of  Westenberg  and  Favintu, 
Sarli  et  alas,  and  Wilde87,88.  Westenberg  and  Favin  did  a  comprehensive  study  of 
nonequilibrium  nozzle  flows  for  two  important  types  of  mixtures:  one  typical  of  rocket 
exhaust  products,  containing  H2  ,  H20  ,  CO  ,  C02  ,  H  ,  and  OH  and  described  kine- 
tically  by  four  elementary  reversible  reactions;  and  a  second  typical  of  ramjet  exhaust 
products  after  hydrogen-air  combustion,  containing  H?  ,  H20  ,  02  ,  H  ,  OH  ,  and  0 
and  described  kinetically  by  eight  elementary  reversible  reactions.  Particularly,  for 
the  ramjet  example,  the  results  demonstrate  the  importance  of  bimolecular  reactions 
such  as 


(a) 

H20  +  H  H2  +  OH  , 

(c) 

h2  +  0 

OH 

+  H 

(b) 

o2  +  H  5=±0H  +  0  , 

(d) 

h2o  +  0 

i=t  OH 

+  H 

which  are  typically  fast  compared  to  the  several  dissociation-recombination  reactions 
involved. 

For  the  example  of  H2-Air  combustion,  Figure  4-7  shows  results  from  Reference  84 
for  calculated  equilibrium  and  nonequilibrium  distributions  of  H?  ,  H  ,  and  OH 
along  a  25°  conically  divergent  nozzle.  The  mixture  expands  from  supersonic,  initial 
equilibrium  conditions  of  3000°K  and  5  x  105  cm/sec  velocity  at  the  initial  pressures 
and  H2/air  stoichiometric  equivalence  ratios  (E. R. )  shown.  The  departures  from  equili¬ 
brium  composition  are  seen  to  be  substantial.  Increased  pressure  level  produces  a 
flow  closer  to  equilibrium,  as  in  the  case  of  the  pure  diatomic  gas.  Also,  the  gas- 
dynamic  consequences  of  nonequilibrium  are  similar  to  the  pure  diatomic  gas  case  - 
local  pressure,  translational  temperature,  and  velocity  are  reduced  below  corresponding 
equilibrium  values.  The  faster  than  equilibrium  rate  of  decrease  of  H?  and  OH  at 
low  pressure,  and  the  downstream  increase  in  H  at  high  pressure,  are  interesting 
features  which  are  dile  to  the  dominance  of  the  above  bimolecular  reactions.  As  will 
be  discussed,  bimolecular  reactions  play  a  similarly  important  role  in  the  case  of  air, 
where  analogous  reactions  to  (b)  and  (c)  must  be  considered  involving  NO  . 

In  the  study  by  Sarli  el  al ?5,  machine  calculations  of  a  similar  nature  were  done 
for  H2-air  combustion  products  expanding  in  a  subsonic- entry  ramjet  exhaust  nozzle. 

In  this  case  the  chemical  kinetic  model  also  included  the  bimolecular  and  dissociation- 
recombination  reactions  involving  N  and  NO  which  had  previously  been  demonstrated 
to  be  important  in  high-temperature  pure  air  expansions98, 97 .  However,  in  the  example 
considered,  the  inlet  temperature  before  expansion  (about  2800°K)  was  really  not  high 
enough  for  N  and  NO  to  play  a  significant  role.  The  nozzle  entrance  conditions 
in  this  example  were  near-equilibrium  conditions  resulting  from  a  separate  machine 
calculation  of  the  H?-air  combustion  process  in  a  subsonic  constant-area  combustor. 

The  same  kinetic  system  and  calculation  methods  were  applied  by  the  authors  to  estimate 
the  relative  importance  of  thrust  losses  due  to  skin  friction  and  chemical  nonequilibrium 
for  a  simplified  supersonic-combustion  ramjet  configuration  employing  a  conical  exhaust 
nozzle86.  The  results  indicate  that  skin  friction  losses  are  dominant,  but  that  for 
optimum  nozzle  expansion  angles  and  high  thrust  levels  the  nonequilibrium  losses  arc 
of  comparable  magnitude  at  the  flight  conditions  considered  (flight  Mach  number  of 
15  at  150,000  ft  altitude). 
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Wilde  has  performed  careful  numerical  studies  to  assess  nozzle-flow  nonequilibrium 
effects  for  several  complex  kinetic  mixtures  which  give  high  rocket  impulse.  The 
study  of  Reierence  87  concerns  H2  -  F2  mixtures  while  Reference  88  considers 
UDMH/N2H4  -  N204  and  also  CH2  -  02  systems.  Particularly  in  the  case  of  hydrogen- 
fluorine.  the  combustion  temperature  is  sufficiently  high  that  appreciable  energy  can 
be  stored  in  molecule  dissociation  prior  to  flow  expansion  in  th.  nozzle.  Thus  the 
thrust  loss  from  chemical  freezing  becomes  of  greater  importance  with  such  high  per¬ 
formance  systems11.  In  both  studies,  impulse  values  are  calculated  using  the  assumed 
kinetic  models  and  rate  constants  for  comparison  with  selected  experimental  data. 

For  CH2  -  02  ,  the  calculated  exhaust  temperatures  agree  well  with  experimental  values 
reported  by  Boynton123  for  the  RP-l/02  (gas)  propellant  system  over  a  range  of  oxidizer/ 
fuel  ratios  giving  appreciable  nonequilibrium  effects.  The  calculated  vacuum  impulse 
exceeds  experimental  values  of  Fortini,  Hendrix,  and  Huff124  obtained  for  RP-l/02 
(liquid)  by  amounts  up  to  five  or  six  percent.  This  discrepancy  in  impulse  is  attri¬ 
buted  to  unaccounted  experimental  losses,  rather  than  inadequacy  of  the  chemical 
kinetic  model  used  in  the  calculations.  In  the  case  of  the  H2  -  F2  system,  the 
calculated  vacuum  impulse  agrees  very  closely  with  experimental  data  (leported  by 
Olson11)  at  higher  chamber  pressures  for  near  stoichiometric  mixtiues,  but  falls  a 
few  percent  below  experimental  values  with  increased  amounts  of  H2  .  At  near  stoichio¬ 
metric  ratios,  the  impulse  loss  due  to  nonequilibrium  effects  is  about  5  to  6%  for  a 
chamber  pressure  of  20  atm,  a  conical  nozzle  half-angle  of  15°,  an  expansion  area 
ratio  of  25,  and  a  throat  diameter  of  10  cm. 

Wilde  interprets  the  relatively  high  experimental  values  of  impulse  with  hydrogen 
rich  H2  -  F2  mixtures  in  terms  of  HF  being  unusually  effective  as  a  catalytic 
third  body  in  H-atom  recombination.  On  this  basis,  the  three-body  recombination  rate 
constant  with  HF  as  the  third  body  is  interpreted  as  the  order  of  1000  times  larger 
than  with  H2  as  the  third  body.  The  author  also  speculates  that  the  inferred  rate 
enhancement  may  be  due  to  formation  of  the  stable  intermediate  species  H2F  and  a 
consequent  rapid  two-body  path  for  H-atom  removal,  i.e. 

H  +  HF  5ZS  H2F 

H  +  H2F  H2  +  HF  . 

The  interesting  possibility  of  deliberately  introducing  such  two-body  chain  reactions 
(as  a  faster  path  for  atom  recombination  than  three-body  reactions)  by  addition  of 
small  amounts  of  a  catalyst  body  to  H.  has  been  considered  by  several  investigators 
(Refs. 92,  93,  125).  Nonequilibrium  nozzle-flow  calculations  to  assess  such  possible 
gas-phase  catalysis  of  H-atora  recombination  are  described  in  Reference  92  for  seeding 
with  carbon,  in  Reference  93  for  seeding  with  oxygen  and  oxygen-nitrogen  mixtures,  and 
in  Reference  87  for  seeding  with  fluorine.  In  general,  the  results  of  these  numerical 
studies  suggest  that  such  catalytic  seeding  of  dissociated  hydrogen  offers  only 
marginal  gains  in  impulse  because  of  the  molecular  weight  penalty  that  must  also  be 
incurred.  More  definitive  theoretical  assessment  of  this  technique  for  reducing  nozzle- 
flow  dissociation  levels  requires  more  exact  rate  data  for  the  elementary  reactions 
involved. 

Before  turning  to  studies  of  high-enthalpy  airflows,  further  mention  may  be  made 
of  nozzle-flow  experiments  with  complex  mixtures  of  interest  in  propulsion  systems. 

In  general,  the  emphasis  in  such  experiments  is  on  measurement  of  the  gross  thermo- 
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gasdynamic.  characteristics  as  pressure,  temperature,  and  nozzle  thrust,  for  example. 
Representative  studies  of  recent  years  include  the  experiments  of  Hoglund,  Carlson, 
and  Byron90,  on  recombination  lag  in  nozzle  flows  of  the  RP- 1/0  2  (gas)  system,  the 
experiments  of  Lezberg  and  Franciscus82  on  nozzle  flows  of  hydrogen-air  and  also 
methane-air  combustion  products,  and  the  experiments  of  Boynton123  with  RP-l/02  already 
mentioned.  Typically,  the  observed  departures  of  the  gasdynamic  quantities  from  their 
equilibrium  behavior  are  predicted  fairly  well  by  exact  finite-rate  calculations,  in 
spite  of  the  considerable  complexity  of  the  kinetic  models  and  uncertainties  in  many 
of  the  rate  constants  involved.  While  this  is  encouraging,  improved  rate  constant 
data  are  certainly  needed,  as  are  additional  nozzle- flow  experiments  wherein  species 
concentrations  are  also  determined.  On  the  difficult  matter  of  determining  species 
concentrations,  the  recent  experiments  of  Hill,  Unger,  and  Dickens91  may  be  noted. 

These  authors  used  infrared  emission  and  absorption  spectroscopy  to  determine  the 
concentration  of  CO  ,  as  well  as  temperature,  in  nozzle  expansions  of  solid-propellant 
combustion  products  containing  CO  ,  C02  ,  H20  .  H2  ,  II  ,  and  OH  .  The  concentrations 
of  C02  along  the  nozzle  axis  determined  from  the  absorption  measurements  are  in 
good  agreement  with  the  authors'  calculated  results  obtained  from  the  kinetic  model 
used  by  Westenberg  and  Favin84. 

4.4.2  Nonequi librium  Sozzle  Flows  of  Air 

The  chemically-active  neutral  species  of  importance  in  dissociated  air  at  high 
temperatures  are  usually  taken  to  be  N2  ,  02  ,  N  ,  0  ,  and  NO  (Section  2).  This 
composition  neglects  more  complex  molecules,  such  as  high  oxides  of  nitrogen,  which 
are  usually  unimportant  at  high  temperatures  because  of  their  small  concentrations. 

As  discussed  in  Section  3,  an  important  kinetic  role  is  played  by  NO  ,  although  it 
may  only  be  present  in  relatively  small  amounts,  through  the  fast  bimolecular  exchange 
or  "shuffle"  reaction  system 

1.  0  +  N2  NO  +  N 

2.  N  +  02  NO  +  0  . 

These  reactions,  in  addition  to  the  dissociation-recombination  reactions 

3.  02  +  M  <-*  20  +  M 

4.  N2  +  M  2N  +  M 

5.  N0  +  M?±N+0  +  M 

constitute  the  important  elementary  reactions  governing  the  neutral  species  in  high- 
temperature  air,  as  was  first  demonstrated  in  the  normal  shock  wave  calculations  of 
Duff  pnd  Davidson 12 6 ,  and  Lin  and  Teare127. 

Numerical  solutions  for  nonequilibrium  nozzle  flows  of  dissociated  air  based  on  the 
coupled  rate  equations  for  the  above  kinetic  model  were  first  obtained  by  Eschenroeder, 
Boyer,  and  Hall97,98  for  a  range  of  conditions  of  interest  in  hypersonic  shock-tunnel 
testing.  These  solutions  included  equilibrium  reservoir  temperatures  up  to  8000°K 
at  reservoir  pressures  from  10  to  1000  atm;  the  calculations  were  subsequently  extended 
by  Lordi  and  Mates100,  101  to  15,000°K  and  9000  atm.  These  studies  were  essentially 
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limited  to  determination  of  neutral  species  concentrations  and  the  flow  variables. 

For  the  conditions  considered,  the  energy  involved  in  ionization  was  sufficiently 
small  that  the  ionization  kinetics  could  be  uncoupled  from  the  fluid  dynamic  and 
neutral-species  kinetics  equations.  Corresponding  nonequilibrium  solutions  for  the 
ionized  species  were  subsequently  obtained  by  Eschenroeder61,  10?. 

The  nozzle-flow  solutions  for  the  above  kinetics  reveal  an  important  role  of  the 
bimolecular  shuffle  reactions  1  and  2  in  delaying  N-atom  freezing.  On  flow  expansion 
from  the  nozzle  reservoir  state,  reaction  1,  and  also  reaction  2  at  higher  temperatures 
and  lower  pressures,  depart  from  equilibrium  in  the  direction  of  N-atom  removal  because 
of  the  corresponding  low  activation  energies.  Even  at  relatively  low  concentrations 
of  0?  and  NO  ,  these  bimolecular  paths  are  much  more  efficient  for  recombining  N 
atoms  during  the  flow  expansion  than  the  three-body  recombination  path  via  reaction  4. 

As  a  consequence,  the  freezing  of  N  atoms  can  be  greatly  delayed  beyond  the  point 
where  three-body  recombination  from  reaction  4  has  essentially  ceased.  Practically 
all  N  is  thereby  recombined  to  N2  by  this  bimolecular  path  in  dissociated  air 
expansions  from  intermediate  (100  atm)  or  higher  reservoir  pressures.  This  is  in 
contrast  to  the  very  rapid  freezing  of  N  which  would  occur  if  three-body  recombination 
through  reaction  4  were  the  only  path  for  N-atom  removal. 

The  foregoing  remarks  are  illustrated  by  Figure  4-8  from  Reference  97  which  shows 
(calculated)  equilibrium  and  nonequilibrium  species  mass  concentrations  versus  nozzle 
expansion  ratio  for  air  initially  at  8000°K  and  100  atm  pressure  expanding  in  a  nozzle 
of  geometry  stated  in  Section  4.2  with  L/a  =  1  cm  .  At  these  reservoir  conditions 
the  nitrogen  is  appreciably  dissociated,  and  the  oxygen  essentially  completely  so. 

The  initial  energy  in  nitrogen  dissociation  actually  exceeds  that  in  oxygen  dissociation, 
being  about  25%  of  the  total  enthalpy.  The  mass  concentration  of  0  atoms  is  seen 
to  “freeze”  essentially  at  the  nozzle  throat.  The  amount  of  0  actually  increases 
somewhat  with  flow  expansion  due  to  O-atom  production  via  the  bimolecular  reactions  1 
and  2  as  three-body  recombination  of  0  via  reaction  3  ceases.  By  contrast,  the  con¬ 
centration  of  N  drops  rapidly  with  flow  expansion  due  to  the  bimolecular  path  for 
N-atom  removal.  When  N-atom  freezing  finally  sets  in,  the  N-atom  concentration  is 
only  about  1%  of  the  reservoir  value. 

Corresponding  calculated  results  from  Reference  97  for  reservoir  pressures  of  10 
atm  and  1000  atm  at  8000°K  are  shown  in  Figures  4-9  and  4-10,  respectively.  At  low 
pressures  (Fig.  4-9),  the  initial  N-atom  concentration  is  very  large.  Because  of  the 
low-density  level  not  even  the  fast  bimolecular  reactions  are  effective  in  recombining 
N  atoms  in  the  flow  expansion  in  this  case.  Hence,  freezing  of  both  0  and  N 
occurs  at  near  reservoir  levels.  At  high  pressures  (Fig.  4-10)  the  N-atom  concentration 
is  initially  small  and  energetically  unimportant,  and  is  equilibrated  during  the  flow 
expansion.  Three-body  recombination  is  greatly  enhanced  at  this  high  pressure,  and, 
as  a  result,  0-atom  freezing  is  appreciably  delayed.  Similar  trends  with  increase  in 
reservoir  pressure  are  shown  in  the  calculations  of  References  100  and  101  for  higher 
reservoir  temperatures. 

It  may  be  noted  that  the  calculated  entropy  changes  occurring  in  these  nonequilibrium 
air  expansions  are  typically  a  very  small  fraction  of  the  reservoir  entropy;  i.e.  the 
expansions  are  closely  isentropic.  The  use  of  the  reservoir  entropy  to  correlate 
certain  frozen  properties  of  diatomic  gas  expansions  for  given  nozzle  geometry  was 
discussed  in  Section  4.2.  Lordi  and  Mates100,101  confirmed  that  a  similar  correlation 
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applied  to  the  exact  solutions  obtained  for  air  expansions.  Figure  4-11  shows  such  a 
correlation  of  the  calculated  enthalpy  frozen  in  the  expansion  (hF)  for  the  range  of 
air  reservoir  conditions  noted  and  the  nozzle  geometry  previously  stated  (L/a  =  1  cm). 

The  rather  marked  change  in  slope  at  about  S  =  2.7  cal/gm  °K  is  due  to  the  onset  of 
significant  freezing  of  N  atoms.  At  higher  values  of  S  ,  oxygen  is  essentially  all 
dissociated  and  frozen.  Other  nonequilibrium  effects  as  frozen  species  mass  fractions 
and  reductions  in  pressure  and  temperature  are  also  correlated  by  the  entropy100,  101'104. 

Experimental  studies  of  nozzle  flows  of  dissociated  air  wherein  neutral  species 
concentrations  were  determined  appear  to  be  virtually  nonexistent  to  date.  However, 
a  number  of  experimenters  have  measured  gasdynamic  quantities,  such  as  pressure,  in 
dissociated  nozzle  airflows  and  demonstrated  the  expected  magnitude  of  departures 
from  calculated  equilibrium  values.  Measurements  of  axial  static  pressure  distributions 
have  been  reported,  for  example,  by  Nagamatsu  and  co-workers33,  105,  by  Duffy106,  and 
by  Zonars107.  In  addition  to  static  pressure  distributions  along  the  nozzle  wall, 

Zonars  also  measured  the  transverse  distributions  of  pitot  pressure,  stagnation-point 
heat  transfer,  and  local  mass  flow  at  the  nozzle  exit.  Measurements  of  local  molecular 
speed  ratio  (proportional  to  Ui  (u/T)  .  where  U  is  local  free-stream  velocity,  fJ.  is 
molecular  weight,  and  T  is  translational  temperature)  have  been  made  by  Vidal,  Skinner, 
and  Bartz108  in  shock-tunnel  nozzle  expansions  of  dissociated  air  using  small-scale 
probe  techniques  which  depend  on  free-mole  :ule  heat  transfer. 

In  general,  these  experimental  studies  typically  show  a  fairly  good  agreement  between 
measured  quantities  and  nonequilibrium  values  calculated  on  the  basis  of  the  chemical 
kinetic  model  for  air  already  discussed.  Figure  4-12,  for  example,  shows  a  comparison 
given  by  Duffy106  of  measured  and  predicted  axial  static  pressure  distributions  in  a 
30°  total  angle  conical  nozzle  for  air  reservoir  conditions  of  7000°K  and  250  lb/in2 
absolute.  The  experimental  pressure  data  were  carefully  corrected  for  boundary -layer 
and  static-pressure  probe  size  effects.  The  reduction  of  the  measured  pressures  below 
equilibrium  values  is  seen  to  be  very  large  at  area  ratios  exceeding  100.  The  theore¬ 
tical  pressure  distribution  shown  for  nonequilibrium  flow,  calculated  by  Lordi106  for 
the  chemical  kinetic  model  discussed  above,  is  in  close  agreement  with  the  experimental 
data.  It  may  be  noted  that  for  the  conditions  of  this  experiment,  the  bimolecular 
exchange  reactions  previously  discussed  play  an  important  role  in  recombining  N 
atoms  and  thus  determining  the  static  pressure.  Approximately  equal  amounts  of  energy 
are  involved  initially,  in  the  reservoir,  in  oxygen  and  nitrogen  dissociation.  Whereas 
the  energy  in  nitrogen  dissociation  is  substantially  recovered  during  the  expansion, 
that  in  oxygen  dissociation  is  largely  lost  due  to  the  freezing  of  0  atoms  close  to 
the  nozzle  throat. 

While  these  experiments  do  provide  at  least  some  confidence  as  regards  the  prediction 
of  gasdynamic  effects  due  to  nonequilibrium,  there  is  clearly  need  for  a  more  critical 
assessment  of  theoretical  predictions  by  direct  measurement  of  species  concentrations 
in  nozzle  flows  of  dissociated  air.  As  discussed  in  Section  4.2,  static  pressure  and 
translational  temperature  mainly  reflect  the  energy  balance,  and  therefore  tend  to  be 
relatively  insensitive  to  detailed  features  of  the  kinetics  which  do  not  involve  much 
energy. 

The  foregoing  discussion  has  been  concerned  only  with  the  nonequilibrium  behavior 
of  the  neutral  species  in  nozzle  flows  of  dissociated  air.  As  already  mentioned,  for 
the  conditions  considered  the  energy  in  ionization  is  sufficiently  small  that  ionization 
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phenomena  do  not  influence  the  thermo-gasdynamic  behavior  of  the  flow  significantly. 

Thus  ionization  nonequilibrium  phenomena  can  be  considered  separately  on  the  basis  of 
the  thermo-gasdynamic  history  of  the  flow  determined  from  the  neutral  species  only. 
Estimates  of  nonequilibrium  ionization  phenomena  occurring  in  shock-tunnel  nozzle 
airflows  for  the  regime  of  reservoir  conditions  previously  mentioned  (i.e.  reservoir 
temperatures  up  to  8000°K  and  reservoir  pressures  in  the  range  100  to  1000  atm)  were 
first  reported  by  Eschenroeder  and  Daiber65.  These  authors  considered  the  following 
two  elementary  reactions  for  electron  (e‘)  removal,  or  deionization,  in  the  expanding 
flow: 

(1)  NO*  +  e'  -  N  +  0 

(2)  NO*  +  e‘  +  M  -  NO  +  M  (M  is  any  neutral  particle). 

In  most  of  the  range  of  air  reservoir  conditions  considered,  the  dominant  positive 
and  negative  species  are  NO*  and  e"  whose  concentrations  were  assumed  equal.  The 

bimolecular  reaction  (1),  termed  dissociative  recombination,  involves  dissociation  of 

the  NO  molecule  by  the  energy  given  up  in  electron  recombination,  whereas  reaction 
(2)  involves  the  classical  Thomson-type  three-body  recombination  with  a  third  and 
neutral  body  M  required  to  take  up  the  energy  of  recombination.  On  the  assumption 
of  thermochemical  equilibrium  of  the  neutral  species,  the  authors  employed  these  ion 
kinetics  tc  estimate  the  approximate  freezing  behavior  of  free  electrons  in  shock- 
tunnel  air  expansions.  These  estimates  indicated  appreciable  freezing  of  the  electron 
mass  concentration,  and  also  that  the  dissociative-recombination  reaction  (1)  was  the 
dominant  process  in  removing  electrons,  except  possibly  at  quite  high  densities. 

Despite  the  appreciable  freezing  indicated  for  the  electron  fraction,  the  electron 
number  densities  for  nozzle  expansions  to  typical  hypersonic  test  conditions  were 
predicted  to  be  sufficiently  low  to  avoid  significant  attenuation  of  electromagnetic 
propagation  for  microwave  frequencies  above  about  10s*  c/s.  This  consideration  is 
important  in  connection  with  shock  tunnel  applications  involving  microwave  transmission 
through  the  ionized  shock  layer  formed  about  a  test  model. 

Subsequent  studies  by  the  same  authors  included  analysis  of  the  effects  of  seeding 
air  expansions  with  cesium,  in  order  to  increase  the  electron  density61;  machine  solu¬ 
tions  for  nozzle  airflow  ionization  considering  chemical  nonequilibrium  as  well  as  a 
variety  of  elementary  ion  kinetic  mechanisms  in  addition  to  reactions  (1)  and  (2)61,63 
and  experimental  determination  of  electron  densities  and  collision  frequencies  in 
shock-tunnel  nozzle  airflows  by  means  of  microwave  interferometry63.  These  further 
studies  essentially  confirmed  the  dominance  of  reaction  (1)  as  regards  the  most 
important  ion  kinetic  mechanism  for  electron  removal  in  the  regime  of  shock-tunnel 
airflows  considered.  The  more  exact  solutions  also  indicated  that  the  electron  mass 
fraction  does  not  typically  freeze  fully,  but  may  continue  to  decrease  with  flow 
expansion.  The  typical  nonequilibrium  behavior  of  electron  mass  fraction  and  the 
neutral  species  as  predicted  by  machine  calculation  is  shown  in  Figure  4-13.  The  lack 
of  complete  freezing  of  the  electron  fraction  is  a  consequence  of  the  dominant  bi¬ 
molecular  recombination  reaction  and  the  inverse  3/2  power  dependence  of  the  associated 
rate  coefficient  on  temperature.  As  discussed  by  Eschenroeder61,  and  more  recently  by 
Cheng  and  Lee23’24,  the  asymptotic  freezing  behavior  also  depends  on  the  frozen  neutral- 
species  chemistry  and  the  nozzle  geometry. 

In  the  past  few  years,  intensive  effort  has  continued  on  the  study  of  the  complex 
ionization  phenomena  occurring  in  high-temperature  airflows.  Much  of  this  effort  has 
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concerned  Ionization  nonequilibrium  In  external  flows  about  hypervelocity  bodies, 
particularly  the  problea  of  electron  density  decay  In  wake  flows.  Considerable  progress 
has  been  made  In  further  Identification  of  important  elementary  Ion  and  electron  re¬ 
actions  and  determination  of  specific  rate  coefficients,  as  discussed  In  recent  surveys 
of  Dunn12',  and  Carnicom12*.  for  example.  As  regards  ionization  phenomena  In  hypersonic 
nozzle  flows  of  air,  the  more  recent  studies  of  Ionization  kinetics  do  not  sppear  to 
offer  any  basis  for  serious  conflict  with  the  general  picture  Indicated  above  for 
reservoir  conditions  in  the  range  noted.  In  particular,  the  conclusion  remains  that 
the  dissociative  recombination  of  NO*  dominates  the  removal  of  free  electrons  In 
nozzle  airflows  typical  of  current  shock-tunnel  techniques.  However.  It  must  be 
emphasized  that  there  Is  still  a  serious  lack  of  definitive  experimental  investigation 
of  nozzle  flows  of  ionized  air;  the  need  here  is  perhaps  even  greater  than  the  need  for 
data  on  the  neutral  species  chemistry.  One  important  basic  question  that  remains  to 
be  clarified  is  that  of  the  electron  temperature  behavior.  Numerous  experiments  in¬ 
dicate  that  the  electron  temperature  does  not  equilibrate  with  the  translational 
temperature  in  nozzle  expansions,  but  tends  to  remain  relatively  hlgh,,,s2.  Certain 
specific  rate  coefficients  of  elementary  reactions  involving  free  electrons  depend 
strongly  on  this  temperature.  Of  course,  there  is  a  continuing  need  for  data  on 
specific  rate  coefficients.  Various  elementary  ionization  reactions  have  been  con¬ 
sidered  theoretically  for  which  little  or  no  rate  data  have  been  obtained  in  air. 

One  example  is  that  of  electron  recombination  via  a  three-body  collision  process,  with 
a  free  electron  serving  as  a  catalytic  third  body.  On  the  basis  of  present  estimates, 
which  suggest  a  relatively  large  rate  coefficient,  this  process  could  be  of  cowarable 
importance  to  two-body  dissociative  recombination  for  electron  removal  at  shock-tunnel 
reservoir  conditions  giving  very  high  electron  densities. 

4.4 .3  Sonequi librium  Effects  in  Hypersonic  Testing 

The  occurrence  of  such  nonequilibriua  phenomena  as  discussed  in  the  preceding  section 
creates  a  serious  problem  as  regards  model  testing  in  hypervelocity  airflows  obtained 
by  nozzle  expansion  from  high  enthalpy  conditions.  The  problem  is  now  well  recognized, 
if  not  solved.  Current  shock-tunnel  anc.  arc-jet  test  facilities  already  operate  at 
sufficiently  high  reservoir  temperatures  to  encounter  significant  nonequilibrium 
effects  at  lower  pressures;  the  effects  will  become  more  severe  with  future  facilities 
having  higher  enthalpy  performance.  While  it  has  been  suggested  that  nozzle  flow 
nonequilibrium  might  be  favorably  exploited  in  hypersonic  testing  (e.g.  Reference  112 
as  regards  the  increase  in  flow  Mach  number  and  Reynolds  number  due  to  freezing),  the 
customary  view  is  that  such  effects  are  undesirable  and  to  be  eliminated  if  possible. 

01  course,  this  view  is  in  the  vein  of  duplicating  in  the  wind-tunnel  model  test  the 
ambient  air  chemistry  through  which  the  corresponding  flight  vehicle  flies. 

In  recent  years,  considerable  effort  has  been  devoted  to  studying  possible  means 
of  alleviating,  or  at  least  minimizing,  the  nonequilibrium  test-flow  problem  and  also 
to  assessing  the  consequents  for  typical  hypersonic  test  bodies  of  interest.  Although 
most  of  these  studies  have  been  theoretical,  their  predictions  must  be  given  general 
credence  in  view  of  the  improved  status  of  the  knowledge  of  air  kinetics  today,  and 
the  nozzle  flow  experiments  previously  discussed. 

In  essence,  the  studies  reported  to  data  have  not  indicated  any  method  by  which 
nonequi librium  effects  can  be  essentially  eliminated  in  conventional  nozzle  expansion 
flows.  Studies  to  determine  optimum  nozzle  shapes  giving  sufficiently  slow  expansions 
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to  maintain  near-equilibrium  flow,  and  duct  lengths  required  for  equilibration  of 
frozen  flow,  indicate  these  approaches  to  be  of  limited  usefulness  because  the  rather 
extreme  lengths  required  would  entail  excessive  boundary- layer  growth17,  ie.  Likewise, 
within  practical  limitations  the  increase  of  nozzle  scale,  or  size,  for  a  given  shape 
offers  very  limited  benefits  because  of  the  relative  insensitivity  of  chemical  freezing 
to  scale17,  18,  100,  101.  The  possibility  of  seeding  dissociated  air  with  small  amounts 
of  additive  to  introduce  fast  bimolecular  paths  for  0-atom  removal  (similar  to  the 
seeding  of  H2  expansions  previously  mentioned)  has  received  limited  study.  The  cal¬ 
culations  of  Reference  130  for  such  seeding  of  air  expansions  with  0.01  mole  fraction 
of  H2  suggest  that  this  technique  probably  offers  quite  limited  benefits  in  a 
practical  sense.  Thus  it  seems  that  the  most  effective  method  for  suppressing  nozzle 
airflow  nonequilibrium  still  remains  that  of  using  high  reservoir  pressures,  as 
earlier  calculations  of  nonequilibrium  nozzle  airflows  have  indicated12,99,100.  Of 
course,  this  approach  is  compatible  with  the  need  for  high  reservoir  pressures  in 
order  that  the  test  flow  duplicate  ambient  density  conditions  appropriate  to  hyper- 
velnrlty  flight  at  lower  altitudes. 

Before  reviewing  the  general  magnitude  of  nozzle  test-flow  nonequilibrium  effects, 
it  may  be  noted  that  even  if  such  effects  were  completely  eliminated,  and  ambient 
flight  conditions  duplicated,  a  problem  would  still  remain  as  regards  the  proper 
tunnel  simulation  of  nonequilibrium  effects  occurring  in  the  flow  about  the  flight 
vehicle.  For  similar  geometries  of  the  test  model  and  flight  vehicle,  proper  simulation 
requires,  in  general,  that  the  model  be  full  scale.  As  is  well  kn>'wn,  this  is  necessary 
in  order  that  the  various  time  scales  governing  nonequilibrium  phenomena  remain  in  the 
correct  ratio  with  fluid  particle  transit  times  characterized  by  L/U^  ,  where  L  is 
the  body  length,  say,  and  is  ambient  velocity  in  flight  (and  tunnel).  Of  course, 
the  traditional  usefulness  of  wind  tunnels  is  in  testing  models  of  reduced  scale.  For 
ideal  gas  flows,  this  can  be  done  without  sacrificing  correct  similitude,  at  least  in 
principle.  However,  for  simulation  of  flight  situations  where  nonequilibrium  effects  are 
important,  the  constraints  are  rather  severe  if  correct  similitude  is  to  be  attained.  Of 
course,  other  than  full-scale  testing,  success  here  depends  on  sufficient  a  priori  under¬ 
standing  of  the  flight  problem  that  limited  similitudes  or  partial  simulation  techniques 
can  be  exploited. 

As  a  final  remark  in  this  digression,  it  may  be  mentioned  that  limited  similitudes 
have  been  established  which  permit  this  requirement  of  a  full-scale  model  to  be 
relaxed  for  certain  nonequilibrium  flow  situations.  Obviously,  the  two  limiting 
situations  of  either  completely  equilibrium  or  completely  frozen  flow  about  the  flight 
vehicle  permit  model  scale  to  be  es;  entially  arbitrary.  This  requires,  of  course, 
that  the  same  limiting  behavior  obtains  in  the  model  flow  field  (which  is  generally 
difficult  to  accomplish  for  the  equilibrium  case)  and  that  the  appropriate  fluid 
dynamic  similitude  is  maintained  (i.e.  M^,  Re,*,,  etc.).  An  example  here,  which  also 
offers  flexibility  as  regards  relaxing  duplication  of  ,  is  the  equilibrium-flow 
inviscid  similitude  for  blunt-noze  1  slender  bodies  due  to  Cheng  (Reference  90,  Appendix 
B).  This  similitude,  which  entails  the  assumptions  of  hypersonic  small-disturbance 
theory,  was  also  extended  to  include  nonequilibrium  flow.  In  the  nonequilibrium  case, 
still  need  not  be  duplicated,  but  the  scale  L  is  no  longer  arbitrary  but  rather 
fixed  according  to  the  choice  of  in  order  to  preserve  the  flow  transit  time 

L/Uo,  .  Inger15  has  given  the  complete  formalism  for  this  type  of  similitude.  A  second 
special  similitude  which  may  be  noted  for  nonequilibrium  flows  is  that  usually  referred 
to  as  binary  scaling.  This  entails  an  assumption  about  the  nature  of  the  nonequilibrium 
rate  processes  rather  than  about  the  fluid  dynamics,  namely,  that  all  the  rate  processes 
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involve  only  two-body  or  binary  collisions.  Thus  three-body  recombination  is  excluded. 
In  this  case,  a  simplifying  similitude  results  in  terms  of  the  product  pL  ,  where  p 
is  a  characteristic  density  level  in  the  flow  field  and  L  is  the  body  scale.  Detailed 
considerations  and  applications  of  binary  scaling  are  given  in  References  111,  118,  and 
References  37-40  and  53  or  Appendix  B. 

Returning  to  the  question  of  nozzle-flow  nonequilibrium,  it  is  instructive  to  review 
the  magnitude  of  predicted  effects  in  the  expanded  test  gas  for  a  range  of  reservoir 
conditions  representative  of  the  limits  of  current  shock-tunnel  testing  techniques. 

Table  I  shows  machine  calculation  results  from  Reference  111,  for  complete  air  kinetics, 
for  the  mass  fractions  of  0  ,  N  ,  and  NO  frozen  in  nozzle  air  expansions  from  6000 
and  8000°K  reservoir  temperature  and  various  reservoir  pressures.  The  corresponding 
frozen  enthalpy  is  also  given.  These  results,  which  apply  for  the  nozzle  geometry 
A/A^  -  l  +  (Lx/a)2  with  L/a  -  1  cm  ,  illustrate  the  need  for  high  reservoir  pressures 
to  suppress  nonequilibrium  effects.  This  is  particularly  the  case  as  regards  oxygen- 
atom  recombination  which  depends  on  three-body  collisions.  The  effectiveness  of  the 
bimolecular  reactions  discussed  for  recombining  nitrogen  atoms  is  very  evident  in  the 
results  for  8000°K  and  100  atm  where  essentially  all  N  atoms  are  recombined.  At 
the  extreme  of  high  reservoir  temperatures  and  low  reservoir  pressures,  the  gasdynamic 
effects  of  freezing  are  very  large;  about  50£  of  the  total  enthalpy  is  frozen  in  the 
dissociated  species  at  8000°K  and  10  atm.  Whereas  lower  reservoir  temperatures  and 
higher  reservoir  pressures  reduce  the  frozen  mass  fractions  of  0  and  N  atoms,  and 
also  reduce  the  frozen  enthalpy  fraction,  it  will  be  noted  in  Table  I  that  this  trend 
is  reversed  for  the  frozen  fraction  of  NO  .  Thus  at  low  temperatures  and  high  pres¬ 
sures,  NO  tends  to  become  the  dominant  frozen  "non-equilibrium”  species,  as  illus¬ 
trated  in  Table  I  for  G000°K  and  1000  atm.  At  the  extreme  of  low  temperatures  and 
high  pressures  in  this  respect,  higher  oxides  of  nitrogen  not  included  in  the  present 
calculations  might  also  be  frozen  at  levels  of  some  importance. 

For  expansions  giving  equilibrium  airflow  Mach  numbers  of  20,  Table  II  shows  results 
corresponding  to  Table  I  for  the  calculated  ratios  of  actual  (i.e.  nonequilibrium)  to 
equilibrium  values  of  the  test-flow  gasdynamic  quantities  of  interest.  The  Mach 
number  in  Table  II  is  u/[(dp/dx)/(dp/dx)] «  ,  and  the  ratios  of  stream  Reynolds  numbers 
are  based  on  a  linear  viscosity-temperature  relation.  The  tabulated  ratios  are  rather 
weakly  dependent  on  the  flow  or  nozzle  expansion  ratio  in  each  case,  so  that  the 
results  are  representative  of  an  appreciable  range  of  equilibrium  Mach  number.  As 
expected  from  pure  diatomic  gas  examples,  the  effects  of  freezing  are  largest  on  the 
translational  temperature  and  pressure,  and  much  smaller  on  the  density 'and  velocity. 

The  reductions  in  temperature  and  pressure  below  corresponding  equilibrium  values 
become  very  large  at  8000°K  and  100  atm  reservoir  conditions.  The  reduction  in  tem¬ 
perature  strongly  increases  the  Mach  number  and  more  strongly  the  stream  Reynolds 
number  in  this  case*. 

The  generally  large  magnitude  of  predicted  nonequilibrium  effects  in  the  test  gas 
for  expansions  from  high  enthalpies  emphasizes  the  need  for  much  more  complete  measure¬ 
ment  of  test  flow  characteristics  than  has  traditionally  been  necessary  in  wind-tunnel 
testing.  Clearly  the  need  here  exceeds  current  state-of-the-art  diagnostic  techniques 
as  regards  routine  measurement.  This  is  most  obvious  for  the  determination  of  the 
frozen  chemical  composition  of  the  test  gas,  but  also  applies  to  determination  of  the 

•  Results  similar  to  Tables  I  and  II,  but  for  reservoir  temperatures  to  15,000°K,  are  given  by 
Lordi  and  Mates100. 
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static  temperature  and  pressure  of  a  hypersonic  test  flow.  The  customary  measurement 
of  pitot  pressure  is  of  very  limited  usefulness  as  an  indication  of  nonequilibrium 
effects  since  this  quantity  is  essentially  and  therefore  is  quite  insensitive 

to  nonequilibrium  in  hypersonic  nozzle  flows. 

A  number  of  studies  have  now  been  made,  mostly  theoretical,  of  the  effects  of  test- 
flow  nonequilibrium  on  the  flow  fields  about  typical  hypersonic  test  bodies109"118. 
Typically  a  comparison  is  made  of  model  flow-field  quantities  for  equilibrium  and  non¬ 
equilibrium  test  streams  expanded  to  the  same  area  ratio  from  the  same  reservoir 
state.  For  equilibrium  flow  fields  behind  strong  shock  waves,  the  effects  of  test- 
flow  nonequilibrium  tend  to  be  quite  limited.  An  important  example  of  this  is  the 
equilibrium  stagnation- region  flow  on  a  blunt  body.  The  reduction  in  pitot  pressure 
p'  due  to  nozzle  freezing  is  roughly  of  the  order 

Pp  ~  P0  EQ  ^  Uoo  ~  Ucd  EQ  ^  1  hF 

P0  EQ  EQ  ^  hQ 

for  hypersonic  frozen  test  flows,  where  subscript  EQ  denotes  values  for  the  corres¬ 
ponding  equilibrium  test  stream.  Thus  p'  cannot  be  very  strongly  reduced  (see 
Table  I  for  hp/h0),  and  since  hQ  is  unchanged,  the  entire  equilibrium  state  determined 
by  hQ  and  p'  must  remain  relatively  insensitive.  In  particular,  the  gas  temperature 
and  composition  tend  to  be  quite  insensitive,  while  the  density  is  reduced  roughly 
the  same  order  as  the  pressure.  Local  equilibration  behind  a  strong  shock  thus  tends 
to  erase  the  nonequilibrium  history  of  the  nozzle  expansion  process. 

For  nonequilibrium  flow  fields  behind  strong  shocks,  the  local  flow  field  chemistry 
can,  of  course,  be  strongly  affected  by  test-flow  nonequilibrium.  Also,  effects  can 
be  appreciable  on  local  temperature  and  density.  However,  the  pressure  will  tend  to 
remain  rather  insensitive,  particularly  for  blunt  bodies  where  the  flow  is  Newtonian 
in  character.  As  an  example,  figures  4-14  and  4-15  from  Reference  111  show  calculated 
local  temperatures  and  species  concentrations  in  the  nonequilibrium  flow  field  of  an 
axisymmetric  blunt  nose  for  test-flow  conditions  attained  by  equilibrium  and  non¬ 
equilibrium  expansions  of  air  from  a  reservoir  state  of  6000°K  and  100  atm.  The  air 
kinetic  models  for  this  example  included  the  coupled  reactions  previously  discussed. 
Figure  4-14  shows  the  temperature  level  of  most  of  the  nonequilibrium  shock  layer  to 
be  considerably  reduced  by  the  test  stream  nonequilibrium  due  to  the  energy  initially 
frozen  in  dissociation.  The  shock  layer  density  level  is  also  reduced,  as  evidenced 
by  the  increase  in  shock  standoff  distance.  Figure  4-15  shows  the  very  large  effect 
of  nozzle  freezing  in  greatly  increasing  the  concentrations  of  atomic  oxygen  and 
nitric  oxide  along  a  ray  at  7°  to  the  body  axis. 

In  the  limiting  situation  where  the  chemistry  is  frozen  throughout  the  model  flow 
field,  as  tends  typically  to  be  the  case  with  slender  models  of  small  scale,  then  the 
flow  field  is  that  of  an  ideal  gas  with  a  (frozen)  specific  heat  ratio  determined  by 
the  free-stream  composition.  For  moderate  degrees  of  nozzle  freezing,  the  frozen 
specific  heat  ratio  will  be  rather  close  to  that  of  undissociated  air.  This  is  parti¬ 
cularly  true  when  diatomic  NO  dominates  the  frozen  "nonequilibrium”  species.  Such 
frozen-flow  model  tests  can  provide  useful  data  provided  the  free-stream  properties 
are  known.  Of  course,  the  key  assumption  involved,  that  of  frozen  flow  about  the  test 
model,  is  usually  difficult,  if  not  impossible,  to  verify  by  direct  measurement. 
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Usually  theoretical  estimates  of  the  inviscid  and  boundary-layer  flow  behavior  must 
be  relied  on  in  this  respect.  The  possibility  of  surface  catalytic  effects  arises  in 
connection  with  surface  heat  transfer  measurements  in  frozen  boundary- layer  flows 
(Refs.  35,  131,  132). 

As  a  final  remark,  it  should  perhaps  be  emphasized  that  current  hypersonic  test 
facilities  are  being  used  for  research  studies  in  a  performance  range  where  appreciable 
nozzle  freezing  does  occur.  For  example,  recent  shock-tunnel  applications  at  the 
Cornell  Aeronautical  Laboratory  include  studies  of  hypersoiic  airflows  over  blunt  and 
slender  bodies  at  conditions  where  nonequilibrium  effects  are  significant133'136. 

These  applications  have  included  extension  of  shock-tunnel  techniques  to  measurements 
of  electromagnetic-microwave  transmission  through  the  ionized  shock  layer  on  the 
model,  including  determination  of  radiation  patterns  from  antennae  located  in  the 
model133,  13u.  In  conjunction  with  such  studies,  methods  of  accounting  for  the  effects 
of  nozzle  freezing  are  being  developed,  and  this  trend  will  undoubtedly  continue.  A 
promising  approach  to  more  general  simulation  relationships  appropriate  to  the  problem 
has  already  been  made  by  Gibson116.  It  is  clear  that  the  more  sophisticated  applica¬ 
tions  of  wind-tunnel  testing  methods  in  this  respect  will  require  much  more  complete 
measurement  of  local  flow  properties,  as  previously  noted. 
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TABLE  I 

Percentage  of  0,  N,  NO,  and  Reservoir  Enthalpy 
Frozen  in  Hypersonic  Nozzle  Airflow 

l  -  1  cm 


T0  “K  - - 

6000 

8000 

Pq  atm  - ► 

100 

1000 

10 

100 

1000 

FROZEN  MASS  7  OF  0 

9.5/. 

1.127 

23.2  7 

22.5  7 

6.207 

FROZEN  MASS  %  OF  N 

0 

0 

30.2 

0.120 

0 

FROZEN  MASS  %  OF  NO 

6.75 

9.36 

0.033 

0.771 

7.29 

FROZEN  ENTHALPY  100  hp/h0 

15.5 

4.80 

50.8 

8.13 

TABLE  II 

Ratios  of  Actual  to  Equilibrium  Airflow  Variables  for 
Coupled-Reaction  Expansions  Giving  Equilibrium  Mach  Number  20 

l  -  1  cm 


To 

6000  °K 

8000  °K 

po 

100  atm 

1000  atm 

100  atm 

1000  atm 

p/p  EQ 

0.425 

0.875 

0.159 

0.544 

t'te« 

0.359 

0.855 

0.121 

0.497 

^EQ 

1.095 

1.014 

1.110 

1.045 

0.913 

0.987 

0.901 

0.957 

m/mE)? 

1.455 

1.06 

2.30 

1.305 

Re  /Rey 
x  xE(J 

2.79 

_ 

1.17 

8.25 

2.01 

nnoti  oiiMCKt 


Numerical  solution  for  translational  temperature  in  nozzle  flow  of 
dissociated  oxygen  (from  Reference  18) 


4-3  Numerical  solution  for  pressure  in  nozzle  flow  of  dissociated  oxygen 

(from  Reference  18) 


5.0 


Numerical  solution  for  velocity  in  nozzle  flow  of  dissociated  oxygen 

(from  Reference  18) 


counts  it  act  too  (limes 
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Fig.  4-12  Comparison  of  measured  and  predicted  pressure  Fig. 4-13  Numerical  solution  for  ionized  and  neutral 

distributions  for  nozzle  flow  of  air.  T„  =  7000°K,  species  concentrations  in  nozzle  flow  of  air. 

P0  =  250  lb/in2  absolute,  l  =  2.75  cm  T0  =  5800°K,  PQ  =  657  atm  (from  Reference  134) 

(from  Reference  106) 
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Pig. 4- 14  Shock  location  and  temperature  distribution  in  Fig. 4-15  Concentrations  of  0  and  NO  in  none qui librium 
nonequilibrium  blunt-nosed  flowfields  for  nonequilibrium  blunt-nose  flowfields  for  nonequilibrium  and  equilibrium 

and  equilibrium  nozzle  airflows.  T0  =  6000°K,  nozzle  airflows.  T0  =  6000°K,  PQ  =  100  atm.  A/A*  =  18. 190, 

P0  =  100  atm.  A/ A*  =  18,190,  Rg  =  0.493  ft  Rg  =  0.493  ft  (from  Reference  111) 

(from  Reference  111) 
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SECTION  5 
CONCLUSION 


The  experimental  and  theoretical  research  in  the  field  of  nonequilibrium  expanding 
flows  has  been  reviewed.  Recent  work  concerning  the  basic  physical-chemical  processes 
involved  are  also  summarized  and  separate  bibliographies  for  various  aspects  of  shock- 
heated  flows  are  included.  There  is  a  strong  emphasis  throughout  on  air  flows,  although 
some  discussion  of  other  gas  mixtures  is  included. 

The  basic  equilibrium  equations  are  summarized  in  Section  2  and  several  tables  of 
thermodynamics  properties  of  gases  are  described.  The  capabilities  of  available 
computer  programs  for  the  solution  of  the  equilibrium  equations  are  discussed.  It  is 
pointed  out  that  the  computer  calculations  are  now  so  rapid  as  to  provide  efficient 
sub-routines  for  other  aerodynamic-flow  calculations. 

Current  problems  associated  with  the  description  of  nonequilibrium  gases  are  reviewed 
in  Section  3.  The  methods  of  treatment  of  vibrational  relaxation  is  emphasized,  and  the 
work  associated  with  translational  and  rotational  equilibration  is  omitted  entirely. 

The  several  methods  available  for  the  calculation  of  vibrational  transition  probabi¬ 
lities  are  compared,  and  the  limitations  of  these  calculations  are  discussed.  Methods 
of  treatment  of  the  master  equations  for  vibrational  relaxation  are  then  described, 
both  for  cases  where  exchange  effects  are  unimportant  and  for  cases  of  pure  molecular 
gases  and  mixtures  of  gases  where  exchange  effects  are  dominant.  It  is  pointed  out 
that  the  non-Boltzmann  effects  that  can  be  obtained  when  vibrational  energy  exchange  is 
important  can  greatly  increase  the  vibrational  relaxation  rates.  For  electronic 
excitation,  the  limited  information  available  for  rates  and  mechanisms  in  the  high- 
temperature  range  is  summarized.  The  basic  problems  in  understanding  simple  chemical 
dissociation  processes  are  delineated,  and  the  models  that  have  been  proposed  for 
treatment  of  vibration -dissociation  coupling  are  reviewed.  Available  experimental 
data  are  discussed  for  dissociation  reactions  and  for  the  nitric  oxide  shuffle  reactions. 
Recent  reviews  of  chemical  problems  associated  with  reentry  are  referenced,  as  are 
reviews  of  relevant  ionization  and  deionization  rates. 

Supersonic-nozzle  flows  are  reviewed  in  Section  4.  First,  the  pseudo-one-dimensional 
formulation  is  described  briefly,  and  the  particular  difficulties  associated  with 
numerical  nonequilibrium  calculations  of  such  flows  are  reviewed.  The  various  methods 
that  have  been  used  to  solve  for  the  steady-state  flow  of  a  reacting,  diatomic  gas  in 
a  nozzle  are  then  described.  The  relevant  parameters  for  the  separation  of  the  non¬ 
equilibrium  phenomena  in  the  flow  field  are  defined,  and  exact  and  approximate  solutions 
are  compared.  The  problem  of  vibrational  relaxation  in  expanding  flows,  without  disso¬ 
ciation,  is  treated  separately.  Recent  data  demonstrating  the  rapid  relaxation  of 
diatomic  gases  in  such  flows  are  reviewed  and  compared  with  theory. 
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The  use  of  machine  programs  to  treat  more  complicated  chemical  flows  is  next  dis¬ 
cussed  in  terms  of  gas  mixtures  of  interest  in  propulsion.  Graphical  results  of  these 
calculations  for  different  test  conditions  are  compared  and  experiments  are  discussed. 
Similar  calculations  and  experiments  are  then  described  for  air,  including  the  effects 
of  seeding  to  increase  the  electron  concentration.  There  is  clearly  a  lack  of  de¬ 
finitive  chemical  experiments  for  expanding  flows  in  air,  and  the  need  for  measurements 
of  chemical  composition  in  such  flows  is  demonstrated.  Finally,  the  results  of  the 
calculations  are  reviewed  in  relation  to  the  effect  of  nonequilibrium  on  hypeiocnic 
testing,  and  the  complications  that  this  introduces  in  aerodynamic  scaling  are  discussed. 
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